Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,263,501 G→C 7.6% P54A (CCG→GCG) ‡ tdcD ← propionate kinase/acetate kinase C, anaerobic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,263,5010GC 7.6% 379.1 / 18.6 119P54A (CCG→GCG) ‡tdcDpropionate kinase/acetate kinase C, anaerobic
Reads supporting (aligned to +/- strand):  ref base G (49/61);  new base C (7/2);  total (56/63)
Fisher's exact test for biased strand distribution p-value = 8.17e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.46e-01

Only 100 of 119 total aligned reads displayed.

TCATCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACC  >  NC_000913/3263353‑3263646
                                                                                                                                                    |                                                                                                                                                    
tCATCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGt                                                                                                                                                     <  1:1336306/150‑1 (MQ=255)
 cATCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTg                                                                                                                                                    <  2:161864/150‑1 (MQ=255)
  agcGGAAATAATGGCTGACTCTGTAAAAATACTGCCGCCGTGATCGAGTCGGTGGCCAATTTAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGc                                                                                                                                                   <  2:301753/148‑1 (MQ=255)
    cGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTg                                                                                                                                                 >  1:646632/1‑150 (MQ=255)
     ggTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTgg                                                                                                                                                >  1:1121505/1‑150 (MQ=255)
      gTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGCTGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCCATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGc                                                                                                                                               <  1:1485690/150‑1 (MQ=255)
      gTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGc                                                                                                                                               <  1:16134/150‑1 (MQ=255)
      gTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGc                                                                                                                                               <  2:646632/150‑1 (MQ=255)
      gTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGGGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGc                                                                                                                                               <  1:529403/150‑1 (MQ=255)
          taatGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGcggtc                                                                                                                                           >  1:134618/1‑146 (MQ=255)
           aatGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑Tcc                                                                                                                                         <  1:1549465/150‑1 (MQ=255)
           aatGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCACCGCTGGCTC‑Acc                                                                                                                                         <  1:1073929/150‑1 (MQ=255)
            atGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑Tccc                                                                                                                                        <  1:1435953/150‑1 (MQ=255)
             tGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCTAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAACGGTGCTGGatcaaccc                                                                                                                                        >  2:1456856/1‑145 (MQ=255)
                  gACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCATTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGCGCCAGCCGTTGTGGGAC‑TGCCCC‑Attt                                                                                                                                 >  2:1586679/1‑150 (MQ=255)
                    cTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTAc                                                                                                                               <  1:212578/150‑1 (MQ=255)
                    cTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCACCCGTGCTGGCTC‑TCCCCC‑ATTTAc                                                                                                                               >  2:506301/1‑150 (MQ=255)
                      cGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACgg                                                                                                                             <  2:1552305/150‑1 (MQ=255)
                        gTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGATCCCGCGGTGCTGGAAC‑TCCCCCAATTTACgtt                                                                                                                            >  2:568247/1‑148 (MQ=255)
                           aaaaTACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATACCGGATTTCCAGTTCAAATGCAATAGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCAAGCGGTGCTGGCTC‑ACCCAC‑ATTTACGGATAAg                                                                                                                        >  1:54956/1‑150 (MQ=255)
                            aaaTACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGa                                                                                                                       >  2:552565/1‑150 (MQ=255)
                            aaaTAATGCCGCCGTGAGCGATGCGTGGGCCAATTAGGGCCACCATGGCTTTTAAATGCGTTTTTGCCATTTAAATTGAAATTGCATTCAATGCCCCTTCTTAGCGTTGGTGAGCAAGCGGTCATGGATC‑TCCCCC‑ATTTACGGATacta                                                                                                                       >  2:1216430/1‑147 (MQ=255)
                                aCTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGGTGGGTC‑GCCCCC‑AATTACGGATAAGAATgg                                                                                                                   >  2:374799/1‑149 (MQ=255)
                                    ccgccgTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCCCCGGGGGGGGGGG‑TCCCCC‑ATTTTAGGATTAGAATGGAt‑tt                                                                                                              >  1:1063222/1‑150 (MQ=255)
                                    ccgccgTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCGA‑TCCCCC‑CATTAACGTTAATAATGGAt‑tt                                                                                                              >  2:597037/1‑150 (MQ=255)
                                     cgccgTGCGCGATGCGGTGGCCAATTAAGGCCACACTGTCAATTAAATTCCGTTGTGCCAGTTCAACTGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑GCCCCC‑ATTTACGGATAAGAATGCAt‑ttt                                                                                                             <  1:287197/150‑1 (MQ=255)
                                     cgccgTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAt‑ttt                                                                                                             >  1:665156/1‑150 (MQ=255)
                                      gccgTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTc                                                                                                            <  1:990538/150‑1 (MQ=255)
                                       ccgTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTcc                                                                                                           >  2:57764/1‑150 (MQ=255)
                                         gTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGCGCCAGGGGTGCTGGGGC‑TCCCCC‑AATTACGGGTTAGAATGCAT‑TTTGCgg                                                                                                         >  1:414324/1‑149 (MQ=255)
                                         gTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGa                                                                                                         >  2:22533/1‑150 (MQ=255)
                                         gTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGa                                                                                                         <  2:162955/150‑1 (MQ=255)
                                            aGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCTTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGtt                                                                                                      <  1:1453424/150‑1 (MQ=255)
                                            aGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGtt                                                                                                      <  1:1111166/150‑1 (MQ=255)
                                                   cGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCg                                                                                               <  1:1019000/150‑1 (MQ=255)
                                                    ggTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCCGCGGGGGTGGCTC‑TACCCC‑ATTTTCGGATTATAATGCAT‑TTTACGAGTTAATAACgg                                                                                              >  1:1703449/1‑149 (MQ=255)
                                                      tGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCg                                                                                            >  1:1015107/1‑150 (MQ=255)
                                                      gggCCAAGTCAGGCCTCACTGTCATTTAAATTCCGTTTTTCCAGTTCAACGGCAATTGCCTTCAATGCACCATCGTAGCTGTGGTGCGCCAGCGGTGCTGGCTC‑TCCGCC‑ATTTACGCATAAGAATGCAT‑TTTCCGAGTTAATACCGGCg                                                                                            <  2:314684/149‑1 (MQ=255)
                                                       ggCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCgg                                                                                           >  1:920587/1‑150 (MQ=255)
                                                        gCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGc                                                                                          <  2:195188/150‑1 (MQ=255)
                                                        gCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGc                                                                                          <  1:825709/150‑1 (MQ=255)
                                                         ccAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑GTCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGCGTTAATACCGTCGGcc                                                                                         >  1:37240/1‑149 (MQ=255)
                                                           aaTTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAAt                                                                                       >  1:36153/1‑150 (MQ=255)
                                                                aGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTg                                                                                  <  2:842179/150‑1 (MQ=255)
                                                                 ggCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGa                                                                                 <  1:1448372/150‑1 (MQ=255)
                                                                  gCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTTGAGGCTCAAACCCC‑ATTTAC‑GATTAGAAATCAA‑TTTACGAGATTATACCCTAGGGAAAACCTGaa                                                                                >  2:1477218/1‑149 (MQ=255)
                                                                  gCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGAc                                                                                <  2:161830/150‑1 (MQ=255)
                                                                  gCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGAc                                                                                <  2:1453338/150‑1 (MQ=255)
                                                                  gCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGAc                                                                                <  1:606572/150‑1 (MQ=255)
                                                                     acacTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCTGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAAAGCAT‑TTTCCGAGTTAATACCGTCGGCAATTCCTGACAtt                                                                             >  1:242259/1‑150 (MQ=255)
                                                                      cacTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGCGCCACCGGTGGGTGCGG‑ACCCCC‑AGTTAAAGAAAAGAATGCATGTTAACAGGTTAAAACCG‑GGGGAATACCTGACagtt                                                                            >  2:871129/1‑147 (MQ=255)
                                                                        cTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAAt                                                                          <  2:1121505/150‑1 (MQ=255)
                                                                            cATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACtt                                                                      >  1:1020699/1‑150 (MQ=255)
                                                                             aTTTAAATTCCGTTTTTCCAGTTCAAATTCAATTGCCTTCAATGCAACTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCACCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTc                                                                     <  1:261206/150‑1 (MQ=255)
                                                                              tttAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGAAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTca                                                                    <  2:723895/150‑1 (MQ=255)
                                                                               ttAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTcac                                                                   >  2:236613/1‑150 (MQ=255)
                                                                                tAAATTCCGTTTTTCCAGTTCAAATGAAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCGCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCAGCAATACCTGACATTAATACTTcaca                                                                  <  2:319279/150‑1 (MQ=255)
                                                                                  aaTTCCTTTTTTCCAGTTCAAATGCAATTGCCTTCAATGAACCTTCGTAGCTGTTGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGTCAATACCTGACATTAATACTTcacatt                                                                <  2:642474/150‑3 (MQ=255)
                                                                                       cGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGCGCCAGCGGTGCTGGCTA‑TACCCC‑AATTAAGGATAAGAATGAAT‑TTTACGAGTTAATACCGTCGGCAAAACATGACAATAATAAATAACAGTCGCTg                                                           >  1:1034137/1‑150 (MQ=255)
                                                                                       cGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTg                                                           >  1:986154/1‑150 (MQ=255)
                                                                                        gTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTgg                                                          <  1:980228/150‑1 (MQ=255)
                                                                                        gTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTgg                                                          <  1:583996/150‑1 (MQ=255)
                                                                                          ttttCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTAGCTGGCa                                                        >  1:1144220/1‑150 (MQ=255)
                                                                                           tttCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCAt                                                       <  1:1271494/150‑1 (MQ=255)
                                                                                              ccAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGa                                                    <  2:217589/150‑1 (MQ=255)
                                                                                              ccAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGa                                                    <  1:1674948/150‑1 (MQ=255)
                                                                                               cAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAg                                                   <  1:1477218/150‑1 (MQ=255)
                                                                                                aGTTCAAATGCAATTGCCTTCAATGCACCTTCTTATCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATAAATACTTCACAGTCGCTGGCATCGAGc                                                  <  1:297380/150‑1 (MQ=255)
                                                                                                aGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCACCGGGGCTGGCTC‑TCCCCC‑CATTAAGGATAATAAATCAA‑TTTAAGAGTTAATACCGTCGGCCATACCGGACAATAATAAATAACAATAGCAGGCATCGAGc                                                  >  2:385737/1‑150 (MQ=255)
                                                                                                  ttCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTTGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAAAACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCAc                                                <  2:748133/150‑1 (MQ=255)
                                                                                                   tCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGCGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAAAACCGTCGGCAATACCAGACCTTAAAACAGCACAGTCGCAGGCATCGAGCaag                                               >  2:894179/1‑148 (MQ=255)
                                                                                                       aTGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑AATTACGGATAAGAATGCAT‑TTTACGAGTTAATACCGTCGGCAATAACTGACATTAATAAATAAAACTCGCTGGCATCGAGCACGGaaa                                           >  2:1609291/1‑150 (MQ=255)
                                                                                                         gCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCCGCGGTGGTGGCTC‑TCCCCC‑ATTTAAGGATAAGAATGCAT‑TTTCCGAGTTAAAACCGTGGGCAATTCCTGACATTAAAACGTCACAGTAGCTGGCCTAGAGCCCGGAAAAc                                         >  2:141996/1‑150 (MQ=255)
                                                                                                          cAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑AATTAAGGATAAGAAAGCAT‑TTTACGAGTTAAAACCCTCGGCAAAACCTTAAAATAATAAATAACAGTAGCTTGCCTCGAGCACGGAAAAca                                        >  2:453589/1‑149 (MQ=255)
                                                                                                              tGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAAt                                    <  1:1456856/150‑1 (MQ=255)
                                                                                                               gCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATc                                   <  2:665156/150‑1 (MQ=255)
                                                                                                                 cTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGa                                 <  1:236613/150‑1 (MQ=255)
                                                                                                                   tCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAg                               <  1:788308/150‑1 (MQ=255)
                                                                                                                            cTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACcac                      >  1:1124063/1‑150 (MQ=255)
                                                                                                                             ttCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACcaca                     <  2:1015107/150‑1 (MQ=255)
                                                                                                                                      gtggtgAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAAc            <  1:453589/150‑1 (MQ=255)
                                                                                                                                      gtggtgAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAAc            <  1:552565/150‑1 (MQ=255)
                                                                                                                                      gtggtgAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAAc            >  1:362307/1‑150 (MQ=255)
                                                                                                                                      ggggagAGCCAGCACCGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCGGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAAc            >  2:278721/6‑150 (MQ=255)
                                                                                                                                       gggagAGCCAGCACCGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGa                              >  2:1273074/5‑131 (MQ=255)
                                                                                                                                       gggagAGCCAGCACCGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGa                              <  1:1273074/127‑1 (MQ=255)
                                                                                                                                       gggagAGCCAGCACCGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAAcc           >  1:635309/5‑150 (MQ=255)
                                                                                                                                         gtgCGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCaa         <  2:1501750/150‑1 (MQ=255)
                                                                                                                                         gagAGCCAGCACCGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCaa         >  1:633492/3‑150 (MQ=255)
                                                                                                                                         gagAGCCAGCACCGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCaa         >  1:730257/3‑150 (MQ=255)
                                                                                                                                         gagAGCCAGCACCGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACACTCGCTGGCATCGAGCACGGAAAACTAAATCGAAGACGAACCACAGTTAATAACCaa         >  2:1281983/3‑150 (MQ=255)
                                                                                                                                          tgAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGCATGCAT‑TTTCCGAGTTAATACCGTCGGCAATTCCTGCCATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCaaa        >  1:77555/1‑150 (MQ=255)
                                                                                                                                          tgAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCaaa        <  1:491064/150‑1 (MQ=255)
                                                                                                                                          agAGCCAGCACCGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCaaa        >  1:122473/2‑150 (MQ=255)
                                                                                                                                             gCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAaaca     <  2:134618/150‑1 (MQ=255)
                                                                                                                                             gCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAaaca     <  2:886878/150‑1 (MQ=255)
                                                                                                                                              ccAGCGGTGCTGGATC‑TCCCCC‑ATTTACGTATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAaacaa    <  1:1110782/150‑1 (MQ=255)
                                                                                                                                                aGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAAc            <  1:867777/140‑1 (MQ=255)
                                                                                                                                                aGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAAc            >  2:867777/1‑140 (MQ=255)
                                                                                                                                                aGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATAAATACTTCACAGTCGCTCGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAAcc  <  1:57764/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCATCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTC‑TCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACC  >  NC_000913/3263353‑3263646

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑