Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 163,418 A→T 6.6% K438N (AAA→AAT hrpB → putative ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913163,4180AT 6.6% 220.7 / 3.0 78K438N (AAA→AAThrpBputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base A (38/34);  new base T (3/3);  total (41/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCGGTGAGGCAGACAGTTCGCTTA  >  NC_000913/163270‑163566
                                                                                                                                                    |                                                                                                                                                    
cgAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCaaaa                                                                                                                                                      <  2:430526/149‑1 (MQ=255)
cgAAACGGCTGTTACAAATGCCGGGGGGCCTGGAGGGGGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGGTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCaaaa                                                                                                                                                      <  2:1057657/149‑1 (MQ=255)
 gAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCaaaaa                                                                                                                                                     >  1:167510/1‑149 (MQ=255)
  aaaCGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAAt                                                                                                                                                    >  1:1100140/1‑149 (MQ=255)
    aCGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGTACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTg                                                                                                                                                  >  2:293772/1‑149 (MQ=255)
    aCGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTg                                                                                                                                                  >  1:931432/1‑149 (MQ=255)
    aCGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTg                                                                                                                                                  >  2:719635/1‑149 (MQ=255)
      gccTGTTACAAATGCTGGGGGCACTAGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTgcc                                                                                                                                                <  1:648420/147‑1 (MQ=255)
      gTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTgcc                                                                                                                                                <  2:355682/149‑1 (MQ=255)
           ttACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCAt                                                                                                                                           >  1:1079387/1‑149 (MQ=255)
            tACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTATAGCTCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCATTGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAATTTTTgccgcggta                                                                                                                                          <  2:124177/149‑5 (MQ=255)
             aCAAATGCTGGGGGCACTGGAGGGTGACCGGCTTAGTGCGCTAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATtc                                                                                                                                         >  1:713098/1‑149 (MQ=255)
               aaaTGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATtctc                                                                                                                                       >  2:282537/1‑149 (MQ=255)
                aaTGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGAAGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCCTTCACg                                                                                                                                      >  2:807566/1‑149 (MQ=255)
                aaTGCTGGGGGCAATGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAATATGGGAGCGCTGGGTATCGATCCGCGGTTAGCGGCAATGGTGGTTAGCGCGAAGATCGACGACGAAGCTGCTACAGCGGCAAAATTTGCCGCCATTCTCg                                                                                                                                      <  2:177970/149‑1 (MQ=255)
                 aTGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGa                                                                                                                                     >  2:256630/1‑149 (MQ=255)
                 aTGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGa                                                                                                                                     >  2:863536/1‑149 (MQ=255)
                  ggCTGGGGGCACTGGAGGGTGAACGGCGTAGTGCGCAAGGGCAAAAAATGGCAGCGCGGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGACAAAAATTGCCGCCATTCTCGaa                                                                                                                                    <  1:1005790/148‑1 (MQ=255)
                      gggggCACTGGAGCGTGAACGGCTTAGTGCTCAAGGGCAAAAAATGGCCGCGCTGGGTAACGATCCGCGTTTAGCGGAAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAgc                                                                                                                                <  1:307978/149‑1 (MQ=255)
                            aCTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGACGCCATTCTCGAAGAGCCGCAAc                                                                                                                          >  2:1164106/1‑149 (MQ=255)
                             cTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTTCCGCCATTCTCGAAGAGCCGCCACg                                                                                                                         >  2:306773/1‑149 (MQ=255)
                                  gggTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATgg                                                                                                                    >  2:590594/1‑149 (MQ=255)
                                   ggTTAACGGCTTAGTGCGCAAGGGCAAAAAATCGCACCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATggg                                                                                                                   <  1:43133/149‑1 (MQ=255)
                                   ggTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATggg                                                                                                                   <  1:1135328/149‑1 (MQ=255)
                                     ggAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCa                                                                                                                 <  2:720588/148‑1 (MQ=255)
                                            cTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCGACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGa                                                                                                          <  1:66928/149‑1 (MQ=255)
                                               aGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCt                                                                                                       >  1:943360/1‑149 (MQ=255)
                                                 tGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTgg                                                                                                     >  2:556410/1‑149 (MQ=255)
                                                  gcgcAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTggg                                                                                                    >  1:817043/1‑149 (MQ=255)
                                                    gcAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGggcg                                                                                                  >  2:181232/1‑149 (MQ=255)
                                                         ggCAAAACAGGGCAGCGCTGGGTAACGATCCGCGTGTAGCGGCAATGCTGGTTAGCGCGAAGACCGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGggcgtggcg                                                                                             <  1:132149/149‑1 (MQ=255)
                                                          gCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGggcgtggcgt                                                                                            <  2:1043887/149‑1 (MQ=255)
                                                                actGGCAGCGCTGGGTAACGATCCGCGGTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTcgc                                                                                      <  1:68244/147‑1 (MQ=255)
                                                                aaTGGCAGCGCTGGGTAACGATCCGCGGTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTcgc                                                                                      <  1:420622/149‑1 (MQ=255)
                                                                    gCAGCGCGGGGCAACGATCCGCGTTTAGCGGCAATGCTCGTAAGCGCGAAGAACCACGCCGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCCTAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCaa                                                                                  <  2:491374/149‑1 (MQ=255)
                                                                     cAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATTGTGACCTGGGCGTGGCGTTTTCGCGCAAt                                                                                 >  2:1029584/1‑149 (MQ=255)
                                                                         gcTGGGTAACGATCCGCGTTTAGCGGAAATGCTGGTTAGCGCGAAGAACGACTACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGAAATAGTGACCTGGGCGTGGCGTTTTCGCGCTATCAAc                                                                             <  2:911916/149‑1 (MQ=255)
                                                                                 aCGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTgg                                                                     <  1:339180/149‑1 (MQ=255)
                                                                                  cGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGgc                                                                    <  2:54509/149‑1 (MQ=255)
                                                                                   gATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGgca                                                                   >  2:750106/1‑149 (MQ=255)
                                                                                   gATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCAGGGCGTGGCGTTTTCGCGCAATCACCCCGCATGgcc                                                                   >  1:1104575/1‑148 (MQ=255)
                                                                                   gATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGCATGGGCAATAGTGACCTGGGCGTGGAGTTTTCGCGCAATCAACCAGCCTGgca                                                                   >  2:992032/1‑149 (MQ=255)
                                                                                    aTCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGgcag                                                                  >  2:367869/1‑149 (MQ=255)
                                                                                     tCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGgcagc                                                                 >  2:606467/1‑149 (MQ=255)
                                                                                       cgcgTTTTGCGGGATTGCTGGTTAGCGCGATGAACGACGAGGATGCTGATACCGCGGCATAAATTGCCGCCATTCTGGAAGAGCCTCAACGGATGGGAAATAGTGACCTGGGCGTGGCGTTTTTGCGCAATCAACCAGCCTGGCAGCaa                                                               <  2:1147821/149‑1 (MQ=255)
                                                                                       cgcgTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGTATCTGGTACCGCGGCAAAAATTGCCGCCATTCTCGACGTGCCGCCACGGATGTGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCaa                                                               <  2:744394/149‑1 (MQ=255)
                                                                                       cgcgTATAGCGGCAATGCTGCTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCaa                                                               <  1:697050/149‑1 (MQ=255)
                                                                                        gcgTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAAc                                                              <  2:167510/149‑1 (MQ=255)
                                                                                        gcgTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAAc                                                              >  1:239847/1‑149 (MQ=255)
                                                                                          gTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGGTTTCGCGCAATCAACCAGCCTGGCAGCCACGt                                                            >  1:941179/1‑149 (MQ=255)
                                                                                           gttAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTa                                                           <  1:37822/148‑1 (MQ=255)
                                                                                              aGCGGCAATGCTGGTAAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTc                                                        <  1:90826/149‑1 (MQ=255)
                                                                                                 ggCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGc                                                     >  2:151736/1‑149 (MQ=255)
                                                                                                  gCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCa                                                    >  1:18546/1‑149 (MQ=255)
                                                                                                      tGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGAAATAGGGACCTGGGACTGGCGTTTTCGCGCAAACAACCAGCCTCGCAGCAGCGTAGTCAGCAAcgg                                                >  2:113003/1‑147 (MQ=255)
                                                                                                       gCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTTGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGt                                               <  2:1056877/149‑1 (MQ=255)
                                                                                                        gtGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGtt                                              <  2:943360/148‑1 (MQ=255)
                                                                                                          ggTTAGCGCGAAGAACGTCGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCACCCAGCCTGGCAGCAACGTAGTCAGCAACTGTTaa                                            <  1:73837/149‑1 (MQ=255)
                                                                                                           gTTAGCGCGAAGAACGACGACGAGGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTaaa                                           >  1:889292/1‑149 (MQ=255)
                                                                                                           gTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCCCGGATGGGCAATAGTGACATGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTaaa                                           >  2:365674/1‑149 (MQ=255)
                                                                                                             tAGCGCGAAGAACGACGACGTAGCTGCTACCGCGGCAAATATTGCCGGCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAAc                                         <  2:782031/149‑1 (MQ=255)
                                                                                                             tAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAAc                                         <  1:494432/149‑1 (MQ=255)
                                                                                                                  cgAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTa                                    >  1:922934/1‑149 (MQ=255)
                                                                                                                    aagaaCGCCGACGAAGCTGCTACCGCGGCAATAATTGCCGCCATTCTCGCAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCCGCAACGTAGTCAGCAACTGTTAAAACGCTTaaa                                  >  1:968531/1‑149 (MQ=255)
                                                                                                                     agaaCGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAac                                 <  2:431207/149‑1 (MQ=255)
                                                                                                                     agaaCGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAac                                 <  1:284117/149‑1 (MQ=255)
                                                                                                                     agaaCGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAAACAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAac                                 <  1:295160/149‑1 (MQ=255)
                                                                                                                         cgacgacgaAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAacgtac                             >  1:379425/1‑149 (MQ=255)
                                                                                                                         cgacgacgaAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAacgtac                             <  2:503033/149‑1 (MQ=255)
                                                                                                                             gacgaAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTgg                         <  2:698821/149‑1 (MQ=255)
                                                                                                                                gaAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCgg                      >  1:1032542/1‑149 (MQ=255)
                                                                                                                                gaAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCgg                      >  1:356020/1‑149 (MQ=255)
                                                                                                                                gaAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGGTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAAGTGTTAATACGCTTAAACGTACGTGgggg                      >  1:967527/1‑149 (MQ=255)
                                                                                                                                      agaatggcgGCAATTTTTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCCGCAACTGTTAAAACGCTTAAACGTACGTGGCGGTGAGGc                >  1:1054335/7‑149 (MQ=255)
                                                                                                                                      agaatggcgGCAATTTTTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCGGTGAGGc                >  2:438805/7‑149 (MQ=255)
                                                                                                                                      agaatggcgGCAATTTTTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCGGTGAGGc                >  2:413886/7‑149 (MQ=255)
                                                                                                                                         aCCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCGGTGAGGCAGa             <  1:282537/149‑1 (MQ=255)
                                                                                                                                         aCCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCGGTGAGGCAGa             >  1:41864/1‑149 (MQ=255)
                                                                                                                                           cgcgGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAAGGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCGGGGAGGCAGACa           >  1:912227/1‑149 (MQ=255)
                                                                                                                                              ggCAAATATTGCCGCCATTCTCGCAGAGCCGCCCCGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCATAAACGTACGTGGCGGTGAGGCAGACAGtt        <  1:225181/149‑1 (MQ=255)
                                                                                                                                                cAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCGGGGCGGGGCGTTTTCGCGCAATCAACCAGCCAGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCGGTAAGGCAGACAGTTCg      <  2:713098/149‑1 (MQ=255)
                                                                                                                                                 aaataTTTCCGCCATTCTAGAAGAGCCTCCACTGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCGGTGAGGCAGACAGTTCGc     <  1:187617/149‑1 (MQ=255)
                                                                                                                                                 aaaaaTTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCGGTGAGGCAGACAGTTCGc     <  1:556410/149‑1 (MQ=255)
                                                                                                                                                 aaaaaTTGCCGCCATTCTCGAAGAGCCGCCACGGAAGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCGGTGAGGCAGACAGTTCGc     <  2:1144321/149‑1 (MQ=255)
                                                                                                                                                  aaaaTTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCGGTGAGGCAGACAGTTCGCt    >  2:20675/1‑149 (MQ=255)
                                                                                                                                                   aaaTTGCCGCCATTCTCGAAGAGCCGCCACGGATGGTCAATAGTGACCTGGGCGTGTCGATTTAGAGCATTAAACCAGCATGTCAGCAACGTAGTCAGCAAATGTTAAAACTATTAAACGTACGTGGCGGTGAGGCAGACAGTTCGCtt   >  2:516262/1‑149 (MQ=255)
                                                                                                                                                    aaTTCCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGGGGCGGTGAGGCAGACAGTTCGCTTa  <  2:1085590/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCGGTGAGGCAGACAGTTCGCTTA  >  NC_000913/163270‑163566

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: