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.

TCATCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACC  >  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)
           aatGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTcc                                                                                                                                         <  1:1549465/150‑1 (MQ=255)
           aatGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCACCGCTGGCTCAcc                                                                                                                                         <  1:1073929/150‑1 (MQ=255)
            atGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTccc                                                                                                                                        <  1:1435953/150‑1 (MQ=255)
             tGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCTAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAACGGTGCTGGatcaaccc                                                                                                                                        >  2:1456856/1‑145 (MQ=255)
                  gACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCATTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGCGCCAGCCGTTGTGGGACTGCCCCAttt                                                                                                                                 >  2:1586679/1‑150 (MQ=255)
                    cTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTAc                                                                                                                               <  1:212578/150‑1 (MQ=255)
                    cTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCACCCGTGCTGGCTCTCCCCCATTTAc                                                                                                                               >  2:506301/1‑150 (MQ=255)
                      cGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACgg                                                                                                                             <  2:1552305/150‑1 (MQ=255)
                        gTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGATCCCGCGGTGCTGGAACTCCCCCAATTTACgtt                                                                                                                            >  2:568247/1‑148 (MQ=255)
                           aaaaTACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATACCGGATTTCCAGTTCAAATGCAATAGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCAAGCGGTGCTGGCTCACCCACATTTACGGATAAg                                                                                                                        >  1:54956/1‑150 (MQ=255)
                            aaaTACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGa                                                                                                                       >  2:552565/1‑150 (MQ=255)
                            aaaTAATGCCGCCGTGAGCGATGCGTGGGCCAATTAGGGCCACCATGGCTTTTAAATGCGTTTTTGCCATTTAAATTGAAATTGCATTCAATGCCCCTTCTTAGCGTTGGTGAGCAAGCGGTCATGGATCTCCCCCATTTACGGATacta                                                                                                                       >  2:1216430/1‑147 (MQ=255)
                                aCTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGGTGGGTCGCCCCCAATTACGGATAAGAATgg                                                                                                                   >  2:374799/1‑149 (MQ=255)
                                    ccgccgTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCCCCGGGGGGGGGGGTCCCCCATTTTAGGATTAGAATGGAttt                                                                                                              >  1:1063222/1‑150 (MQ=255)
                                    ccgccgTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCGATCCCCCCATTAACGTTAATAATGGAttt                                                                                                              >  2:597037/1‑150 (MQ=255)
                                     cgccgTGCGCGATGCGGTGGCCAATTAAGGCCACACTGTCAATTAAATTCCGTTGTGCCAGTTCAACTGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCGCCCCCATTTACGGATAAGAATGCAtttt                                                                                                             <  1:287197/150‑1 (MQ=255)
                                     cgccgTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCAtttt                                                                                                             >  1:665156/1‑150 (MQ=255)
                                      gccgTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTc                                                                                                            <  1:990538/150‑1 (MQ=255)
                                       ccgTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTcc                                                                                                           >  2:57764/1‑150 (MQ=255)
                                         gTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGCGCCAGGGGTGCTGGGGCTCCCCCAATTACGGGTTAGAATGCATTTTGCgg                                                                                                         >  1:414324/1‑149 (MQ=255)
                                         gTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGa                                                                                                         >  2:22533/1‑150 (MQ=255)
                                         gTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGa                                                                                                         <  2:162955/150‑1 (MQ=255)
                                            aGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCTTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGtt                                                                                                      <  1:1453424/150‑1 (MQ=255)
                                            aGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGtt                                                                                                      <  1:1111166/150‑1 (MQ=255)
                                                   cGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCg                                                                                               <  1:1019000/150‑1 (MQ=255)
                                                    ggTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCCGCGGGGGTGGCTCTACCCCATTTTCGGATTATAATGCATTTTACGAGTTAATAACgg                                                                                              >  1:1703449/1‑149 (MQ=255)
                                                      tGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCg                                                                                            >  1:1015107/1‑150 (MQ=255)
                                                      gggCCAAGTCAGGCCTCACTGTCATTTAAATTCCGTTTTTCCAGTTCAACGGCAATTGCCTTCAATGCACCATCGTAGCTGTGGTGCGCCAGCGGTGCTGGCTCTCCGCCATTTACGCATAAGAATGCATTTTCCGAGTTAATACCGGCg                                                                                            <  2:314684/149‑1 (MQ=255)
                                                       ggCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCgg                                                                                           >  1:920587/1‑150 (MQ=255)
                                                        gCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGc                                                                                          <  2:195188/150‑1 (MQ=255)
                                                        gCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGc                                                                                          <  1:825709/150‑1 (MQ=255)
                                                         ccAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCGTCCCCATTTACGGATAAGAATGCATTTTCCGCGTTAATACCGTCGGcc                                                                                         >  1:37240/1‑149 (MQ=255)
                                                           aaTTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAAt                                                                                       >  1:36153/1‑150 (MQ=255)
                                                                aGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTg                                                                                  <  2:842179/150‑1 (MQ=255)
                                                                 ggCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGa                                                                                 <  1:1448372/150‑1 (MQ=255)
                                                                  gCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTTGAGGCTCAAACCCCATTTACGATTAGAAATCAATTTACGAGATTATACCCTAGGGAAAACCTGaa                                                                                >  2:1477218/1‑149 (MQ=255)
                                                                  gCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGAc                                                                                <  2:161830/150‑1 (MQ=255)
                                                                  gCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGAc                                                                                <  2:1453338/150‑1 (MQ=255)
                                                                  gCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGAc                                                                                <  1:606572/150‑1 (MQ=255)
                                                                     acacTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCTGTGCTGGCTCTCCCCCATTTACGGATAAGAAAGCATTTTCCGAGTTAATACCGTCGGCAATTCCTGACAtt                                                                             >  1:242259/1‑150 (MQ=255)
                                                                      cacTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGCGCCACCGGTGGGTGCGGACCCCCAGTTAAAGAAAAGAATGCATGTTAACAGGTTAAAACCGGGGGAATACCTGACagtt                                                                            >  2:871129/1‑147 (MQ=255)
                                                                        cTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAAt                                                                          <  2:1121505/150‑1 (MQ=255)
                                                                            cATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACtt                                                                      >  1:1020699/1‑150 (MQ=255)
                                                                             aTTTAAATTCCGTTTTTCCAGTTCAAATTCAATTGCCTTCAATGCAACTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCACCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTc                                                                     <  1:261206/150‑1 (MQ=255)
                                                                              tttAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGAATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTca                                                                    <  2:723895/150‑1 (MQ=255)
                                                                               ttAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTcac                                                                   >  2:236613/1‑150 (MQ=255)
                                                                                tAAATTCCGTTTTTCCAGTTCAAATGAAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCGCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCAGCAATACCTGACATTAATACTTcaca                                                                  <  2:319279/150‑1 (MQ=255)
                                                                                  aaTTCCTTTTTTCCAGTTCAAATGCAATTGCCTTCAATGAACCTTCGTAGCTGTTGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGTCAATACCTGACATTAATACTTcacatt                                                                <  2:642474/150‑3 (MQ=255)
                                                                                       cGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGCGCCAGCGGTGCTGGCTATACCCCAATTAAGGATAAGAATGAATTTTACGAGTTAATACCGTCGGCAAAACATGACAATAATAAATAACAGTCGCTg                                                           >  1:1034137/1‑150 (MQ=255)
                                                                                       cGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTg                                                           >  1:986154/1‑150 (MQ=255)
                                                                                        gTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTgg                                                          <  1:980228/150‑1 (MQ=255)
                                                                                        gTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTgg                                                          <  1:583996/150‑1 (MQ=255)
                                                                                          ttttCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTAGCTGGCa                                                        >  1:1144220/1‑150 (MQ=255)
                                                                                           tttCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCAt                                                       <  1:1271494/150‑1 (MQ=255)
                                                                                              ccAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGa                                                    <  2:217589/150‑1 (MQ=255)
                                                                                              ccAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGa                                                    <  1:1674948/150‑1 (MQ=255)
                                                                                               cAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAg                                                   <  1:1477218/150‑1 (MQ=255)
                                                                                                aGTTCAAATGCAATTGCCTTCAATGCACCTTCTTATCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATAAATACTTCACAGTCGCTGGCATCGAGc                                                  <  1:297380/150‑1 (MQ=255)
                                                                                                aGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCACCGGGGCTGGCTCTCCCCCCATTAAGGATAATAAATCAATTTAAGAGTTAATACCGTCGGCCATACCGGACAATAATAAATAACAATAGCAGGCATCGAGc                                                  >  2:385737/1‑150 (MQ=255)
                                                                                                  ttCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTTGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAAAACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCAc                                                <  2:748133/150‑1 (MQ=255)
                                                                                                   tCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGCGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAAAACCGTCGGCAATACCAGACCTTAAAACAGCACAGTCGCAGGCATCGAGCaag                                               >  2:894179/1‑148 (MQ=255)
                                                                                                       aTGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCAATTACGGATAAGAATGCATTTTACGAGTTAATACCGTCGGCAATAACTGACATTAATAAATAAAACTCGCTGGCATCGAGCACGGaaa                                           >  2:1609291/1‑150 (MQ=255)
                                                                                                         gCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCCGCGGTGGTGGCTCTCCCCCATTTAAGGATAAGAATGCATTTTCCGAGTTAAAACCGTGGGCAATTCCTGACATTAAAACGTCACAGTAGCTGGCCTAGAGCCCGGAAAAc                                         >  2:141996/1‑150 (MQ=255)
                                                                                                          cAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCAATTAAGGATAAGAAAGCATTTTACGAGTTAAAACCCTCGGCAAAACCTTAAAATAATAAATAACAGTAGCTTGCCTCGAGCACGGAAAAca                                        >  2:453589/1‑149 (MQ=255)
                                                                                                              tGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAAt                                    <  1:1456856/150‑1 (MQ=255)
                                                                                                               gCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATc                                   <  2:665156/150‑1 (MQ=255)
                                                                                                                 cTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGa                                 <  1:236613/150‑1 (MQ=255)
                                                                                                                   tCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAg                               <  1:788308/150‑1 (MQ=255)
                                                                                                                            cTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACcac                      >  1:1124063/1‑150 (MQ=255)
                                                                                                                             ttCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACcaca                     <  2:1015107/150‑1 (MQ=255)
                                                                                                                                      gtggtgAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAAc            <  1:453589/150‑1 (MQ=255)
                                                                                                                                      gtggtgAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAAc            <  1:552565/150‑1 (MQ=255)
                                                                                                                                      gtggtgAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAAc            >  1:362307/1‑150 (MQ=255)
                                                                                                                                      ggggagAGCCAGCACCGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCGGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAAc            >  2:278721/6‑150 (MQ=255)
                                                                                                                                       gggagAGCCAGCACCGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGa                              >  2:1273074/5‑131 (MQ=255)
                                                                                                                                       gggagAGCCAGCACCGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGa                              <  1:1273074/127‑1 (MQ=255)
                                                                                                                                       gggagAGCCAGCACCGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAAcc           >  1:635309/5‑150 (MQ=255)
                                                                                                                                         gtgCGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCaa         <  2:1501750/150‑1 (MQ=255)
                                                                                                                                         gagAGCCAGCACCGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCaa         >  1:633492/3‑150 (MQ=255)
                                                                                                                                         gagAGCCAGCACCGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCaa         >  1:730257/3‑150 (MQ=255)
                                                                                                                                         gagAGCCAGCACCGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACACTCGCTGGCATCGAGCACGGAAAACTAAATCGAAGACGAACCACAGTTAATAACCaa         >  2:1281983/3‑150 (MQ=255)
                                                                                                                                          tgAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGCATGCATTTTCCGAGTTAATACCGTCGGCAATTCCTGCCATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCaaa        >  1:77555/1‑150 (MQ=255)
                                                                                                                                          tgAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCaaa        <  1:491064/150‑1 (MQ=255)
                                                                                                                                          agAGCCAGCACCGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCaaa        >  1:122473/2‑150 (MQ=255)
                                                                                                                                             gCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAaaca     <  2:134618/150‑1 (MQ=255)
                                                                                                                                             gCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAaaca     <  2:886878/150‑1 (MQ=255)
                                                                                                                                              ccAGCGGTGCTGGATCTCCCCCATTTACGTATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAaacaa    <  1:1110782/150‑1 (MQ=255)
                                                                                                                                                aGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAAc            <  1:867777/140‑1 (MQ=255)
                                                                                                                                                aGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAAc            >  2:867777/1‑140 (MQ=255)
                                                                                                                                                aGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATAAATACTTCACAGTCGCTCGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAAcc  <  1:57764/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCATCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACC  >  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: