Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 163,419 A→T 5.8% I439F (ATT→TTT)  hrpB → putative ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913163,4190AT 5.8% 214.8 / 2.6 78I439F (ATT→TTT) hrpBputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base A (38/35);  new base T (3/2);  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 = 8.83e-01

GAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCGGTGAGGCAGACAGTTCGCTTA  >  NC_000913/163271‑163566
                                                                                                                                                    |                                                                                                                                                   
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:41864/1‑149 (MQ=255)
                                                                                                                                        aCCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCGGTGAGGCAGa             <  1:282537/149‑1 (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)
                                                                                                                                                    |                                                                                                                                                   
GAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCGGTGAGGCAGACAGTTCGCTTA  >  NC_000913/163271‑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: