Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 162,603 A→T 6.7% M167L (ATG→TTG)  hrpB → putative ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913162,6030AT 6.7% 181.6 / 2.9 69M167L (ATG→TTG) hrpBputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base A (35/29);  new base T (2/3);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.14e-01

ACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGAAAGCGGATCATTAC  >  NC_000913/162458‑162750
                                                                                                                                                 |                                                                                                                                                   
aCTGGTGATCCTTGATGAATTTAATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGAAGTGCAACAAGGTCTGCGTGATGACCTTAAACTGATGATTATGTCGGCTACCCTGGACCACGACCGCTTGCAGCAAAtgc                                                                                                                                                  >  2:31229/1‑149 (MQ=255)
    gTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGcc                                                                                                                                              >  1:138779/1‑149 (MQ=255)
    gTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGcc                                                                                                                                              >  1:155297/1‑149 (MQ=255)
    gTGATCCTTGAGGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGATGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTCCTGATTATGTCGTCTACCCTGGACACCGACCGCTTGCTGCCAATGCTGcc                                                                                                                                              >  2:1017019/1‑149 (MQ=255)
      gATCCTTGATGCATTTCATGAGCGCAGCTTGCAGGCGGATTTGGAGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAg                                                                                                                                            <  1:1008102/149‑1 (MQ=255)
      gATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCATTTGCTGCAAg                                                                                                                                            <  2:1096097/149‑1 (MQ=255)
         ccTTGAAGAATTTCATGCGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGGGCAACAAGGTCTGCGGGATGACCTTAAACTGCTGATTATGGCGGCTACCCTGGACAACGACCGCTTGCAGCATTTGCTGCAAgcgg                                                                                                                                         <  1:1072476/149‑4 (MQ=255)
          cTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAgc                                                                                                                                        >  1:132982/1‑149 (MQ=255)
           ttGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAgcg                                                                                                                                       >  2:372140/1‑149 (MQ=255)
                gaATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAATGGCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTgt                                                                                                                                  <  2:1172469/149‑1 (MQ=255)
                 aaTTTGATGAGAGCTGCTTGCTGGCGGATTTGGCGTTGGAGCTGTTTCTCGATGTGGAACTTGGTCTCCGTTATGACCTTAAACTGCTGATTATGTCGGCTAACCTGGACAACGACCGCTTTGAGCAAATGCTGCCAGAAGCGCCTgtc                                                                                                                                 <  2:155443/149‑1 (MQ=255)
                 aaTTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTgtc                                                                                                                                 <  1:755332/149‑1 (MQ=255)
                  aTTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGCCCGCTTGCAGCAAATGCTGCCAGAAGCGCCTgtcg                                                                                                                                >  1:493749/1‑149 (MQ=255)
                  aTTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGACGCAAATGCTGCCAGAAGCGCCTgtcg                                                                                                                                >  2:23983/1‑149 (MQ=255)
                     tCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCGGCGGGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCGTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCa                                                                                                                             <  1:17240/149‑1 (MQ=255)
                      cATGAGCGCAGCTTGTAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGACTGCGTGATGACCTTAAACTGGTGATTATGTCGGCTACCCTGGACAACGACCGATTGCAGCAAATCATGCCAGAAGCGCCTGTCGTCAt                                                                                                                            >  1:594512/1‑149 (MQ=255)
                       aTGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCAtc                                                                                                                           <  1:24937/149‑1 (MQ=255)
                         gAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCAtctc                                                                                                                         <  2:155297/149‑1 (MQ=255)
                               aGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAACGTCTGCGTGATGACCTAAAACTGCTGATTATGTCGGATACCCAGGACAACGAACGCTTGCAGCAAATGCTTCAAGAAGCGGCTGTCCTCATCTCAGAAgg                                                                                                                   >  2:513375/1‑149 (MQ=255)
                                 cTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGgc                                                                                                                 >  2:133185/1‑149 (MQ=255)
                                 cTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGgc                                                                                                                 <  2:822776/149‑1 (MQ=255)
                                  ttGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGGCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAATGgcg                                                                                                                <  2:349405/149‑1 (MQ=255)
                                          ggATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACGGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTc                                                                                                        <  1:748708/149‑1 (MQ=255)
                                          ggATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTc                                                                                                        >  1:518841/1‑149 (MQ=255)
                                           gATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTcc                                                                                                       <  1:621447/149‑1 (MQ=255)
                                               tggcgttggcgCTGTTACTCGATGTGCAACAAGTTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGtt                                                                                                   <  1:315451/149‑1 (MQ=255)
                                               gggcgttggcgCTGTTACGCGATGTGCAACAAGGTCTGCGTGAGGACCTTAAACTGCTGATTATGGCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGtt                                                                                                   <  1:1038021/148‑1 (MQ=255)
                                                         gcTGTTACTCTATGTGCAACAAGGTCTGCGGGATGACCTTAAACTGCTGATTATGTCTGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGGCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGtt                                                                                         <  1:962763/149‑1 (MQ=255)
                                                                  cGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGc                                                                                >  2:996827/1‑149 (MQ=255)
                                                                     tgtgCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCGCGTTTCCGGTTGAACGCCGTTATTTACCGCTGc                                                                             <  2:9656/149‑1 (MQ=255)
                                                                       tgCTACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGccc                                                                           <  2:132982/149‑1 (MQ=255)
                                                                         caacaaGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGACCGCCGGTATTTACCGCTGCCcgc                                                                         >  2:159928/1‑149 (MQ=255)
                                                                         caacaaGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGAGCCTGTCGTCATCTCAGACGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCcgc                                                                         >  2:781025/1‑149 (MQ=255)
                                                                          aacaaGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCcgcg                                                                        >  1:573016/1‑149 (MQ=255)
                                                                               ggTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGCCCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCa                                                                   >  1:645560/1‑149 (MQ=255)
                                                                                 gcTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGc                                                                 <  2:304311/148‑1 (MQ=255)
                                                                                     cGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAAGGCCGTTATTTACCGCTGCCCGGGGAGCAGCGttt                                                             >  2:656634/1‑149 (MQ=255)
                                                                                      gggatgaCCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGtttt                                                            <  2:1135748/147‑1 (MQ=255)
                                                                                         atgaCCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGAc                                                         >  1:905124/1‑149 (MQ=255)
                                                                                              cTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGc                                                    >  1:882883/1‑149 (MQ=255)
                                                                                              cTTAAACTGCTGATTATGGCGGCTACCCTGGACAACGACCGGTTGCAGCAAATGCTGCAAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTTAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGc                                                    <  1:23983/149‑1 (MQ=255)
                                                                                                  aaCTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGtt                                                <  2:54841/149‑1 (MQ=255)
                                                                                                    ctgctgATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCCGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGc                                              <  2:783914/149‑1 (MQ=255)
                                                                                                       ctgATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGt                                           <  1:309936/149‑1 (MQ=255)
                                                                                                          aTTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGc                                        <  2:1053339/149‑1 (MQ=255)
                                                                                                           ttATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCATAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGcc                                       <  1:825560/149‑1 (MQ=255)
                                                                                                           ttATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGcc                                       >  1:1013561/1‑149 (MQ=255)
                                                                                                               gTCGGCTACCCTGTACAACGACCGCTGGCAGCAAATGCTGCCAGACGCGCCTGGCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGGCGTTATTTACCGCTGCCCGCGCATCTGCGTTTTGACGATGCCGTTGCGGTAGCCACCg                                   <  2:932859/149‑1 (MQ=255)
                                                                                                                     tACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGaaa                             >  2:173060/1‑149 (MQ=255)
                                                                                                                         cTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAAtgct                         >  2:346546/1‑149 (MQ=255)
                                                                                                                           ggACAACGACCGCTTGCAGCAAATGCGGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTAGCCGGGTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAAagctgc                       <  2:464136/149‑1 (MQ=255)
                                                                                                                            gACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCg                      >  1:998685/1‑149 (MQ=255)
                                                                                                                            gACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCGCGTTTCCGTTTGAACGCCGTTATTTACCGCTGCCCGCGCAGCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCg                      <  2:183527/149‑1 (MQ=255)
                                                                                                                            caccccGACCGCGTGCAGCACAGGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGGAATTTACCGCTGCCCGCGCATCAGCGTTTTGCCGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCg                      <  2:647637/144‑1 (MQ=255)
                                                                                                                              cAACGACCGCTTGCAGCAATTGCTTCCAGAAGCGCCTGTCGTCATCTCCGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTc                    <  1:556334/149‑1 (MQ=255)
                                                                                                                              cAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGGCGTCATCGCAGAAGGGCGCTCGTTTCCGGTTGAACTCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCAGTTGCGGTAGCCACCGCTGAAATGCTGCGTc                    >  2:640521/1‑149 (MQ=255)
                                                                                                                                  gACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTAATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTCAAATGCTGCCTCAGGa                >  1:837294/1‑149 (MQ=255)
                                                                                                                                    ccGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGaaa              >  1:438953/1‑149 (MQ=255)
                                                                                                                                    ccGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGCCATCTCAGAAGGGCGCTCGTTTCCGGATGAACGCCGTTATTTACCGCTGCCCGCGCATCAGGCTTTTGACGATGCCGTTGCGGAAGCCCCCGCTGAAATGCTGCGTCAGGaaa              >  2:789976/1‑149 (MQ=255)
                                                                                                                                     cGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGATAg             >  1:601663/1‑149 (MQ=255)
                                                                                                                                      gCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGAAAGc            >  2:3527/1‑149 (MQ=255)
                                                                                                                                       ctggcagcaTTTGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGAAAGCg           >  1:596602/4‑149 (MQ=255)
                                                                                                                                       cTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGAAAGCg           >  2:563727/1‑149 (MQ=255)
                                                                                                                                        tttcagcaAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGATATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGAAAGCgg          >  2:1143979/4‑149 (MQ=255)
                                                                                                                                         ggcagcaTTTGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGAAAGCGGa         >  2:223136/2‑149 (MQ=255)
                                                                                                                                          gcagcaAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGAAAGCGGAt        <  1:372140/149‑1 (MQ=255)
                                                                                                                                           cagcaAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGAAAGCGGATc       >  1:557722/1‑149 (MQ=255)
                                                                                                                                             gcaAATGCTGCCTGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGAAAGCGGATCAt     >  2:105693/1‑149 (MQ=255)
                                                                                                                                              caAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCTTCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGAATGCGGATCAtt    >  1:1048995/1‑149 (MQ=255)
                                                                                                                                                aaTGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGAAAGCGGATCATTAc  >  2:232633/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
ACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATTATGTCGGCTACCCTGGACAACGACCGCTTGCAGCAAATGCTGCCAGAAGCGCCTGTCGTCATCTCAGAAGGGCGCTCGTTTCCGGTTGAACGCCGTTATTTACCGCTGCCCGCGCATCAGCGTTTTGACGATGCCGTTGCGGTAGCCACCGCTGAAATGCTGCGTCAGGAAAGCGGATCATTAC  >  NC_000913/162458‑162750

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: ‑