Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,401,607 C→T T312T (ACC→ACT hflX → ribosome rescue factor HflX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,401,6070CT98.4% 173.5 / NA 63T312T (ACC→ACThflXribosome rescue factor HflX
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (25/36);  minor base G (1/0);  total (27/36)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

ACCCTCGACCCGACGTTGCGGCGTATTGACGTTGCAGATGTCGGTGAAACCGTACTTGCAGATACCGTAGGGTTTATTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGTACAAGAAAACATCGAAGCGGTGAATACGGTTCTTGAAGAGATCGACGCTCACGAGATCCCAACCCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCGCGTATTGATCGGGACGAAGAGAACAAACCGAACCGTGTCTGGCTTTCCGCACAGACCGGAGCGGGGATACCACAGCTTTTTCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGT  >  NC_000913/4401362‑4401848
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 
ccccTCGACCCGACTTTGCGGCGTATTGCCGTTGCAGATGTCGGTGAAACCGTACTTGCAGATACCGTAGGGTTTATTCGCCACCGGCCGCACGATCTGGTGGCTGCATTTAAAGCCACGTTACAAGAGTCGCGGCAAGCCACATTACTGCTGCACGGCATTGATGCGGCTGATGTGCGTGTACAAGAAAACATCGAAGCGGGGAATACGGTTCTTGAAGAGATCGACGCTCACGAGATCCCAACTctgct                                                                                                                                                                                                                                              <  1:382586/250‑1 (MQ=255)
 cccTCGACCCGACGTTGCGGCGTATTGACGTTGCAGATGTCGGTGAAACCGTACTTGCAGATACCGTAGGGTTTATTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGTACCAGGAAAAATCGAAGCGGTGAATAAGGATCTTGAAAAGATCGACGCTCACGAGATCCCAACTctgctg                                                                                                                                                                                                                                             >  1:264755/1‑251 (MQ=255)
  ccgcGACCCGGCTTTGCGGCGTATTGACGTTGCAGAGGTCGGTGAAACCGTACTTGCAGATACCGTAGTGTTTATTTGCCTCCTGCCGCACGATCTGTTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGTACAAGAAAACATCGAAGCGGTGAAGACGGTTCTTGAAGAGATCGACGCTCACGAGATCCCAACTCTGCTgg                                                                                                                                                                                                                                            <  2:375737/248‑1 (MQ=255)
  ccTCGACCCGCAGTTGCGGCGTATTGACGTTGCAGATGTCGGTGAAACCGTACTTGCAGCTACCGTAGGGATTATTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGTACAAGAAAAAATCGCAGCGGTGAATACGGTTCTTGCAGTGATCGACGCTAACGAGATCCCAACTCTGCTgg                                                                                                                                                                                                                                            >  2:53647/1‑251 (MQ=255)
                  gggcgTGTTGACGTTTCAGAGGTCGGTGAAAACGTACTTGCAGATCCCGTAGGGTTTTTTCGCCACTGGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACAGTACTGCTGCACGTCATTGATGCGGCGGAGGTGCGTGTACAAGAAAACATCGAAGCGGTGAATACGGGTCTTGAAGAGATCGACGCTCACGAGATCCTAACTCGGCTGTTGATGAACACGATCGa                                                                                                                                                                                                                            <  2:272167/250‑1 (MQ=255)
                   ggcgTATTGACGTTGCATATGTCGGCGAAACCGTACTTGCAGATCCCGTAGTTGTTATTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTTCAAGAGACGCGACAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGTACAAGAAAACATCGAAGCGGTGAATACGGTTCTTGAAGACATCGACGCTCACGAGATCCCAACTCTGCTGGTGATGAACATGATCGat                                                                                                                                                                                                                           <  2:32043/251‑1 (MQ=255)
                    gcgTATTGACGTTGCAGATGTCGGTGAAACCGTACTTGCAGATACCGTAGGGTTTATTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGTACAAGAAAACATCGAAGCGGTGAATACGGTTCTTGAAGAGATCGACGCTCACGAGATCCCAACTCTGCTGGTGATGAACAAGATCGata                                                                                                                                                                                                                          <  1:369779/251‑1 (MQ=255)
                     cgTATTGACGTTGCAGATGTCGGTGAAACCGTACTTGCAGATACCGTAGGGTTTATTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGTACAAGAAAACATCGAAGCGGTGAATACGGTTCTTGAAGAGATCGACGCTCACGAGATCCCAACTCTGCTGGTGATGAACAAGATCGatat                                                                                                                                                                                                                         <  1:174305/251‑1 (MQ=255)
                       tATTGACGTTGCAGATGTCGGTGAAACCGTACTTGCAGATACCGTAGGGTTTATTCGCCACCTGCAGCACGATCTGGAGGCGGCATTTAAAGCAACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGTACAAGAAAACCTAGCAGCGGTGAATACGGTTCCTGAAGAGATCGACGCTCACACGATCCCCACTCTGCTGGTGATGAACAAGATCGATATgt                                                                                                                                                                                                                       >  2:170938/1‑250 (MQ=255)
                        aTTGACGTTGCAGATGTCGGTGAAACCGTACTTGCAGATACCGTAGGGTTTATTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGTACAAGAAAACATCGAAGCGGTGAATACGGTTCTTTAAGAGATCGACGCTCACTCGATCCCAACTCCGCTGGTGATGAACAAGATCGATATGCt                                                                                                                                                                                                                      >  1:319910/1‑251 (MQ=255)
                           gACGTTGCAGATGTCGGTGAAACCGTACTTGCAGATTCCGTAGGGTTTATTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGTACAAGAAAACATCGAAGCGGTGAATACGGTTCTTGAAGAGATCGACGCTCACGAGATCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGa                                                                                                                                                                                                                   <  2:323176/251‑1 (MQ=255)
                             cGTTGCAGATGTCGGTGAAACCGTACTTGCAGATACCGTAGGGTTTATTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGTACAAGAAAACATCGAAGCGGTGAATACGGTTCTTGAAGAGATCGACGCTCACGAGATCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAAg                                                                                                                                                                                                                 <  1:115708/251‑1 (MQ=255)
                                 tcAGATGACGGTGAAAACGTACTTGCAGATAACGTAGGGTTTATTCGCCACCTGCCGCACGATCTGGTGGCGGCCTTTAAAGCCACGTTCAAAGAGACGCGGCACGCAACATTACTGCTGCACGTCATTGATGCGGCGGAGGTGATTGTACAAAAAAACATCGAAGCGGTGAAAACGGGTCTTGAAGAGATTGGCGAGAACGAGATCCCAACCCTGCCGGTGATGGACAAGaactcaatacaggacgattt                                                                                                                                                                                                             >  2:232341/2‑232 (MQ=255)
                                   aGATGTCGGTGAAACCGTACTTGCAGATACCGTAGGGTTTATTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGTACAAGAAAACATCGAAGCGGTGAATACGGTTCTTGAAGAGATCGACGCTCACGAGATCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCg                                                                                                                                                                                                           <  2:100144/251‑1 (MQ=255)
                                   aGATGTCGGTGAAACCGTACTTGCAGATACCGTAGGGTTTATTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGGGTACAAGAAAACATCGAAGCGGTGAAAACGGGTCTTGAAGAGATCGAAGCTAACGAGATCCCAACTCTTATGGGGATGAACAAgcacatgctggtggaagatttct                                                                                                                                                                                                           >  2:312660/1‑229 (MQ=255)
                                           gTGAAACCGTACTTGCAGATACCGTAGGGTTTATTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGCTGCGCGTGTACAAGAAAACATCGAAGCGGTGAATACGGTTCTTGAAGGGATCGACGCTCACGCGCTCCCACCTCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCGCGt                                                                                                                                                                                                   >  1:59559/1‑251 (MQ=255)
                                             gaatccGTACTTGCAGATACCGTAGGGTTTATTCGCCCCCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCGCGTCATTGATGCGGCGTATGTGCGTGTACAAGTAAACATCGAAGCGGTGAATACGGTTCTTGAAGAGATCGACGCTCACGAGATCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCGCGTAt                                                                                                                                                                                                 <  1:203473/247‑1 (MQ=255)
                                               aaCCGTACTTGCAGATACCGTAGGGTTTATTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGTACAAGAAAACATCGAAGCGGTGAATACGGTTCTTGAAGAGATCGACGCTCACGGGATCCCAACTCTGCTGGGGATGAAAAAGATCGATATGCTGGAAGATTTCGAACCGCGTATTg                                                                                                                                                                                               >  1:127524/1‑251 (MQ=255)
                                                 ccGTACTTGCAGATACCGTAGGGTTTATTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGTACAAGAAAACATCGAAGCGGTGAATACGGTTCTTGAAGAGATCGACGCTCACGAGATCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCGCGTATTGAt                                                                                                                                                                                             <  1:77321/251‑1 (MQ=255)
                                                   gTACTTGCAGATACCGTAGGGTTTATTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGTACAAGAAAACATCGAAGCGGTGAATACGGTTCTTGAAGAGATCGACGCTCACGAGATCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCGCGTATTGATCg                                                                                                                                                                                           >  1:192382/1‑251 (MQ=255)
                                                        ggCAGATACCGTAGGGTTTAGTCGCCACCTGCCTCACGATCGGGGGGCGGCATGTAAAGCCACGTGACAAGAGACGCGGCAAGCCACATTACTGCTGCACGGCATTGATGCGGCGGATGTGCGTGTACAAGAAAACATCGAAGCGGTGAATACGGTTCTTGAAGAGAACGACGCTCACGAGATCCCAACTCTGCTGGTGATGAACAAGATCTATATGCTGGAAGAGTTCGAACCGCGTATTGATCGGGaa                                                                                                                                                                                       <  2:296829/249‑2 (MQ=255)
                                                                   tagggaatATTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCAACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGCTGTGCGTGTACCAGAAAACATCGAAGCGGTGAATACGGTTCTTGAAGGGATCAACGCTCACGAGATCCCCACTCTGCTGGTGATGAAAAAGATAGATATGCTGGAAGATTTCGAACCGCGTATTGATCGGGACGAAGAGACCaa                                                                                                                                                                            >  1:261982/8‑250 (MQ=255)
                                                                         ttATTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGTACAAGAAAACATCGAAGCGGTGAATACGGTTCTTGAAGAGATCGACGCTCACGAGATCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCGCGTATTGATCGGGACGAAGAGAACAaaccgaa                                                                                                                                                                      <  2:127524/250‑1 (MQ=255)
                                                                                 cACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGTACAAGAAAACATCGAAGCGGTGAATACGGTTTTTGAAGAGATCGACGCTCACGAGATCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCGCGTATTGATCGGGACGAAGAGAACAAACCGAACCGTGTCTg                                                                                                                                                             <  1:353406/251‑1 (MQ=255)
                                                                                      gCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGGGCGTGTACCAGAAAACATCGAAGCGGGGAATAAGGTTCTTGAAGAGATCGACGGTCACGAGATCCCAACTCTGATGGTGATGAACAAGATCGATATGCTGGGAGATTTCAAACCGCGTATTGATCGGGACTAAAAGAACAAACCGAACCGGGTCTGGCttt                                                                                                                                                        >  2:293771/1‑251 (MQ=255)
                                                                                                    tggCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGTACAAGAAAACATCGAAGCGGTGAATACGGTTCTTGAAGAGATCGACGCTCACGAGATCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCGCGTATTGATCGGGACGAAGAGAACAAACCGAACCGTGTCTGGCTTTCCGCACAGACCggg                                                                                                                                          >  1:53139/1‑250 (MQ=255)
                                                                                                          cATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGTACAAGAAAACATCGAAGCGGTGAATACGGTTCTTGAAGAGATCGACGCTCACGAGATCCCAACTCTGCTGGTGATGAACAAGATCGATAGGCTGGAAGATTTCGAACCGCGTATTGATCGGGACGAAGAGAACAAACCGAACCGTGTCTGGCTTTCCGCACAGACCGGAGCgggg                                                                                                                                    <  2:261039/251‑1 (MQ=255)
                                                                                                                   ccACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGTACAAGAAAACATCGAAGCGGTGAATACGGTTCTTGAAGAGATCGACGCTCACGAGATCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCTCGTATTGATCGGGACGAAGAGAACAAACCGAACCGTGTCTGGCTTTCCGCACAGACCGGAGCGGGGATACCACAg                                                                                                                           <  1:65195/251‑1 (MQ=255)
                                                                                                                     aCGTTACAAGAGACGGGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGTACAAGAAAACATCGAAGCGGTGAATACGGTTCTTGAAGAGATCGACGCTCACGAGATCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCGCGTATTGATCGGGACGAAGAGAACAAACCGAACCGTGTCTGGCTTTCCGCACAGACCGGAGCGGGGATACCACAGCt                                                                                                                         <  2:312457/251‑1 (MQ=255)
                                                                                                                       gTTACAAGAGACGCGGCAAGCCACAATACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGTACAAGAAAACATCGAAGCGGTGAATACGGTTCTTGAAGAGATCGACGCTCACGAGATCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCGCGTATTGATCGGGACGAAGAGAACAAACCGAACCGTGTCTGGCTTTCCGCACAGACCGGAGCGGGGATACCACAGCttt                                                                                                                       <  2:192382/251‑1 (MQ=255)
                                                                                                                          aCAAGAGACGCGGAAAGCCTCATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGGACAAGAAAACATCGAAGCGGTGAATACGGTTCTTGAAGAGATCGACGCTCACGAGATCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCGCGTATTGATCGGGACGAAGAGAACAAACCGAACCGTGTCTGGCTTTCCGCACAGACCGGAGCGGGGATACCACAGCTTTTTc                                                                                                                    <  2:59559/251‑1 (MQ=255)
                                                                                                                                              caTTACTGCTGCACGTCATTGGTGCGGCGGATGTGCGTGTACAAGAAAACAGCGAAGCGGTGAATTCTGTTCTTGAAGAGATCGACGCTCTCGAGATCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAATAGTTCGAACCGCGTATTGATCTGGATGAAGAGAACAAGCCGAACCGTGTCTGGCTTGTCGCACAGACCGGAGCGGGGATGCCACAGCTTTTGCAGGcatggactgtggggctg                                                                                                <  2:33515/251‑17 (MQ=255)
                                                                                                                                                  aCTGCTGCACGTCATTGATGCGGCGGATGTGCGTGTACAAGAAAACATCGAAGCGGTGAATACGGTTCTTGAAGAGATCGACGCTCACGAGATCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCGCGTATTGATCGGGACGAAGAGAACAAACCGAACCGTGTCTGGCTTTCCGCACAGACCGGAGCGGGGATACCACAGCTTTTTCAGGCTTTGACGGAGCGGCTTTcc                                                                                             <  1:142099/250‑1 (MQ=255)
                                                                                                                                                        gcACGTCATTGATGCGGCGGATGTGCGTGTACAAGAAAACATCGAAGCGGTGAATACGGTTCTTGAAGAGATCGACGCTCACGAGATCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCGCGTATTGATCGGGACGAAGAGAACAAACCGAACCGTGTCTGGCTTTCCGCACAGACCGGAGCGGGGATACCACAGCTTTTTCAGGCTTTGACGGAGCGGCTTTCCGGCGAgg                                                                                      <  1:347731/251‑1 (MQ=255)
                                                                                                                                                                  gATGCGGCGGATGTGCGTGTACAAGAAAAAATCGAAGCGGTGAATACGGTTCTTGAAGAGATCGACGCTCACGAGATCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCGCGTATTGGTCGGGACGAAGAGAACAAACCGAACCGTGTCTGGCTTTCCGCACAGACCGGAGCGGGGATACCACAGCTTTTTCCGGCTTTGACGGAGCGGCTTTCCGGCGAGGGGGCgcagc                                                                             >  1:48184/1‑250 (MQ=255)
                                                                                                                                                                       ggcggATGTGCGTGTACAAGAAAAAATCGAAGCGGTGAATACGGTTCTTGAAGAGATCGACGCTCACGAGATCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCGCGTATTGATCGGGACGAAGAGAACAAACCGAACCGTGTCTGGCTTTCCGCACAGACCGGAGCGGGGATACCACAGCATTTTCTGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATaaat                                                                       >  1:68677/1‑248 (MQ=255)
                                                                                                                                                                                    gtACAAGAAAACATCGAAGCGGTGAATACGGTTCTTGAAGAGATCGACGCTCACGAGATCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCGCGTATTGATCGGGACGAAGAGAACAAACCGAACCGTGTCTGGCTTTCCGCACAGACCGGAGCGGGGATACCACAGCTTTTTCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCAcc                                                          >  1:283131/1‑251 (MQ=255)
                                                                                                                                                                                       cAAGAAAACATCGAAGCGGTGAATACGGTTCTTGAAGAGATCGACGCTCACGAGATCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCGCGTATTGATCGGGACGAAGAGAACAAACCGAACCGTGTCTGGCTTTCCGCACAGACCGGAGCGGGGATACCACAGCTTTTTCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGc                                                        >  1:10940/1‑250 (MQ=255)
                                                                                                                                                                                         agcaaaCATCGAATCGATGACTACGGTTCTTGAAGAGATCGAAGAACACGAGATCACAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCGCGTATTGATCGGGACGAAGAGAACAAACCGAACAGTGTCTGGCTTTCCGCACAGACCGGAGCGGGGATACCACAGCTTTTTCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGccgg                                                     >  2:118239/4‑248 (MQ=255)
                                                                                                                                                                                          gAAAACATCGAAGCGGTGAATTCGGTTCTTGAAGAGAAGGCCGCATACGAGATCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCGCGTATTGATCGGGACGAAGAGAACAAACCGAACCGTGTCTGGCTTTCCGCACAGACCGGAGCGGGGATACCACAGCTTTTTCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGa                                                    <  1:150911/251‑1 (MQ=255)
                                                                                                                                                                                           aaaaCATCGAAGCGGTGAATACGGTTCTTGAAGAGATCGACGCTCACGAGATCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCGCGTATTGATCGGGACGAAGAGAACAAACCGAACCGTGTCTGGCTTTc                                                                                                                                                       <  1:120160/151‑1 (MQ=255)
                                                                                                                                                                                           aaaaCATCGAAGCGGTGAATACGGTTCTTGAAGAGATCGACGCTCACGAGATCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCGCGTATTGATCGGGACGAAGAGAACAAACCGAACCGTGTCTGGCTTTc                                                                                                                                                       >  2:120160/1‑151 (MQ=255)
                                                                                                                                                                                             aaCATCGAAGCTGTGAATACGGTTCTTGAAGAGATCGACGCTCACGAGATCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCGCGTATTGATCGGGACGAAGAGAACAAACCGATCCGTGTCTGGCTTTCCGCACAGACCGGAGCGGGGATACCACAGCTTGTTCAGGCTTTGACGGAGCGGCTGGCCGTCGAGGTGGCGCAGCATACATGGCGTCTGCCACCGCAGGAAgg                                                 <  2:264755/251‑1 (MQ=255)
                                                                                                                                                                                               cATCGAAGCGGTGAATACGGTTCTTGAAGAGATCGACGCTCACGAGATCCCAACTCTGCTGGTGATGAACCCGATCGATATGCTGGAAGATTTCGAACCGCGTATTGATCGGGACGACGAGAACAAACCGAACCGTGTCTGGCTTTCCGCACAGACCGGAGCGGGGATACCACAGCTTTTTCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGc                                               >  1:361876/1‑251 (MQ=255)
                                                                                                                                                                                                   gAAGCGGTGAATACGGTTCTTGAAGAGATCGACGCTCACGAGATCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCGCGTATTGATCGGGACGAAGAGAACAAACCGAACCGTGTCTGGCTTTCCGCACAGACCGGAGCGGGGATACCACAGCTTTTTCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCt                                           <  1:168881/251‑1 (MQ=255)
                                                                                                                                                                                                   gAAGCGGTGAATACGGTTCTTGAAGAGATCGACGCTCACGAGATCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCGCGTATTGATCGGGACGAAGAGAACAAACCGAACCGTGTCTGGCTTTCCGCACAGACCGGAGCGGGGATACCACAGCTTTTTCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCt                                           <  1:290576/251‑1 (MQ=255)
                                                                                                                                                                                                   gAAGCGGTGAATACGGTTCTTGAAGAGATCGACGCTCACGAGATCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCGCGTATTGATCGGGACGAAGAGAACAAACCGAACCGTGTCTGGCTTTCCGCACAGACCGGAGCGGGGATACCACAGCTTTTTCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCt                                           <  1:265515/251‑1 (MQ=255)
                                                                                                                                                                                                    tagagatGAATACGGTTCTTGAAGAGATCGACGCTCACGACATCCCAACGCTACTGGTGATGAACAAGAACGATATGCTGGAAGATTTCGTACAGCGTATTGCTCGGGACGAAGAGAACACACCGAACCGCGTCTGGCTTTCCGCACAGACCGGAGCGGGGATAAAACAGCTATTTCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGAGGCGCAGCGTACATTGCGTCTGCCCCCGCAGGCAGGGCGTCTg                                          >  2:18337/7‑251 (MQ=255)
                                                                                                                                                                                                    aaGCGGAGAATACGGTTCTTGAAGAGATCGACCCTCACGAGATCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCGCGTATTGATCGGGACGAAGAGAACAACCCGAACCGTGTCTGGCTTTCCGCACAGACCGGCGCGGGGATACCACAGCTTTTTCAGGCTTTGACGGAGCGGCTGTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTg                                          >  2:249686/1‑251 (MQ=255)
                                                                                                                                                                                                        ggTGAATACGGTTCTTGAAGAGATCGACGCTCACGAGATCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCGCGTATTGATCGGGACGAAGAGAACAAACCGAACCGTGTCTGGCTTTCCGCACAGACCGGAGCGGGGATACCACAGCTTTTTCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGaa                                      <  1:223072/251‑1 (MQ=255)
                                                                                                                                                                                                               aCGGTTCTTGAAGAGATCGACGCTCACGAGATCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCGCGTATTGATCGGGACGAAGAGAACAAACCGAACCGTGTCTGGCTTTCCGCACAGACCGGAGCGGGGATACCACAGCTTTTTCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGtt                                <  2:100667/250‑1 (MQ=255)
                                                                                                                                                                                                                 ggTTCTTGAAGAGATCGACGCTCACGAGATCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCGCGTATTGATCGGGACGAAGAGAACAAACCGAACCGTGTCTGGCTTTCCGCACAGACCGGAGCGGGGATACCACAGCTTTTTCAGGCTTTGCCGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCCCCGCCGGACGGGCGTCTGAGAAGCCGttttt                             >  2:387341/1‑251 (MQ=255)
                                                                                                                                                                                                                 ggTTCTTGAAGAGATCGACGCTCACGAGATCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCGCGTATTGATCGGGACGAAGAGAACAAACCGAACCGTGTCTGGCTTTCCGCACAGACCGGAGCGGGGATACCACAGCTTTTTCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGttttt                             >  1:292215/1‑251 (MQ=255)
                                                                                                                                                                                                                       tGAAGAGATCGACGCTCACGAGATCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCGCGTATTGATCGGGACGAAGAGAACAAACCGAACCGTGTCTGGCTTTCCGCACAGACCGGAGCGGGGATACCACAGCTTTTTCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATc                          <  1:152919/248‑1 (MQ=255)
                                                                                                                                                                                                                           gagaTCGACGCTCACGAGATCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCGCGTATTGATCGGGACGAAGAGAACAAACCGAACCGTGTCTGGCTTTCCGCACAGACCGGAGCGGGGATACCACAGCTTTTTCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCa                   <  1:110054/251‑1 (MQ=255)
                                                                                                                                                                                                                             gaTCGACGCTCGCGAGATCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCGCGTATTGATCGGGACGAAGAGAACAATCCGAACCGTGTCTGGCTTTCCGCACAGACCGGAGCGGGGATACCACAGCTGTTTCAGGCGTTGACGGAGCGGCTTTCCTGCGAGGTGGCGCAGCATACATGGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTGTTTATCAGCGTCAgg                 <  2:161410/251‑1 (MQ=255)
                                                                                                                                                                                                                                   cGCTCACGAGATCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCGCGTATTGATCGGGACGAAGAGAACAAACCGAACCGTGTCTGGCTTTCCGCACAGACCGGAGCGGGGATACCACAGCTTTTTCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAg           <  2:339284/251‑1 (MQ=255)
                                                                                                                                                                                                                                    gCTCACGCGATCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCGCGTATTGATCGGGACGAAGAGAACAAACCGAACCGTGTCTGGCTTTCCGCACAGACCGGAACGGGGATACCACAGCTTTTTCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCACCCGCAGGAAGGGCGTCTGAGACGTCGTTTTTATCAGCTTCAGGCAAt             >  2:37319/1‑248 (MQ=255)
                                                                                                                                                                                                                                       cACGAGATCCGAACTCTGCTGGTGATGAACCCGATCGATATGCTGGAAGATTTCGAACCGCGCATTGATCGGGACGAAGAGAACAAACCGCACCGTGTCTGGCTTTCCGCACAGACCGGCGCGGGGATACCACAGCTTTTTCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGCGGCGCAGCATACATTGCGTCTGCCACCGCCGGAAGGGCGTCTGGGAAGTCGTTTTTATCAGCTTCAGGCAATAGaaa        >  2:150600/1‑250 (MQ=255)
                                                                                                                                                                                                                                       cACGAGATCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCGCGTATTGATCGGGACGAAGAGAACAAACCGAACCGTGTCTGGCTTTCCGCACAGACCGGAGCGGGGATACCACAGCTTTTTCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGaaa        >  1:3533/1‑250 (MQ=255)
                                                                                                                                                                                                                                        aCGAGATCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCGCGTATTGATCGGGACGAAGAGAACAAACCGAACCGTGTCTGGCTTTCCGCACAGACCGGAGCGGGGATACCACAGCTTTTTCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGaaaaa      <  1:53647/251‑1 (MQ=255)
                                                                                                                                                                                                                                          gagaTCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCGCGTATTGATCGGTACGAAGAGAACAAACCGAACCGTGTCTGGCTTTCCGCACAGACCTGAGCGGGGATACCACAGCTTTTTCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTTTCAGCTTCAGGCAATAGAAAAaga    <  2:308506/251‑1 (MQ=255)
                                                                                                                                                                                                                                            gaTCCCAACTCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCGCGTATTGATCGGGACGAAGAGAACAAACCGAACCGTGTCTGGCTTTCCGCACAGACCGGAGCGGGGATACAACAGCTTTTTCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCCGGCAGGGCGTCTGAGAAGTCGTTTTTATCCGCTTCAGGCACTAGAAAAAGAGt  >  2:258323/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 
ACCCTCGACCCGACGTTGCGGCGTATTGACGTTGCAGATGTCGGTGAAACCGTACTTGCAGATACCGTAGGGTTTATTCGCCACCTGCCGCACGATCTGGTGGCGGCATTTAAAGCCACGTTACAAGAGACGCGGCAAGCCACATTACTGCTGCACGTCATTGATGCGGCGGATGTGCGTGTACAAGAAAACATCGAAGCGGTGAATACGGTTCTTGAAGAGATCGACGCTCACGAGATCCCAACCCTGCTGGTGATGAACAAGATCGATATGCTGGAAGATTTCGAACCGCGTATTGATCGGGACGAAGAGAACAAACCGAACCGTGTCTGGCTTTCCGCACAGACCGGAGCGGGGATACCACAGCTTTTTCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGT  >  NC_000913/4401362‑4401848

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: