Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,708,8800CA27.7% 277.7 / 35.6 184A654A (GCC→GCArsxCSoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex
Reads supporting (aligned to +/- strand):  ref base C (103/30);  new base A (49/2);  total (152/32)
Fisher's exact test for biased strand distribution p-value = 2.00e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 199 total aligned reads displayed.

ACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCT  >  NC_000913/1708730‑1709030
                                                                                                                                                      |                                                                                                                                                      
 cAGCAACATGCTATTGCTTAACAATAAGATAATCTCTTTCCGTGCAATTCCGTCTTCGTATCTGCTATTGCCCGTGCCAAAGCTCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGccgccg                                                                                                                                                       <  4:548488/151‑1 (MQ=2)
        aGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCGGCAGTGGCAGCGGCTATTGCCCGCGCCCACGCACGCAAGCGCGCCCAGACACCGGCAGATCACGAGGCAGAACACCAGGGAGATCCGCGCAAAGCGGCAGGCGacgg                                                                                                                                                >  1:1307401/1‑148 (MQ=9)
            tAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTGCCAGAAGAACCCGGGGATGCGCGCACGGCTGCCGTCGAGGCCGcg                                                                                                                                            >  3:616752/1‑150 (MQ=11)
             aatggggaacttttagaaCTGGTCGTTCCTCGCAATTCAGCTTTAGTAGCCGCTATTGGCTTTGCCAAATCTAGGATGTTGGAAGAGCAACAGGCTAATGCGGATCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTa                                                                                                                                           <  4:679519/137‑1 (MQ=2)
             aaTGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTGGACGCGGCTATTGCCCGAGCGAAGGCCAGCAACCGCGAACAGCACCAGGCTGATCCCGAGCCCGAACAGCAGGCGGATCCGCGCAACTCTGCAGTCGACGCCGCt                                                                                                                                            >  1:1105714/1‑150 (MQ=9)
              attatgAAGCAGAACAAGTGGTCGATCCTCGTAAATCAGCCGTCGAAGCAGCTATTGCCCGGGCCAAAGCTCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTAt                                                                                                                                          <  2:278805/146‑1 (MQ=11)
               tGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGGGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAAGCACGCCACAGCACAACACGGAGCTCCCCGCAACGCTCACGTCGAcgccgcactt                                                                                                                                         >  8:267516/1‑147 (MQ=2)
                         gaagaaCAGGTCGATCCGCGCAAAGCCGCCGTCGAATCCGCTATTGCCCGTGCCAAACCCCCCAAGCTGGACCAGCTACAGGCTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAAGCGGCCGTCGGAGCCGCtctggcgcgcgcac                                                                                                                               >  8:660065/1‑138 (MQ=14)
                           agaaCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCGAAAGCCCGCAAGCTGGAACAGCAACAGGCTAATGCCGTGCCAGAACACCAGGTGGATCCGCGCGAAGCCGCCGTCGACGACGCTATTGCCCgcgcgccg                                                                                                                             >  3:130840/1‑147 (MQ=11)
                                aGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCAGAGCCAACGGCGGCAAGCGGGAACCGCAACAGGCTAATGCCGCGCCAGAACACCAGGTAGATCCGCGCAACGCGGCCGTCGCAGCGGCGCTGGGCCGAGcccacggggg                                                                                                                        >  8:490775/1‑143 (MQ=2)
                                      aTCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTGGATCCGCGCAAAGCGGCAGTTGACGCCGCTATTGCCCgcgcacacgcccccaaccc                                                                                                                  >  6:1160198/1‑136 (MQ=14)
                                      aTCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGGGCCAAAGCACGCTAGGTGGAGCAGCACGCGGCTAAGGAGGCGCCAGACCTAAACGGGGATGCGCGCTAATCCGGCGTAGAAGCCGCTATTGGCGGCGCCCAGgcgcgcgcggg                                                                                                                  >  1:759301/1‑146 (MQ=2)
                                        accccgaacactcccgTAGAAGCCGCTATGCCCCGTGACAAACCTCGAAAGCTGGCAGATCAACATGATAATTACGACCCAGAACCACAGGGCGATCCTCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACgg                                                                                                                 <  6:294984/137‑1 (MQ=9)
                                               aaGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGCCGCGGCTATTGCCCGCGCGCACGCGCGCAAGCGCGCACAGc                                                                                                          >  5:35749/1‑150 (MQ=14)
                                                   cgccgTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCg                                                                                                     <  7:521370/151‑1 (MQ=255)
                                                    gccgTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCGGCTATTGCCCGAGCGAACGCCCGCAAACGCGAACAGCaccagg                                                                                                    >  6:9165/1‑146 (MQ=14)
                                                     ccgTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACACCAACAGGCTAATGCCGAGCCAGAAGAACAGGAAGATCCGCGCAAAGCCGCAGTAGAAGCCGCTATTGCCCGAGCAAAATCCACAAAACCCGAACAGCAACAGGc                                                                                                   >  2:1457442/1‑151 (MQ=11)
                                                     actgggAATCCGATATATCCCGCGCCAAATCGAGCAACGAGGCACAGCAACAGCCAAATCACGACCCAGAACCACATGGCGATCCCCGCAAACCCGCCGTCGAAGCAGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGc                                                                                                   <  1:980469/146‑1 (MQ=255)
                                                          gAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCGGCAGTGGCCGCCGCTATTGCCCGCGCCCAGGCGCGCAAACGGGCACAGCaccagggtcatg                                                                                              >  1:923581/1‑140 (MQ=14)
                                                              ccGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAAGCCGCAGTAGAAGCCGCTATTGCCCGCGACAAATCACGCAACCGGGAACAACAACAGGCTAATgatgt                                                                                          >  4:1198224/1‑147 (MQ=11)
                                                                gCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAACGGGACCAGCCCCCGGCTAATGCGgcgc                                                                                        >  8:72232/1‑151 (MQ=255)
                                                                 cTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAAGCAGCAGTAGAAGCCGCTATTGCCCGCGCAAAAGCGAAAAAACCCGAACAGCAGCAGGATAAATCGGaac                                                                                        >  7:1163134/1‑148 (MQ=14)
                                                                  tATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAATCCGAAGACCAACCGGATATTGAACAAGCGAAAGCAAGCAATATGGAACAGCAACCGGCGAAtcagtaaaa                                                                                       >  5:1834387/1‑142 (MQ=9)
                                                                  tATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAAGCGGCAGTAGAAGCCGCTATTGCCCGAGCAAAAGCCCGCAAACGGGACCAGCAAAAGGCTAATCCgtaacc                                                                                       >  4:633250/1‑145 (MQ=14)
                                                                  tATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAAGCCGCAGTAGAAGCCGCTATTGCCCGAGCCAAAGCGAACAACCGAAAACAGAAAGCGACTAATACGGAgac                                                                                       >  4:1581114/1‑148 (MQ=14)
                                                                     tGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAAGAAATCGTATCACCGACCAAAGCCGAAATCGAAGACGCTATTGACATAGACAAAGCGCGCAAACTGTAACATAAAAAGGCTAATGCGGAGCCaga                                                                                    >  5:521973/1‑150 (MQ=9)
                                                                      gCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCAGTCGAAGCCGCTCTAGCACCAGCCCAAGCCCGCACACTCCCACAGCAAGAGGCGAATGCGGCGCCagaac                                                                                  >  1:1172239/1‑150 (MQ=14)
                                                                      ccccgctccAAAGCCCGCAAGCTGGAACAGCAACAGGATAATGCGGACCCAGAACAACAGGTCGATCCGAGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCagaag                                                                                  <  4:80002/144‑1 (MQ=11)
                                                                      ccccccgcccccccccGCAAGCGGCAACAGCAACATCCTAATTCCCAGCCACATCACCACGTCGATCCTAGGACAGGGGACTGGCCAGCCGCTGTGGGCGGTTCCAATGCGCGCAACCTGGAACAGCAACCGGCTAATGCGGAGCCagaag                                                                                  <  3:1603167/136‑1 (MQ=255)
                                                                       cccGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCGGACGTCGCAGCCGCTCTTGGCGGAGCGAAAGCCCCCCAGCGGGAACAGCTACGGGGTCTTGCGGTTCCaga                                                                                    >  8:928065/1‑148 (MQ=14)
                                                                       cccGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAGCGGGAACCGCACCCGGCTCCTGCGGAGCCagc                                                                                    >  1:229411/1‑147 (MQ=14)
                                                                        ccGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCAGTAGAAGCCGCTATTGCCCGAGCAAAAGCGAGAAAACCCGAAAAGTAACAAGATAATCCGTAGCCagaaga                                                                                 >  6:423486/1‑150 (MQ=14)
                                                                          gTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAAGCCGCAGTAGACGCCGCTATTGCCCGCGCCAAAGCGCGCAAACGGGAACAGCACCAGGCTGATGCGGAGCCagaagacc                                                                               >  1:1443855/1‑148 (MQ=14)
                                                                             ccAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAAGCCGCCGTAGAAGCCGCTATTGCCCGCGCGAACGCGAGAAAACGCGACCAGCACAAGGATATTTAGGAGCCAGAAGAACAgg                                                                            >  8:1467609/1‑150 (MQ=14)
                                                                                acgccctcAAGCTGGAACAGCATCAGGCTAACGCGGACGCAGCACATCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGa                                                                        <  8:1006531/144‑1 (MQ=14)
                                                                                aaGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCAGTCGAAGCCGCTATTGCCCGCGCCCACGCGCGAAAACGGGAACAGCACCAGGCTAATGCGGAGCCAGAAGTCCAGGTTgc                                                                        >  6:659423/1‑150 (MQ=14)
                                                                                   cACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAAGCCGCAATATAAGAAGATATTCCACGAGCCAAAGCGCGGAAGCTGTAAGAGAAAACTGCGAATGCGGAGACAGAAGAAAAAGATGATcc                                                                     >  2:1762869/1‑151 (MQ=9)
                                                                                     cGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCAGTAGAAGCGGCTATTGCCCGCGCGAAGGCCCACAAAGCCGAACAGCAGAAGGATATTTCGGAGGAAGAAAAACAAGtacataagc                                                                   >  7:1834547/1‑143 (MQ=2)
                                                                                       caatctttaaCAGCACCCGCCTAATGCCGAGCCACTCCAAGCGGTCTTTCGGCGATAAGGGGGCTTCGAAGCTGTTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCgcgc                                                                 <  6:1827115/143‑1 (MQ=2)
                                                                                        aaGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCAGTGGCCGCGGCTATTGCCCGAGCGAACGCCAGCAAAGCGGACCAGCACCAGGGTAATGAGCCGGCAGAAGAGCACCGTGCTCCgcgc                                                                 >  6:1798951/1‑150 (MQ=2)
                                                                                            tGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTAGCTACGCGCAAACCACCAGTCGAAGAAGCTATTACCCGAGACAAAGCAAGAAAACGGGAACAGCAACAGGATAATGAAGAGCCAGAAGAACAGGTTggtgcgggcgaggg                                                            >  4:292930/1‑138 (MQ=9)
                                                                                             ggAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAAGCCGCAGTGGCAGCCGCTATTGCCCGCGCAAAAGCCACCAAACGGGAACAGCAACAGGATCATGAGGAGCCAGAAAACCAGGTAGATAcgcctaaatcc                                                           >  3:1113355/1‑143 (MQ=9)
                                                                                               aaCAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTGGAAGCGGCTATTGCCCGAGCCAAGGCCCGCAAACGGGACCAGCACCAGGCTAAAGCCGAGCCAGAAGAACACGTAGATCACCGCACAGCtcc                                                         >  1:1642947/1‑149 (MQ=14)
                                                                                                aCAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCGGCAGTCTAAGCCGCTATTGCCCGCGCCAACGCGAGAAAACCGGAACAGCACCAGGCTAATGCGGAGCCAGAAGACCAGGTTGATAAGCGCACAGcgccc                                                        >  6:1337701/1‑147 (MQ=14)
                                                                                                 cAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCCATCCGCGCAAAGCCGCCGTAGAAGCGGCTATTGCCCGAGCGCACGCCACCAAACCGGCACAGCAACAGGCTTAACACGAGCCCGAAGAAGAGGAAGATCCGCCCAAAGATGccc                                                       >  1:1522024/1‑150 (MQ=255)
                                                                                                  aGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGACCAGGTTGATCCGCGCAAAGCTGCCGt                                                      >  3:448031/1‑151 (MQ=255)
                                                                                                      aCAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACAGGCTAATGCGGAGCCAGAAGACGAGGTTGATCCGCGCAAAGCGGCCGTGGac                                                  >  3:1032027/1‑150 (MQ=255)
                                                                                                        aGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAACGCCCGCAAACGCGCCCAGCACCCGGCTAATGCGGAGCCAGAAGACCAGGTTGCTCCGCGCAAAGCGGCCGTGGCAGc                                                >  6:147620/1‑151 (MQ=14)
                                                                                                         cagctaTGCCGAGCCAGAACAAGTGTTTTATCATCTCAACGCGGCCTTCGATTCCTCTATTGCCGGATACAAAGTTGGCATACTGGAAGATCATCCGTTTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCg                                               <  2:1746846/146‑1 (MQ=255)
                                                                                                          gCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAACGCGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCGGCCGtggccgcgg                                              >  6:1704785/1‑143 (MQ=14)
                                                                                                           cTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGACCCAGAACGAGAGGTCGAGCCGCCCTCTGCAGCCCTCGCAGCggg                                             >  3:579169/1‑150 (MQ=11)
                                                                                                            tAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGAGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCAGTCGAAGCGGCt                                            >  8:841752/1‑151 (MQ=255)
                                                                                                            tAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAACCGGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCAGTAGAAGCGGCt                                            >  1:327498/1‑151 (MQ=255)
                                                                                                             acccccGACCCAGAACACCCCGTCGATCCGCGCACACCCGCACGCGACGCCCCTATTGCACGAGCCAACCCTCGCAATCTGTAACCGCAACAGGCTAACGCGGACCCAGAAGAACAGGTCGATCCGCGCTAAGCTGCCGTCGAAGCGGCTa                                           <  1:1107053/147‑1 (MQ=2)
                                                                                                             aaTGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACACCAACCGGCTAATGCGGACCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTAGCCGCGGCt                                            >  3:1536187/1‑150 (MQ=255)
                                                                                                              aTGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGGTATTGAACGCGACACCGCGCGCAAACGGGCACCGCACCCGGCTAATGCGGCGCGAGACGCCGAGGTTGATCCGCGCCTCGCTGCCGGCgccgcgcggtg                                          >  1:1748798/1‑141 (MQ=9)
                                                                                                              aTGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATCCGAACAAAGAAGACGACGAAAAGGCTAt                                          >  4:1777887/1‑151 (MQ=11)
                                                                                                                 cggaaccAGAAGAACAGGTCGATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTAACCGAACCAAATCGAGCAAGCAAAAACAGCAACCTGCTAATGAGGAGGCAGAAGAACACGTAGATAAGCGCAAAGCAGCCGTCGAAGCGGCTAtggg                                       >  2:1110556/6‑148 (MQ=2)
                                                                                                                  cGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAACCGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGCTCCGCGCAACGCGGCAGTAGCCGCGGCTATTGcc                                      >  5:1134115/1‑150 (MQ=11)
                                                                                                                   gAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCAGTAGCAGCGGCTATTGCAc                                     >  5:1546653/1‑151 (MQ=255)
                                                                                                                     gCCAGAACAACAGGTCGATCCGCGCAAAGCCGCAGTAGCAGCCGCTATTTCGCAAGCCAAACCGCGAAACATCGAACAGAAACCGGCGACTGCCGATCCAGAAGACCAGGATGATACACGAAAAGCTGCCGTGGCCGCGGCTagggcgcgc                                   >  5:118313/1‑143 (MQ=9)
                                                                                                                      ccAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCCCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCATAGCTGCCGAGGAAGCGGCTATTGCCcgcc                                  >  3:1232323/1‑150 (MQ=255)
                                                                                                                        aGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATCAGCGCAAAGCTGCAGTCGCAGAGGCTATTCCCCGCGcc                                >  6:1097229/1‑151 (MQ=14)
                                                                                                                           accccaGCTCGATCCGCGCAAAGCCGCCTTTTAAGCCGCTATTGCCCGCGCCAAAGCGCTCAAACTTGAACATCAACCTTCTAATGCGGTTCCATAATAAGAGGTCGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCaaa                             <  4:927842/147‑1 (MQ=255)
                                                                                                                            caacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCAGTAGACGCGGCTATTGCACGCGCaaacg                            >  8:566655/1‑149 (MQ=17)
                                                                                                                            caacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCCCGCAAACGGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATAAGCGCAAATCGGCCGTAGAAGCGGCTATTGCACGCGCGaagg                            >  5:884614/1‑149 (MQ=14)
                                                                                                                             aacaGGTCGATCCGCGCAAAGCCGCCGTAGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAACGGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCAGTAGCAGCGGCTATTGCACgcgctcacgc                           >  4:1823690/1‑145 (MQ=14)
                                                                                                                              acaGGTCGATCCGCGCAAAGCCGCAGTCGAAGCCGCTATTGCCCGCGCCAAAGCGAGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATAATCGCAAAGCTGCAGAAGAAGCGCCTAATCCACGCGCAAAAGc                           >  7:1499533/1‑150 (MQ=17)
                                                                                                                                ccgccGATCCGCGCAAAGCCGCCGTCGAAGCAGCTATTGCCCGAGCCAAAGCTCGCAAACTGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACTGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACg                        <  7:526247/147‑1 (MQ=255)
                                                                                                                                aGGTCGATCCGCGCAAAGCCGCCGTCGACGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCACCCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTGGACGCGGCTATTGCCCGCGCCACCGCACg                        >  3:251508/1‑151 (MQ=17)
                                                                                                                                aGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCAAAAGCCAGCAAACGGGAACAGCAACAGGCTAATCCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCAGTAGAAGCGGCTATTGCCCGCGCAAAAGccac                        >  2:1624700/1‑148 (MQ=14)
                                                                                                                                aGGTCGATCCGCGCAAAGCCGCAGTAGAAGCCGCTATTGCCCGAGCAAAAGCCAGAAAACGCGAACAGCAACAGGCTAATGAGGAGCAAGAAGAACAGGTTGATAAGCCAAAATCACCAGTATAATCGGCTATTGCAcgagaaaacgcac                         >  4:1409851/1‑139 (MQ=9)
                                                                                                                                 ggTCGATCCGCGCAAAGCCGCAGTCGAAGCCGCTATTGCCCGAGCAAAAGCAAGAAAGCAGGAACAGCAACAGGCTAATGACGAAACAGAAGAACAGCTTGATCCGCGAAAAGCTGCCGTAGACGCGGCTATTGCCCgcgcgcaggcgca                        >  5:1631942/1‑141 (MQ=2)
                                                                                                                                   tCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGCAGCGGCTATTGCCCGCGCCAAcgcacgca                      >  6:1546194/1‑150 (MQ=255)
                                                                                                                                   tCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGACGCGGCTATTGCACGCGCCAAcgcacgca                      >  8:459112/1‑150 (MQ=255)
                                                                                                                                    ccccccccccAAACCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCCCTCTCCATCGAACATAAAGAGGCTTCTCCGGAGCCACAACCACAGGTCTATCGGCGAAAACCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCa                      <  1:445592/140‑1 (MQ=2)
                                                                                                                                    cGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAATCGGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCAGTAGAAGCGGCTATTGCCcgcggctcggcagggcgg                    >  5:87218/1‑137 (MQ=11)
                                                                                                                                    cGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCAGCGCAAAGCTGCCGTGGCCGCGGCTATTGCACGCGCCAAGgcaagcaac                    >  7:1402533/1‑150 (MQ=17)
                                                                                                                                    cGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCAGTAGAAGCGGCTATTGCCCGCGCCACGGCACggagg                    >  4:466011/1‑147 (MQ=17)
                                                                                                                                    cGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCAGCGCAAAGCTGCAGTCGCAGCGGCTATTGCACGCGCCAACGCACGaaac                    >  6:1035538/1‑150 (MQ=255)
                                                                                                                                     gATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGc                   >  6:618325/1‑151 (MQ=255)
                                                                                                                                      aTCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGCAGCGGCTATTGCACGCGCCAACGCAAGCAAGCt                  >  2:708219/1‑151 (MQ=255)
                                                                                                                                      aTCCGCCCAAACCCGCCGTCGAAGCCGCTATTCGCCGAGCCAAAGCGCGCAAACGGCCACATCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCt                  <  6:1634600/151‑1 (MQ=25)
                                                                                                                                        ccGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTgg                >  8:1324632/1‑151 (MQ=255)
                                                                                                                                        ccGCGCAAAGCCGCAGTCGAAGCCGCTATTGCCCGAGCCAAAGCAAGCAATCTGGAACAGCAAAAGGCTAAAGCGGAACCAGAAAAACATGTAGATCCGCGAAAAGCAGCAGTAGACGCAGCTATTACACGAGATAAAGCaacaaagcacg                >  7:643770/1‑141 (MQ=255)
                                                                                                                                         cgcgCAACGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCTCGCAATCTGGAACAGCAACCGGCTAATGCGGATCCAGAAGAACAGGTTGATCCGCGCAAAGCTTCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGa               <  3:525725/151‑1 (MQ=14)
                                                                                                                                         cgcgCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGCCGCGGCTATTGCACGCGCGAAAGCCAGCAAGCTggc               >  5:1270063/1‑150 (MQ=255)
                                                                                                                                          gcgcAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTAGAAGCGGCTATTGCACGCGCGAAAGCAAGAAAGcacga               >  7:776444/1‑146 (MQ=17)
                                                                                                                                           cgcAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCCCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCAGTAGCCGCGGCTATTGCCCGCGCCAAAGCaacaaagcaggcac             >  2:917770/1‑138 (MQ=17)
                                                                                                                                            gcAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGGTGGCAc             >  1:68299/1‑150 (MQ=255)
                                                                                                                                            gcAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTAGAAGCGGCTATTGCCCGCGCAAAAGCACGCAAGCAGGacc             >  2:1776692/1‑148 (MQ=255)
                                                                                                                                            gcAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCAGTCGAAGCGGCTATTGCCCGCGCCCATGCACGCAAGCAGGAAc             >  6:491516/1‑150 (MQ=255)
                                                                                                                                              aaaGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGACGCGGCTATTGCACGCGCGAACGCCCGCAAGCTGGACCAGc          >  2:200397/1‑151 (MQ=255)
                                                                                                                                              aaaGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGc          >  8:1624000/1‑151 (MQ=255)
                                                                                                                                              aaaGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCAGTAGAAGCGGCTATTGCACGCGCGAAAGCAAGCAACCTCGAACAGc          >  2:1266950/1‑151 (MQ=17)
                                                                                                                                              aaaGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACAGGCTAATGCGGAGGAACAAGAGCACGTGGATCAGCGCAAAGCTGGAGTCGAAGCTGGTATTGCACGAGCCAGCGCACGCAAGATGGAACAGc          >  5:667548/1‑151 (MQ=14)
                                                                                                                                                   cgccgTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGCAGCGGCTATTGCCCGCGCGAACGCCCGCAAGCTGGAACAGCacc       >  8:113059/1‑147 (MQ=17)
                                                                                                                                                    gccgTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCAGTCGAAGCGGCTATTGCACGCGCCAAAGCCAGCAAGCTCGACCAGCAACAgg    >  7:366219/1‑151 (MQ=255)
                                                                                                                                                     ccgTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCCCGCAAGCTGGAACAGCACCAgg  >  3:115995/1‑150 (MQ=255)
                                                                                                                                                      cgTCGCAGCCGCGACTGCCCGACCCACCGCGCCGCAACCGCAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTGGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCt  <  3:1764150/151‑1 (MQ=14)
                                                                                                                                                      |                                                                                                                                                      
ACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCT  >  NC_000913/1708730‑1709030

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: