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

GTCTTACCTGAATTACTTACCGCTCTCGATTGTGCGCCGCAGGTATTATTAAGTGCGCCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATT  >  NC_000913/162129‑162566
                                                                                                                                                                                                                             |                                                                                                                                                                                                                          
gttttACCTGAATTGCGTACCGCTCTCGATTGTGCGCCGCAGGTATTATTAAGTGCGCCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCGCa                                                                                                                                                                                                            <  2:763179/233‑1 (MQ=255)
        tGAATTACTTACCGCTCTCGATTGTGCGCCGCAGGTATTATTAAGTGCGCCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCGCAAAACTGCg                                                                                                                                                                                                    >  1:328276/1‑236 (MQ=255)
        tGAATTACTTACCGCTCTCGATTGTGCGCCGCAGGTATTATTAAGTGCGCCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGGTGGGTCCCGGATGCGTGCGCAAAACTTCg                                                                                                                                                                                                    >  1:321661/1‑236 (MQ=255)
           aTTACTTACCGCTCTCGATTGTGCGCCGCAGGTATTATTAAGTGCGCCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCGCAAAACTGCGTc                                                                                                                                                                                                  >  1:400970/1‑235 (MQ=255)
               cttacCGCTCTCGATTGTGCGCCGCAGGTATTATTAAGTGCGCCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGGGCGCAAAACTGCGCCGGGcc                                                                                                                                                                                             >  1:199119/1‑236 (MQ=255)
                     ggtatggATTGTGCCCCGCAGTTATTATTAAGTCCGCAGACCGGGCCCGGGAAATCACCCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATAc                                                                                                                                                                                       <  1:568428/230‑1 (MQ=255)
                        ctcAATTGGGCCCCGAAGGTATTATGAAGTGCGCCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCg                                                                                                                                                                                    <  1:411579/236‑1 (MQ=255)
                              ggtggcccgcAGGTATTATTACGTGCGCCGACCGGGGCCGGGAAATCACCCGGGCTGCCGCTGCAACGGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGa                                                                                                                                                                              <  1:151722/230‑1 (MQ=255)
                                  tcgcgagagttttagtaAGCGCGCCGACGGGGCCCGGAAACTCAACCTGGCTGCCCGTGCAACTGCTGGCGCATCCCGGCATTAGCGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTg                                                                                                                                                                          <  2:604495/225‑1 (MQ=255)
                                            attattAGGTCCGCCGACCGGGGCCGGGAAATCAACCGGGCTCCCGTTGCAACTGCTGGCCCATCCCGGAATTACCGGGAAAATTACCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAg                                                                                                                                                                <  2:199119/236‑1 (MQ=255)
                                                     tGCGCCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCGCAAAACTGCGGCGGGCCGAATACCCGCCTGGAAGTGGTGACCGAAGGCGTGCTGa                                                                                                                                                       >  1:209162/1‑236 (MQ=255)
                                                       cgcCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGAcg                                                                                                                                                     >  1:577105/1‑236 (MQ=255)
                                                       cgcCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGGGGTTACCGAAGGCGTGCTGAcg                                                                                                                                                     >  1:392032/1‑236 (MQ=255)
                                                       cgcCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGAGGCGTGGGGAAAACGGCGACGGGCCGAATAACCGACTTGAAGTGGTTGCCGCAGgcgcgcgtacc                                                                                                                                                     >  2:96009/1‑231 (MQ=255)
                                                        gcCGACCGGGGCCGGAAAATCACCCTGGCTGCCGCTGAACCTGCTGGCCCTTCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGAc                                                                                                                                                      <  2:305684/234‑1 (MQ=255)
                                                           gcccGGTGCCGAGAATTCATCGTGCCTGCCCGTGCAACGGCGGGCCCTTCCAGGCATTAACGGGAAAATTATCCTGCTGGGGCCGCGTCGTCGGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACgcgc                                                                                                                                                  <  2:737770/233‑1 (MQ=255)
                                                           gACCGGGGCCGGGAAACCAACCGGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCa                                                                                                                                                 <  1:96009/236‑1 (MQ=255)
                                                            aCCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCCACTGCTGGCGCCTCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGGGCGGCTTAACGAAAAGCCAGGCGCTACCGTTGGCCACCGGATGGAGGCGTCAAAACGCGTCGGTGCGACTACCCGCCTCGAAGCGGGTACCGGACGCGTGCTTACgcgcct                                                                                                                                                >  2:443600/1‑234 (MQ=255)
                                                                      ggAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATAACGTTGGCCCCCGGATGCGTGCGCAAAAATGCGTCGGGCCGCATACCCGCCTGGAAGTGGTTACCCAATGCGGGCTGAGGCGCACGAACCAGCg                                                                                                                                       >  2:316500/1‑235 (MQ=255)
                                                                          aTCACCCTGGCTGCCGCTGCACCTGGTGGCCAACCCCGGCATTACGGGGAAAATTACCCTGCGGGAGCCGGGTCGTCTGGCGGCGCGAAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATCCCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCg                                                                                                                                       <  2:321661/231‑1 (MQ=255)
                                                                                   gCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCGCCAAACTGCGTCGGGCCGAATACCCGCCTGGGAGTGGGTACCGAAGGCGTGCTGAGCGGCATGATCCAGCGTGACCCGGAACTGa                                                                                                                         >  2:92643/1‑236 (MQ=255)
                                                                                   gCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGGGCATGATCCACGCTGACCCgggaattg                                                                                                                         >  1:739632/1‑230 (MQ=255)
                                                                                    ttcccGCTGAAACTGCTGGCGCATCCCGGAATTACCGGGAAAATTATCCTGCGGGGGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAg                                                                                                                        <  2:134791/233‑1 (MQ=255)
                                                                                    ctgctgCGGCATGTGCTGGCCCACACGGGCGTTAACGTGAAAATTACCCTGCTGGAGCCGCGTCGTCTGGGGGAGCGTACCGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAGGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAg                                                                                                                        <  2:741196/236‑1 (MQ=255)
                                                                                    cTGCCGGTGCAACTGCTGGCTCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAg                                                                                                                        <  1:217379/236‑1 (MQ=255)
                                                                                    cTGCCGCTGCACCTGCGGGCGCATCCCGGCATTAACGGGAAAATTAGCCGGCTGGGGCCGGGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAg                                                                                                                        <  1:379045/236‑1 (MQ=255)
                                                                                    cTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGGGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAg                                                                                                                        <  1:316500/236‑1 (MQ=255)
                                                                                    cTGCCGCTGAAACGGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTATGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAg                                                                                                                        <  1:206377/236‑1 (MQ=255)
                                                                                       ccGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATAGCCGCCTGGAAGTGGTTTCCGAAGGCCGGCTGCAGCGGCGGATCCAGCGGGAGCCGGagccgacgag                                                                                                                     >  2:460295/1‑227 (MQ=255)
                                                                                                       cccATCCCGGCATTAACGGGACAATTATCCTGCTGGACCCGCGTCGTCTGGCGGCGCGTAACGTCGCGAAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGt                                                                                                         <  2:647060/230‑1 (MQ=255)
                                                                                                            ccgggtaTTACCGGAAAAAGTATCCGGCGGGGCCCGTGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTtga                                                                                                <  2:392032/230‑1 (MQ=255)
                                                                                                            cccGGCATTACCGGGAAAATTCTCCTGCGGGAGCCCGTGTGGCGGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTtga                                                                                                <  2:577105/236‑1 (MQ=255)
                                                                                                                         ggAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGAt                                                                                                                                             >  2:627086/1‑178 (MQ=255)
                                                                                                                         ggAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGAt                                                                                                                                             <  1:627086/178‑1 (MQ=255)
                                                                                                                                       ctgGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATAcc                                                                                                                                                                                      <  1:24765/123‑1 (MQ=255)
                                                                                                                                       ctgGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATAcc                                                                                                                                                                                      >  2:24765/1‑123 (MQ=255)
                                                                                                                                         ggAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGGGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCgg                                                                   >  1:439873/1‑236 (MQ=255)
                                                                                                                                                cgtcgtCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAg                                                                                                                        <  2:346698/176‑1 (MQ=255)
                                                                                                                                                cgtcgtCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAg                                                                                                                        >  1:346698/1‑176 (MQ=255)
                                                                                                                                                         ggcggcGCGTAACGTCGCGAAACGGCTGCCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCg                                                                    <  1:581960/220‑1 (MQ=255)
                                                                                                                                                              cgcgTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTa                                               >  1:340597/1‑236 (MQ=255)
                                                                                                                                                                    aCGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGTCGCTGTTACTcca                                          >  1:666557/1‑233 (MQ=255)
                                                                                                                                                                           gcAACGGCTGGCGGAGCTTCTTAACGAAAAGCCGGGCGATACCGTTGGCCCCCGGGTGCGTGCCCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGcaac                                  <  2:674564/236‑1 (MQ=255)
                                                                                                                                                                            cAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGcaac                                  <  1:430567/235‑1 (MQ=255)
                                                                                                                                                                            cAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGCCTCCCGGATGCGGGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGcaaca                                 <  1:385805/236‑1 (MQ=255)
                                                                                                                                                                                 gCTGGCGGACCTGCTTAACGAAAAGCCAGGCGATTCCGTTGGCTCCCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTc                            <  1:386522/236‑1 (MQ=255)
                                                                                                                                                                                         aGCTGCTTAACGAAAGGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGtga                     <  1:92643/235‑1 (MQ=255)
                                                                                                                                                                                          gtcgctAAACGAAACGGCAGGCGATACCGTTGGCTCCCGGATGGGTGCCAAAACCGGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGtgatg                   <  2:177534/233‑1 (MQ=255)
                                                                                                                                                                                                  aCGAAAAGCCAGGCTATACCGTTGGCTACCGGATGGCTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTaa            <  1:742401/235‑1 (MQ=255)
                                                                                                                                                                                                   cgaagcgCCAGGCAGTACCGTTGGCTCCGGGATGGGGGCGGAAAACTGGTTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAAc          <  1:83852/230‑1 (MQ=255)
                                                                                                                                                                                                      aaaGCCAGGCGATACCGTTGGCTCCCGGATGCGTGCGCAAAACTGCGTCGGGCCGAACACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTTTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAActgc       >  1:545041/1‑236 (MQ=255)
                                                                                                                                                                                                             ggtgATACCGTTGGCTCTGGGAGGTGTGCGCAACTCTGGGTGGGCCCGAATTCCCGGCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGAtt  <  2:400970/231‑1 (MQ=255)
                                                                                                                                                                                                                             |                                                                                                                                                                                                                          
GTCTTACCTGAATTACTTACCGCTCTCGATTGTGCGCCGCAGGTATTATTAAGTGCGCCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAACGAAAAGCCAGGCGATACCGTTGGCTACCGGATGCGTGCGCAAAACTGCGTCGGGCCGAATACCCGCCTGGAAGTGGTTACCGAAGGCGTGCTGACGCGCATGATCCAGCGTGACCCGGAACTGAGCGGTGTTGGACTGGTGATCCTTGATGAATTTCATGAGCGCAGCTTGCAGGCGGATTTGGCGTTGGCGCTGTTACTCGATGTGCAACAAGGTCTGCGTGATGACCTTAAACTGCTGATT  >  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: