Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,708,7580GC24.8% 199.3 / 44.0 159E614Q (GAA→CAA) rsxCSoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex
Reads supporting (aligned to +/- strand):  ref base G (42/73);  new base C (4/34);  total (48/111)
Fisher's exact test for biased strand distribution p-value = 2.06e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 168 total aligned reads displayed.

GCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAG  >  NC_000913/1708608‑1708908
                                                                                                                                                      |                                                                                                                                                      
ccccGCGCCAAACCTCGCAAGCTGGAACACCACCAGGCTACTGAGGAACCAGCACCACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCagaag                                                                                                                                                        <  8:111115/150‑1 (MQ=14)
 aaggcgaaaaaaggggcAAGCTGGAAGACCTTCAGGCTACTGCGGAACCAGAACAACATGTCGATCCCCGCAAATCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCagaag                                                                                                                                                        <  4:1590130/135‑1 (MQ=14)
   cgcgACAAATCGCGCAAGCTGGAACAGCAACAGGCTAATCCGGAACCAGCACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAAc                                                                                                                                                     <  2:202113/151‑1 (MQ=255)
    ccgcCACAGCGCGCAAGCTGGAACAGCAACAGGCTAATCCGGAACCAGAACAACAGGGCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAgaagaaaa                                                                                                                                                    <  3:835127/150‑3 (MQ=255)
     ggcCAAAGCGCGCAAGCTGGAACAGCAACAGGATAATGCGGAACCAGAACAACAGGTCGATCCTCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAg                                                                                                                                                   <  3:1465963/150‑1 (MQ=255)
     ggactaatcgcggAAGCTGGAAGAGCAACAGGCTAATGCGGAACCAGAACAACTGGGCGATCCGCGCAAAGCCGCCGTGGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAg                                                                                                                                                   <  4:200289/143‑1 (MQ=14)
      gtcaaatgtggcAAGCTGGAAGATCTTCAGGATAATCAGGAAGCAGAAACACATGTTGATCCTCGTAAAGCGGCCGTTGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAgg                                                                                                                                                  <  1:737312/141‑1 (MQ=9)
       ccAAATCTCGCAACCTGGAAGAGCAACAGCATAATCAGGACGCAGAAACACAGGCCGATCCTCGGAAATCCGCCGTTGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGt                                                                                                                                                 <  1:168731/151‑1 (MQ=14)
       cataatctcggaatcTGGAAGATCATCAGGATACTTCGGACGCAGAAAACGTGGGCGATCCTCGGAAAGCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGt                                                                                                                                                 <  6:1066388/137‑1 (MQ=9)
        cAAAGCGCGGAAGCTGGAACAGCAACAGGCTAATGCGGAAGCAGAACAACAGGGCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTc                                                                                                                                                <  8:291625/151‑1 (MQ=255)
          aaGCGCGCAAGCTGGAACAGCAACAGGATAATGCGGAACCAGAACAACATGGCGATCCGCGCAAAGCGGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGa                                                                                                                                              <  5:1206657/151‑1 (MQ=255)
            gcgcgcAAGCTGGAACATCAACAGGCTAATGCGGAACCAGAACACCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATc                                                                                                                                            <  3:1691248/151‑1 (MQ=255)
            ccgcgcAAGCTGGAAGACCAACAGGATAATGCGGACGCAGAAAAACTGGTGGATCCGCGCAAAACTGCCGTAGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATc                                                                                                                                            <  3:1642938/150‑1 (MQ=14)
            ccccgcccgCTGCAACAGCAACAGGCTAATGCGTATAATGAAGAAGAGGTCGCGACGCGCAAAGGAGCCCTGGCAGCCGATTTCCGCCGTGCGATCTCGCGATATCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATc                                                                                                                                            <  7:408569/143‑1 (MQ=2)
             gtcggaaGCAGGAACACCAACAGGATACTCACGAACCAGCACCACTTGTCGATCCTCGGAAAGCCGCCGTGGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAGTGCGGAACCAGAAGAACAGGTCGATcc                                                                                                                                           <  5:1346975/146‑1 (MQ=14)
              gcggAATCTGGAACATCAACAGGATAATGAGGAACCAGAAAAACATGTCGATCCTCGCAAAGCCGCCGGCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCcg                                                                                                                                          <  2:555364/151‑1 (MQ=14)
              gcgcAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGACCAAGTGGGCGATCCTCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCcg                                                                                                                                          <  7:1015249/151‑1 (MQ=255)
               gggaaGCTGGAAGATCAACAGGCTACTGCGGAAGCAGAACAAGAGGTCGATCCTCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCcgc                                                                                                                                         <  1:1095803/148‑1 (MQ=14)
               ctgaatcTGGAAGAGCTAGAGGCTTATTCGGTACCAGAACATGTGGTCGATGCTCGCATAGCGGCGGTGGAAGCGGGTTTTGCGGGTGGGAAAGGGCGGAAGCTGGAACAGCATCAGGGTTATGGTGAACCAGAAGAACAGGTCGAt                                                                                                                                             <  2:622331/141‑1 (MQ=255)
                gcAATCTGGAACAGCAACAGGATAATGCGGAACCAGCACACGAGGGCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCcgcg                                                                                                                                        <  6:389645/151‑1 (MQ=14)
                gcAAGCTGGAACAGCAACAGGATAATCCGGAACCAGAACACCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGGGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGG‑ACCAGAAGAACAGGTCGATCCgcgc                                                                                                                                       <  4:576504/151‑1 (MQ=14)
                ccaatcaggAACACCAACAGGAAACTCCCGAACCAGAAACCCAGGCCGATCCCCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCcgcg                                                                                                                                        <  6:1125868/144‑1 (MQ=14)
                     cTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTAGATCAGCTCAAAACAGCCGACTAAACCGAAATAGACAATGCCAAAGAGAACAACAAGGAACAGCACCAAGCAAATAAAAAAAAAGAAGAACAGGTGGATCCgcgcgg                                                                                                                                     >  1:409414/1‑147 (MQ=9)
                      tGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTAGATCCGCGCATAGCATCCGGCGAAGCCTCTAGTGCACTCGCCAAAGACCGGAAGCTGGACCAGCAACAGGCTTATGAGGAACCAGAAGAACAGGTGGCTCcgcggagcgg                                                                                                                                  >  3:350410/1‑145 (MQ=9)
                       ggAAGATCATCAGGATAATTAGGAAGCAGAAAAAGTTGCTGATCCTCGTAAATCTGCCGTTGAAGCAGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAgcc                                                                                                                                 <  1:70103/151‑1 (MQ=14)
                       ggAACATCATCAGGATAATCAGGAACCAGAACAAGTTGTGGATCCTCGTAAAGCTGCCGTTGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAg                                                                                                                                   <  1:883797/149‑1 (MQ=14)
                       ggAACAGCAACAGGCTAATGCGGAAGCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAgcc                                                                                                                                 <  2:598298/151‑1 (MQ=255)
                       ggAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAgc                                                                                                                                  <  3:70301/150‑1 (MQ=40)
                       ggAACAGCAACAGGCTAATCCGGAAGCAGAACACCAGGGCGATCCGCGCAAAGCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCaa                                                                                                                                     <  7:1613284/147‑1 (MQ=255)
                         aaCAGCAACAGGATAATTCGGAACCAGAACCACAGGGCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATAGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAgccgc                                                                                                                               <  3:243869/151‑1 (MQ=255)
                          aCATCAATAGGCTAATGTGTTTCCATAAAATAAGGTCTATCGTTGTAAAGCTGCCGTCGAAGCGGCTATTGCACGAGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAgcggca                                                                                                                              <  1:1620203/151‑2 (MQ=14)
                           cAGCATCAGCCTAATGCGGAACCAGAACAAGAGGGCGATCCGCGTAAACCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGccgccg                                                                                                                             <  1:870276/151‑1 (MQ=14)
                            aGCAACAGGCAAATGCGGAACCAGAAAACCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGt                                                                                                                            <  3:1114520/151‑1 (MQ=255)
                             acAACAGGCTACTCACGAACCAGAACAAGATGGCGATCCTCGGAAAGCAGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTc                                                                                                                           <  3:1596242/150‑1 (MQ=14)
                              catcAGGATACTTCGGAAGCAGCACAACTGGGGGATCCCCGTAAATCTGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTc                                                                                                                           <  6:1635973/147‑1 (MQ=11)
                              cAACAGGCTAATGCGGAACCAGAACAACAGGCCGATCCGCGCAAATCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTc                                                                                                                           <  4:228116/150‑1 (MQ=255)
                              cAACAGGATAATGACGACGCAGCACCAGAGGCCGATCCTCGCAAACCCGCCGTGGAAGCCGCTATGGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCg                                                                                                                          <  4:1241069/151‑1 (MQ=11)
                               ttcagcatactccGGAAGCAGAACACGTTGTTGATCCACGCAAAGCTGCCGTTGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGccgccg                                                                                                                             <  6:613878/135‑1 (MQ=11)
                               aaCAGGCTAATCCGGAACCAGCAACCCAGGCGGATCCGCGAAAAACCGCCGTCGAAGCAGCTATCGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGa                                                                                                                         <  1:97444/151‑1 (MQ=14)
                                 caggaaaatcccgAACCAGCAAAACAGGGCGATCCCCGCAAACCCGCCGTCGAAGCCGCTATTGCCCGGGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAg                                                                                                                       <  5:1206812/139‑1 (MQ=14)
                                   gcataatcacgAAGCAGAACACGATGTCGATCCTCGTAAATCTGCCGTTGAAGCAGCTATTGCCCGGGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGcc                                                                                                                     <  6:900202/141‑1 (MQ=11)
                                           ccgAACCAGAACACCATGCCGATCCCCGAAAATCCGCCGTGGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGc                                                                                                             <  5:1298859/149‑1 (MQ=14)
                                              aaCCAGAACACCAGGGCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCg                                                                                                          <  3:639013/151‑1 (MQ=255)
                                                ccAGAACACCAGGTCGATCCGCGAAAAACCGCCGTAGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTg                                                                                                        <  3:1151592/151‑1 (MQ=255)
                                                       aacaGGTCGATCCGCGCAAAGCCGCCGTCGCAGCCGCTATTGCCCGAGCAAAAGCCAGCAAGCTGGAACAGCAGCAGTCTAATGCGCAACCAGAAGACCAGGTCTATCCTCGCAAAGAGGCCGAAGCCGCCGCAATTGccagaccaaa                                                                                                    >  3:558283/1‑140 (MQ=2)
                                                         catgcggATCCGCGTAAAGCGGCCGTAGAAGCCGCTATTGCCCTTGCCATATCGCGCAAGCTGGAACAGCAACAGGATAATCCGGAACCAGAACAACAGGGCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCAc                                                                                               <  4:884543/145‑1 (MQ=11)
                                                          cccgcgctCCCCGCAAACCCGCCGTTGTATCATCTATTGCCCGGGCCAAAGCGCGCAATATTGAACAGCAACATGCTAATGCGTAACCTGATCTACAGGGCTATCCCCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACg                                                                                              <  8:1354282/144‑1 (MQ=2)
                                                           ggTCGATCCGCGCAAAGCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGc                                                                                             <  3:947274/151‑1 (MQ=255)
                                                             gcGATCCAAGCAAAGCCGCCGTAGAAGCATCTATTGCTCGTGCCACAGCGCCCAAGCGGGAACAGCAACAGGATAATCCGGAACCAGAACAACAGGCCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCaa                                                                                           <  3:573128/150‑1 (MQ=11)
                                                                agccgcggtaacctgccgTCGAAGCAGCTATTGCCCGGGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACCACAGGGCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCt                                                                                        <  7:132515/137‑1 (MQ=11)
                                                                    gcgcAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATCCGCGAAAACCCAACATCGAAGACTATATTAACCGTGAAAAAGAAAGAAAGCTAGa                                                                                     >  1:825756/1‑150 (MQ=11)
                                                                     gggaaaaCCGCCTTCGAATCTGATATCTCCCGTGTCAAAGTCCTCTAGCTGGAACAGAAATAGGATAATGAGGAACCAGAACAAGAGGCTGATCCTCGCAAACCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAAc                                                                                   <  6:541093/148‑1 (MQ=2)
                                                                     cgcAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGCAGCCGCTATTGCCCGCGCGAACGccacaaagctggcac                                                                                   >  3:597363/1‑137 (MQ=14)
                                                                      ccAACGCCGCCGGCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACa                                                                                  <  4:192874/150‑1 (MQ=255)
                                                                         aaGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAGAACGCTAAAGCGTAGACAGAAGAAAAGGTGGATCCGATAAAAGAGGACTTCGAAGACGCAATTGCAATAGACAAAGCACGCAAGCTGGAACAGc                                                                                >  8:1160932/1‑150 (MQ=2)
                                                                          aGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACATCCTAATCCCGACCGAGAAGAAAAGGTAGATCCGCTAAAAGAAGCAGTCGAAGCGGCGAGTGACAGTGAAAAAGAACGCAAGCTGGGACGGccg                                                                              >  6:1298413/1‑149 (MQ=2)
                                                                             ggccgTCGAAGCAGCTATTGCCCGGGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACACCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAg                                                                           <  7:74867/150‑1 (MQ=9)
                                                                              gccgTAGAAGAAGCTATTTCCCGGGACAAATCTCGGAAGCTGGAACATCAACAGGCTAATGCGGAAGCAGAAGAAGAGGTCGATCCTCGCAAAGCCGCCGTCGAAGCCGCTAATGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACa                                                                            <  2:1602168/149‑1 (MQ=11)
                                                                                 gTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATCCGCGCAAAGCCGCAGTAGAAACAGCTAAAGCCCGTCCCAAAGCACCCAAGCTGGAACAGAAACAGACt                                                                        >  2:656803/1‑150 (MQ=14)
                                                                                   cGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTGGATCCGCGCAAAGCGGCAGTGGACGCGGCTATTGCCCGCGCCCACGCCCGCAAGCGGGCCCAGCAGCAGGCtctt                                                                     >  3:1709117/1‑148 (MQ=14)
                                                                                   cGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCAGTAGCAGCCGCTATTGCCCGAGCCCAGGCCCGCAAGCCCGCCCAGCACCAGGCTGAt                                                                     >  5:1519242/1‑151 (MQ=14)
                                                                                     aaGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTAGATCCGCGCAAAGCGGCCGTGGACGCGGCTATTGCCCGCGCGAACGCCCGCAAGGGCGAACAGCAGCAGGCTaaacc                                                                   >  5:397169/1‑148 (MQ=14)
                                                                                          gCTATTGCCCGTGCCAAAGCGCGCAATCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAgct                                                                                                                                 <  1:1058861/84‑2 (MQ=11)
                                                                                          gCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAgct                                                                                                                                 >  7:1058362/1‑83 (MQ=255)
                                                                                          gCTATTGCCCGAGACAACGCGCGCAAGCTGGAAGAGCATCAGGATAATCCGGAACCAGAAGACGAGGTCGATCCTCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGc                                                              <  1:553409/151‑1 (MQ=14)
                                                                                             aTTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATCCGCGCAAAGCGGCAGTAGAAGCGGCTATTGCCCGAGCGCACGCAATCAACCCCGACCAGCACCAGGTTCTAGCCGAGCCat                                                           >  7:198040/1‑150 (MQ=2)
                                                                                              ttGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCAAATCCGAAGACAGAAGAAGAGGGAAATACGCTAAGAGAAGCAGTAGAACGCGCTAGTGCCCGTGACAAATAAAAAAATCTGGAACAGCACCCGGcggcgggggggcgcgg                                                          >  1:144307/1‑136 (MQ=2)
                                                                                                    gggacaaaagccgcAAGCTGGAACACCAACAGCATACTCCGGAACCAGAACCACATGCCGATCCCCGCAAACCCGCCGTGGAAGCCGCTATCGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAacaaca                                                    <  7:881815/140‑1 (MQ=11)
                                                                                                    ccgacAAATCCCGCAATCTGGAAGAGCAACAGGCTAATCCGGACCCAGAACAACAGGCCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAacaaca                                                    <  4:354466/147‑1 (MQ=11)
                                                                                                      gCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAAGCCGCAGTAGAAGCCGCTATTGCCCGAGCAAAAGCCAAAAAACCCGAACAGCAAAAGGATAAAAACGAGCCATAaaaacaat                                                  >  1:503818/1‑149 (MQ=9)
                                                                                                       aataatgtgggaatcTGGAAGACCATCAGGAAAATTAGGAACCAGAACAAGATGTTGATCCTCGTAAATCTGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGt                                                 <  2:187368/137‑1 (MQ=2)
                                                                                                         aaaGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATCCGCGCAAAGCAGCAGTAAGATCCGATAATCCACGTGGAAAAGAACTCAAGATGGAACAGAAAAACTCTAATGCGGACCCCGAACAAAAgaaa                                                >  6:364066/1‑147 (MQ=2)
                                                                                                         aaaGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTGGATCCGCGCAAAGCGGCAGTAGAAGCGGCTATTGCCCGCGCGCACGCCAGAAAACCCGAACAGCACCAGGCTCTAGACGAGCCCGaagacctagtcc                                               >  7:1461929/1‑141 (MQ=0)
                                                                                                          aacggcgcAATCTGGAAGAGCACCAGGATAATCCCGAACCAGAACAACATGGCGATCCTCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGa                                              <  4:1485898/147‑1 (MQ=11)
                                                                                                           aGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATCCGCGCAAAGCCGCAGTAGAAGCAGATATTAACAGAGACAAAGAACGAAAGATTGAACACAAAAATGAAAATCACGAGCCAcaacaaaacataaa                                              >  2:1186172/1‑142 (MQ=2)
                                                                                                               cggaatcaggAAGACCAACAGCATAATTCCGAACCAGAACAAGTTGGTGATCCTCGTAAACCTGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGAt                                             <  2:262242/139‑1 (MQ=2)
                                                                                                                ggaatcTGGAACATCTTCAGGATAATCCGGACGCAGCACCAGTGGGCGATCCCCGGAAACCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCcgcg                                        <  8:821284/146‑1 (MQ=11)
                                                                                                                   agcgtgAACAGCAACAGGATAATCAGGACCCAGAACAAGTGGGCGATCCGCGAAAAGCGGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCgcgc                                       <  5:1425833/144‑1 (MQ=11)
                                                                                                                    ccTGGAACACCAACAGGAAAATCCCGAACCAGAACCCCAGGCCGATCCCCGCAAACCGGCCGTCGAAGCCGCTATTGCCCGGGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCaaa                                    <  5:227849/150‑1 (MQ=11)
                                                                                                                     cagtaaCATCATCAGGAAAATTAGGAACCAGAACAACATGCCGATCCTCGCAAATCAGCCGTAGAAGCAGCTATAGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAg                                   <  6:496882/147‑1 (MQ=11)
                                                                                                                       ggaagacctaCAGGATACTTCGGAACCAGCAACACTTGGCGATCCCCGCAAACCAGCCGTCGAAGCAGCTATTGCCCGCGCCAAAGCCCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgcc                                 <  5:1061379/142‑1 (MQ=11)
                                                                                                                       ggAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGGAGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgcc                                 <  2:953227/151‑1 (MQ=255)
                                                                                                                            aGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTGGCCGCCGCTATTGCCCGCGCGAAAGCTCGCAAGCGCGAACAGGAACAGAATAATGGAGGGCAAAAAGAACATGTAGATCCTCGCCAAgagcaagt                            >  6:341235/1‑144 (MQ=2)
                                                                                                                              caccaggaaaaTGCCGAACCAGCACCACATGCCGATCCCCGTAAACCCGCCGTAGAAGCAGCTATCGCCCGCGCCAAAGCTCGCAAGCTGGAACAGCAACAGGCTAATGCCGACCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCg                          <  4:31316/142‑1 (MQ=11)
                                                                                                                              cAACAGCCAAATCCCGACCCAGCACACCAGGCTGATCCCCGCAAACCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCg                          <  7:735550/151‑1 (MQ=11)
                                                                                                                                   ggCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGCCGCCGCTATTGCCCGGGCCAAAGCACGCAAGCGGGACCAGCACCAGGCTCATCCCGAGCCAGAACACCAGGTCGCTCCGCGCACAGCCGCCGttgccgcg                     >  3:136307/1‑144 (MQ=14)
                                                                                                                                   ggCTAATGCGGAACCAGAAGAACAGGTAGATCCGCGCAAAGCAGAACTCAAGCCCGAAGTTTACCGAAACAAAGCACGCAATATGGTACATCTACATGCTAAAGCGGTCCAATAAAAAAAGGTTGATCCGCGCGAAgcggggggccagggg                     >  1:1131013/1‑138 (MQ=255)
                                                                                                                                    gCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCAAACGCCCGCAAGCGCGAACAGCAACAGGCTAATGCCGAGCAAGAACAACACGTCGATCCGCGCAACGCCGCCGTGGCAGccc                    >  7:542015/1‑150 (MQ=14)
                                                                                                                                    gCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCAGTCGAAGCCGCTATTGCCCGAGCCCAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGGACCAGAAGAACAGCGTGCTCCGCGCAAAGCCGCAGTCGAAGCCg                    >  4:767027/1‑151 (MQ=11)
                                                                                                                                    gCTAATGCGGAACCAGAAGAACAGGTAGATACGCGAATAGCAGACGTCTAAGCGGCTATAGCACTAGGCGAAGCCCTAACCCTGCAACAGCAACATCATAAACCCGGGCCAGAAGAAGAGGGGGAGCCGCGCGGCGCCGCCGGCGGAGcgg                    >  4:233478/1‑149 (MQ=255)
                                                                                                                                     cTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCAGTAGAAGCCGCTATTGCCCGAGCAAAAGCACGAAACCTCGAACAGCAACAGGATAATGCCGAGCCAGAAAAACAGGTAGATCCGAGCAAATCCGCAGTAGAAGCCGc                   >  2:1227002/1‑151 (MQ=14)
                                                                                                                                      tAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGCAGCGGCTATTGCCCGCGCCCACGCCCGCAACCACGCCCAGCAGCAGGCTAATGCCGAGCCAGAACACCACGTCGCACCGCGCACAgccgcagagcacgcctc                   >  3:1688391/1‑138 (MQ=9)
                                                                                                                                      tAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCAGTAGAAGCCGCTATTGCCCGCGCAAAAGCAAGCAACCACGAACAGCAACAGGCTACAGCCGACCCAGAACCACACGTCGCTAAGCGCCAAgccgaagtcccaccgcct                  >  7:702686/1‑137 (MQ=9)
                                                                                                                                          gCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCGGCCGTCGACGCCGCTATTGCCCGTGCGAACGCCAGAAAGCGGGAACAGCAGCAGGGTCATGCCGAGACAGCGGAAGAGGACCATCCGCGCAAAGCGGCAGtttacgcgcctactg              >  3:972058/1‑137 (MQ=255)
                                                                                                                                           cGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCGGCCGTCGAAGCGGCTATTGCCCGGGCCAAAGCGCGCACAATGGAACACCAACAGCCTAATGCCGCGCCAGAAGAACAGGTAGCTCCGCGCCACGCGGCAGTAGAAGCGGCTATTGc             >  7:311293/1‑151 (MQ=11)
                                                                                                                                           cGGAACCAGAACAAGAGGTGGATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTg              <  2:390281/150‑1 (MQ=255)
                                                                                                                                                  agaagaACAGGTCGATCCGCGCAAAGCCGCAGTGGACGCCGCTATTGCCCGAGCGAAAGCCAAAAACCGCGCACAGCACCAGGCTAAAGACGAGCCAGAACACGAGTTCTATCAACCCAACGCCGCCGTAGCCGCGGCTATTCCCCgccc       >  8:499766/1‑147 (MQ=2)
                                                                                                                                                   gaagaaCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGGGCCAACGCACGCAACCGGGAACAGCACCAGGGTAATGCCGAGCCCGAACAACAGGTCGCTCCGCGCACAGCCGCCGTGGCCGCCGCTATTGCCCgcgca      >  3:395057/1‑149 (MQ=14)
                                                                                                                                                    aagaaCAGGTCGATCCGCGCAAAGCCGCAGTCGAAGCCGCTATTGCCCGAGCAAAAGCCCGCAACCAGGAACAGCAACAGGCTAAAGCCGAGCCAGAAAAACAGGTAGATCAGCGCAAAGCAGCCGTAGCAGCAGCTATTGCCCGAGCCa     >  2:396290/1‑150 (MQ=14)
                                                                                                                                                      gaaCAGGTCGATCCGCGCAAAGCGGCCGTGCAAGCCGCTATTCCCCGAGCGAAAGCAAGCAAGCTGGAACAGCAGCAGGCTAATGCCGACCCAGAAGAACACGTCCATCCGCCCAACGCCCCAGTTGACCCCGCGATTGcccccgcaacgc  >  4:329360/1‑142 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
GCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAG  >  NC_000913/1708608‑1708908

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