Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 162,349 A→C Y82S (TAC→TCC)  hrpB → putative ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913162,3490AC96.2% 141.0 / ‑5.9 52Y82S (TAC→TCC) hrpBputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base C (19/31);  total (20/32)
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

GTCTTACCTGAATTACTTACCGCTCTCGATTGTGCGCCGCAGGTATTATTAAGTGCGCCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGT‑CGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTACCGGATGCGTGCG‑CAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATT  >  NC_000913/162129‑162566
                                                                                                                                                                                                                             |                                                                                                                                                                                                                          
gttttACCTGAATTGCGTACCGCTCTCGATTGTGCGCCGCAGGTATTATTAAGTGCGCCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGT‑CGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCG‑Ca                                                                                                                                                                                                            <  2:763179/233‑1 (MQ=255)
        tGAATTACTTACCGCTCTCGATTGTGCGCCGCAGGTATTATTAAGTGCGCCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGT‑CGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCG‑CAAAACTGCg                                                                                                                                                                                                    >  1:328276/1‑236 (MQ=255)
        tGAATTACTTACCGCTCTCGATTGTGCGCCGCAGGTATTATTAAGTGCGCCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGT‑CGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGGTGGGTCCCGGATGCGTGCG‑CAAAACTTCg                                                                                                                                                                                                    >  1:321661/1‑236 (MQ=255)
           aTTACTTACCGCTCTCGATTGTGCGCCGCAGGTATTATTAAGTGCGCCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGT‑CGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCG‑CAAAACTGCGTc                                                                                                                                                                                                  >  1:400970/1‑235 (MQ=255)
               cttacCGCTCTCGATTGTGCGCCGCAGGTATTATTAAGTGCGCCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGT‑CGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGGGCG‑CAAAACTGCGCCGGGcc                                                                                                                                                                                             >  1:199119/1‑236 (MQ=255)
                     ggtatggATTGTGCCCCGCAGTTATTATTAAGTCCGCAGACCGGGCCCGGGAAATCACCCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGT‑CGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCG‑CAAAACTGCGTCGGGCCGAATAc                                                                                                                                                                                       <  1:568428/230‑1 (MQ=255)
                        ctcAATTGGGCCCCGAAGGTATTATGAAGTGCGCCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGT‑CGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCG‑CAAAACTGCGTCGGGCCGAATACCCg                                                                                                                                                                                    <  1:411579/236‑1 (MQ=255)
                              ggtggcccgcAGGTATTATTACGTGCGCCGACCGGGGCCGGGAAATCACCCGGGCTGCCGCTGCAACGGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGT‑CGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCG‑CAAAACTGCGTCGGGCCGAATACCCGCCTGGa                                                                                                                                                                              <  1:151722/230‑1 (MQ=255)
                                  tcgcgagagttttagtaAGCGCGCCGACGGGGCCCGGAAACTCAACCTGGCTGCCCGTGCAACTGCTGGCGCATCCCGGCATTAGCGGGAAAATTATCCTGCTGGAGCCGCGT‑CGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCG‑CAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTg                                                                                                                                                                          <  2:604495/225‑1 (MQ=255)
                                            attattAGGTCCGCCGACCGGGGCCGGGAAATCAACCGGGCTCCCGTTGCAACTGCTGGCCCATCCCGGAATTACCGGGAAAATTACCCTGCTGGAGCCGCGT‑CGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCG‑CAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAg                                                                                                                                                                <  2:199119/236‑1 (MQ=255)
                                                     tGCGCCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGT‑CGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCG‑CAAAACTGCGGCGGGCCGAATACCCGCCTGGAAGTGGTGACCGAAGGCGTGCTGa                                                                                                                                                       >  1:209162/1‑236 (MQ=255)
                                                       cgcCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGT‑CGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCG‑CAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGAcg                                                                                                                                                     >  1:577105/1‑236 (MQ=255)
                                                       cgcCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGT‑CGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCG‑CAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGGGGTTACCGAAGGCGTGCTGAcg                                                                                                                                                     >  1:392032/1‑236 (MQ=255)
                                                       cgcCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGT‑CGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGAGGCGTGGG‑GAAAACGGCGACGGGCCGAATAACCGACTTGAAGTGGTTGCCGCAGgcgcgcgtacc                                                                                                                                                     >  2:96009/1‑231 (MQ=255)
                                                        gcCGACCGGGGCCGGAAAATCACCCTGGCTGCCGCTGAACCTGCTGGCCCTTCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGT‑CGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCG‑CAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGAc                                                                                                                                                      <  2:305684/234‑1 (MQ=255)
                                                           gcccGGTGCCGAGAATTCATCGTGCCTGCCCGTGCAACGGCGGGCCCTTCCAGGCATTAACGGGAAAATTATCCTGCTGGGGCCGCGT‑CGTCGGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCG‑CAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACgcgc                                                                                                                                                  <  2:737770/233‑1 (MQ=255)
                                                           gACCGGGGCCGGGAAACCAACCGGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGT‑CGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCG‑CAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCa                                                                                                                                                 <  1:96009/236‑1 (MQ=255)
                                                            aCCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCCACTGCTGGCGCCTCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGT‑CGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGGGCGGCTTAACGAAAAGCCAGGCGCTACCGTTGGCCACCGGATGGAGGCGTCAAAAC‑GCGTCGGTGCGACTACCCGCCTCGAAGCGGGTACCGGACGCGTGCTTACgcgcct                                                                                                                                                >  2:443600/1‑234 (MQ=255)
                                                                      ggAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGT‑CGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATAACGTTGGCCCCCGGATGCGTGCG‑CAAAAATGCGTCGGGCCGCATACCCGCCTGGAAGTGGTTACCCAATGCGGGCTGAGGCGCACGAACCAGCg                                                                                                                                       >  2:316500/1‑235 (MQ=255)
                                                                          aTCACCCTGGCTGCCGCTGCACCTGGTGGCCAACCCCGGCATTACGGGGAAAATTACCCTGCGGGAGCCGGGT‑CGTCTGGCGGCGCGAAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCG‑CAAAACTGCGTCGGGCCGAATCCCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCg                                                                                                                                       <  2:321661/231‑1 (MQ=255)
                                                                                   gCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGT‑CGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCG‑CCAAACTGCGTCGGGCCGAATACCCGCCTGGGAGTGGGTACCGAAGGCGTGCTGAGCGGCATGATCCAGCGTGACCCGGAACTGa                                                                                                                         >  2:92643/1‑236 (MQ=255)
                                                                                   gCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGT‑CGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCG‑CAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGGGCATGATCCACGCTGACCCgggaattg                                                                                                                         >  1:739632/1‑230 (MQ=255)
                                                                                    ttcccGCTGAAACTGCTGGCGCATCCCGGAATTACCGGGAAAATTATCCTGCGGGGGCCGCGT‑CGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCG‑CAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAg                                                                                                                        <  2:134791/233‑1 (MQ=255)
                                                                                    ctgctgCGGCATGTGCTGGCCCACACGGGCGTTAACGTGAAAATTACCCTGCTGGAGCCGCGT‑CGTCTGGGGGAGCGTACCGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAGGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCG‑CAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAg                                                                                                                        <  2:741196/236‑1 (MQ=255)
                                                                                    cTGCCGGTGCAACTGCTGGCTCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGT‑CGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCG‑CAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAg                                                                                                                        <  1:217379/236‑1 (MQ=255)
                                                                                    cTGCCGCTGCACCTGCGGGCGCATCCCGGCATTAACGGGAAAATTAGCCGGCTGGGGCCGGGT‑CGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCG‑CAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAg                                                                                                                        <  1:379045/236‑1 (MQ=255)
                                                                                    cTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGGGCCGCGT‑CGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCG‑CAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAg                                                                                                                        <  1:316500/236‑1 (MQ=255)
                                                                                    cTGCCGCTGAAACGGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGT‑CGTATGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCG‑CAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAg                                                                                                                        <  1:206377/236‑1 (MQ=255)
                                                                                       ccGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGT‑CGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCG‑CAAAACTGCGTCGGGCCGAATAGCCGCCTGGAAGTGGTTTCCGAAGGCCGGCTGCAGCGGCGGATCCAGCGGGAGCCGGagccgacgag                                                                                                                     >  2:460295/1‑227 (MQ=255)
                                                                                                       cccATCCCGGCATTAACGGGACAATTATCCTGCTGGACCCGCGT‑CGTCTGGCGGCGCGTAACGTCGCGAAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCG‑CAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGt                                                                                                         <  2:647060/230‑1 (MQ=255)
                                                                                                            ccgggtaTTACCGGAAAAAGTATCCGGCGGGGCCCGTGT‑CGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCG‑CAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTtga                                                                                                <  2:392032/230‑1 (MQ=255)
                                                                                                            cccGGCATTACCGGGAAAATTCTCCTGCGGGAGCC‑CGTGTGGCGGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCG‑CAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTtga                                                                                                <  2:577105/236‑1 (MQ=255)
                                                                                                                         ggAAAATTATCCTGCTGGAGCCGCGT‑CGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCG‑CAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGAt                                                                                                                                             >  2:627086/1‑178 (MQ=255)
                                                                                                                         ggAAAATTATCCTGCTGGAGCCGCGT‑CGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCG‑CAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGAt                                                                                                                                             <  1:627086/178‑1 (MQ=255)
                                                                                                                                       ctgGAGCCGCGT‑CGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCG‑CAAAACTGCGTCGGGCCGAATAcc                                                                                                                                                                                      <  1:24765/123‑1 (MQ=255)
                                                                                                                                       ctgGAGCCGCGT‑CGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCG‑CAAAACTGCGTCGGGCCGAATAcc                                                                                                                                                                                      >  2:24765/1‑123 (MQ=255)
                                                                                                                                         ggAGCCGCGT‑CGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCG‑CAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGGGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCgg                                                                   >  1:439873/1‑236 (MQ=255)
                                                                                                                                                cgt‑cgtCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCG‑CAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAg                                                                                                                        <  2:346698/176‑1 (MQ=255)
                                                                                                                                                cgt‑cgtCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCG‑CAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAg                                                                                                                        >  1:346698/1‑176 (MQ=255)
                                                                                                                                                         ggcggcGCGTAACGTCGCGAAACGGCTGCCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCG‑CAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCg                                                                    <  1:581960/220‑1 (MQ=255)
                                                                                                                                                              cgcgTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCG‑CAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTa                                               >  1:340597/1‑236 (MQ=255)
                                                                                                                                                                    aCGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCG‑CAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGTCGCTGTTACTcca                                          >  1:666557/1‑233 (MQ=255)
                                                                                                                                                                           gcAACGGCTGGCGGAGCTTCTTAACGAAAAGCCGGGCGATACCGTTGGCCCCCGGGTGCGTGCC‑CAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGcaac                                  <  2:674564/236‑1 (MQ=255)
                                                                                                                                                                            cAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCG‑CAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGcaac                                  <  1:430567/235‑1 (MQ=255)
                                                                                                                                                                            cAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGCCTCCCGGATGCGGGCG‑CAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGcaaca                                 <  1:385805/236‑1 (MQ=255)
                                                                                                                                                                                 gCTGGCGGACCTGCTTAACGAAAAGCCAGGCGATTCCGTTGGCTCCCGGATGCGTGCG‑CAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTc                            <  1:386522/236‑1 (MQ=255)
                                                                                                                                                                                         aGCTGCTTAACGAAAGGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCG‑CAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGtga                     <  1:92643/235‑1 (MQ=255)
                                                                                                                                                                                          gtcgctAAACGAAACGGCAGGCGATACCGTTGGCTCCCGGATGGGTGCC‑AAAACCGGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGtgatg                   <  2:177534/233‑1 (MQ=255)
                                                                                                                                                                                                  aCGAAAAGCCAGGCTATACCGTTGGCTACCGGATGGCTGCG‑CAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTaa            <  1:742401/235‑1 (MQ=255)
                                                                                                                                                                                                   cgaagcgCCAGGCAGTACCGTTGGCTCCGGGATGGGGGCG‑GAAAACTGGTTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAAc          <  1:83852/230‑1 (MQ=255)
                                                                                                                                                                                                      aaaGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCG‑CAAAACTGCGTCGGGCCGAACACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTTTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAActgc       >  1:545041/1‑236 (MQ=255)
                                                                                                                                                                                                             ggtgATACCGTTGGCTCTGGGAGGTGTGCG‑CAACTCTGGGTGGGCCCGAATTCCCGGCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGAtt  <  2:400970/231‑1 (MQ=255)
                                                                                                                                                                                                                             |                                                                                                                                                                                                                          
GTCTTACCTGAATTACTTACCGCTCTCGATTGTGCGCCGCAGGTATTATTAAGTGCGCCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGT‑CGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTACCGGATGCGTGCG‑CAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATT  >  NC_000913/162129‑162566

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑