Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,403,010 C→T R325C (CGT→TGT)  hflK → regulator of FtsH protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,403,0100CT98.6% 220.6 / NA 69R325C (CGT→TGT) hflKregulator of FtsH protease
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (33/35);  minor base A (1/0);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGC  >  NC_000913/4402767‑4403255
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
ccgGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTCCTCGCGCGTGTCTGt                                                                                                                                                                                                                                                  >  2:375060/1‑250 (MQ=255)
   gaagaaGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGACGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTatat                                                                                                                                                                                                                                              >  1:20313/1‑251 (MQ=255)
   gaagaaGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTatat                                                                                                                                                                                                                                              >  2:156523/1‑251 (MQ=255)
   gaagaaGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTatat                                                                                                                                                                                                                                              <  1:254250/251‑1 (MQ=255)
   gaagaaGCAAAAGAGGCGTTTGACGAAGCGCTTGTCGCGCGTGAAAACGACCAGCAATACATTCGTGCAGCAGAAGCGTATAACAACGAAGTTCAGCCGCGTGCGAACGGTAAGGCGCAACGAATCATCGAAGAGGCGCGGGCGTACAAGGCCCAGACCATCCTGGAAGCTAAGGGTGAAGTGGCGCGCTTTGCTCAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAATGTCTGTatat                                                                                                                                                                                                                                              >  2:73144/1‑251 (MQ=255)
      gaagTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCga                                                                                                                                                                                                                                           <  1:190430/251‑1 (MQ=255)
           aaaaGCGGCGGTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACg                                                                                                                                                                                                                                       <  1:266283/250‑1 (MQ=255)
            aaaGCGGCGTTTGACGATGCGAGTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGGGGCGCGCTTTGCGAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACGCGCGAGTGTCTGTATATCGAGACGaa                                                                                                                                                                                                                                     <  2:313255/251‑2 (MQ=255)
                cggcgTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGa                                                                                                                                                                                                                                  >  2:111240/1‑250 (MQ=255)
                 ggcgTTTGACGATGCGATTGCCGCTCGTGAAAACGAACAGCAAGACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTGTCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGAGGCGTGCTTTGCTAAACGTCTGCCGTAATATAAAGCCGCGCCGGGAAGTGCTCGCTAGTGTCTGTATATCGAGACTATGGaaa                                                                                                                                                                                                                                <  2:305993/251‑1 (MQ=255)
                 ggcgTTTGACGATGCGATTGCCGCGCGTGAAAACGAACCGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCCGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGaa                                                                                                                                                                                                                                 >  2:222110/1‑250 (MQ=255)
                      ttGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTatat                                                                                                                                                                                                                                              <  2:53254/232‑1 (MQ=255)
                      ttGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTatat                                                                                                                                                                                                                                              >  1:53254/1‑232 (MQ=255)
                      ttGACGAGGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAgt                                                                                                                                                                                                                            <  2:398886/250‑1 (MQ=255)
                              gCGATTGCCGCGCATGAAAACGAACAGCAATATATTCGTGAAGCAGAAGCGGATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGATGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGTAATATAAGGCCGCGCCGGAGATTACTCTCGAGTGTCTGTGTATCGAGACGAGTGAAAAATTGTTGGGTaa                                                                                                                                                                                                                  <  2:73662/251‑1 (MQ=255)
                                gATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAac                                                                                                                                                                                                                 >  1:150713/1‑250 (MQ=255)
                                 aTTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAacac                                                                                                                                                                                                               >  1:119707/1‑251 (MQ=255)
                                    gcctcgcgTAAAAACTAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCTTGCGAACGTTCAGGCTCAACGTTTCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGTAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAACCCGCGCCGGAAATTACGCGCGAGTTTCTGTATCTCGCGACGATTCAAAAAGTGTTGGTTTACACCCg                                                                                                                                                                                                            <  2:411639/247‑1 (MQ=255)
                                     ccGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGc                                                                                                                                                                                                           >  1:253404/1‑251 (MQ=255)
                                          cgTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGTGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGAGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGt                                                                                                                                                                                                      <  2:315403/251‑1 (MQ=255)
                                                    aaCAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGCGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGCATATAAAGCCGAGCCGGAAATTACTCGCGAGTGTCTGTATCTCGAGACGATGGAAAAAGTGTTGGGTCACACCCGCAAAGTGCTGGTTAAc                                                                                                                                                                                            >  1:461215/1‑251 (MQ=255)
                                                       atcAATACCTTCGTTAAGCAGAAGCGTATTCCAACGAAGGTCTGCCGCGTGCGAACGGTCAGGCTTAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCAAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCGAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGAt                                                                                                                                                                                         <  2:350233/249‑1 (MQ=255)
                                                        gcaatgccttCGTGAAGCAGAAGCGTATACCACCGCAGTTCAGCCGCGTGCGAACGGTCAGTCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAACCTTCTGCCGGACTATCAAGCAGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGGAAAAGGGATGGGTAACCACCGCAAAGTGTTGGttcacact                                                                                                                                                                                         >  2:269553/9‑244 (MQ=255)
                                                         cAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGTGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATa                                                                                                                                                                                        >  2:164613/1‑250 (MQ=255)
                                                             aCATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAgg                                                                                                                                                                                    <  2:20313/250‑1 (MQ=255)
                                                              cATTCGTGAAGCAGAAGAGTATACAAACACAGGTCAGCCGCGTGCGACCGGTCAGGAGCAACGTATCCTCGAAGAGGCGCCTGCGTCCCAGGACAAACCAATCCTGGGAGCTCAGGGCGAAGTGGCGCGCTTTGCCAAACATCTGACGCAATATAACGCCGCGCCGGGAATTACTCGCGCGAGTCTGTAGATAGAGACGCTGGAAGAGGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAgatg                                                                                                                                                                                  >  2:38319/1‑248 (MQ=255)
                                                                     gacgCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGACAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATCTCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCt                                                                                                                                                                           >  2:37990/4‑251 (MQ=255)
                                                                         cAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGa                                                                                                                                                                         <  1:157378/249‑1 (MQ=255)
                                                                         cAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGAAACCTGa                                                                                                                                                                         >  2:157378/1‑249 (MQ=255)
                                                                            aaGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGtt                                                                                                                                                                    >  1:294391/1‑251 (MQ=255)
                                                                               cGTATACCAACGAAGTTCAGCCGCGTGCGAATGGTCAGGCGCAACGGATCCTCGAAGAGGCGTGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGAAGATGGAAAATGTGTTGGGTAACACCCTCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTg                                                                                                                                                                 <  2:401485/251‑1 (MQ=255)
                                                                               cGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAgtgt                                                                                                                                                                                                                          >  1:111001/1‑195 (MQ=255)
                                                                               cGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAgtgt                                                                                                                                                                                                                          <  2:111001/195‑1 (MQ=255)
                                                                                     cccaCGAAGTTCAGCCGCGTTCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAGCCTGATGGTTCTGCCGTTa                                                                                                                                                           <  1:3240/248‑1 (MQ=255)
                                                                                         agAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCGCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTACACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGAcc                                                                                                                                                       <  1:451135/250‑1 (MQ=255)
                                                                                          gAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGAcc                                                                                                                                                       >  1:37905/1‑250 (MQ=255)
                                                                                              ttCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCACAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATg                                                                                                                                                  >  1:408405/1‑251 (MQ=255)
                                                                                                cAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGc                                                                                                                                                 <  1:374016/250‑1 (MQ=255)
                                                                                                  gCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGa                                                                                                                                              <  2:325716/251‑1 (MQ=255)
                                                                                                           gAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTa                                                                                                                                     <  2:345181/251‑1 (MQ=255)
                                                                                                           gAACGGTCAGGCGCAACGTATCCTCGAAGAAGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGAGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAggtggt                                                                                                                                      >  2:59879/1‑250 (MQ=255)
                                                                                                               ggTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGGGTACAATGCCCAGACCATCCTGGAATCTCAGGGTGAAGTTGCGCGCTTTGCTAAACGTCTGCCGGAACATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGGTACGATGGAAAAAGGTTGGGCTAACACCCGCAAAGGGCTGGTTGACGATAAAGCTGGCAGCCTGATGGTTCTGCCGTTAGACCACATGCGGAAAGGTGGTAACGc                                                                                                                                 <  2:450101/251‑1 (MQ=255)
                                                                                                                  cAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGccc                                                                                                                               <  2:119707/250‑1 (MQ=255)
                                                                                                                          aCGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTgcggcg                                                                                                                       >  1:173722/1‑250 (MQ=255)
                                                                                                                                ccTCGAAGAGGCGCGTGCGTACAAGGCCCCGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCg                                                                                                                >  1:139576/1‑251 (MQ=255)
                                                                                                                                ccTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGc                                                                                                                 <  1:156523/250‑1 (MQ=255)
                                                                                                                                           cgcgTGCGTGCAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACGTCTGCCGGAATATAAAGCCGCGCCTGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTTCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGcc                                                                                                     <  1:260151/251‑1 (MQ=255)
                                                                                                                                           cgcgTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGcc                                                                                                     <  2:173722/251‑1 (MQ=255)
                                                                                                                                                 cgtACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAAt                                                                                               >  1:111502/1‑251 (MQ=255)
                                                                                                                                                    aCAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATctg                                                                                            >  1:240381/1‑251 (MQ=255)
                                                                                                                                                     aaagtcccAGACCATCCTGGAAGCTCAGGGTGAGGTGGCGTGCTTTGCTAAACTTCTGCCGGAATATAAATCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGAGGGAAAAAGTGTTGGTTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTTGTAACGCCCCTGCGGCGAAGAGCTATAACGGTGCCAGCAATctgc                                                                                           <  2:327711/246‑1 (MQ=255)
                                                                                                                                                                    tCCTGGAAGCTCAGGGTGAAGTAGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTgccgcc                                                                              >  1:342886/1‑250 (MQ=255)
                                                                                                                                                                           aGCTCAGGGTGAAGGGGCGCGCTTTGCTAAACATCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGGGTTGGGTAACACCCGCAACGTGCTGGTTAACGATAAAGGTGGCACCCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTTGTAACGCCCCTGCGGCGAAGAGCGATAACGGTTCCATCAATCTGCTGCGTCTACCCCCAGCctcgt                                                                      <  2:214872/251‑3 (MQ=255)
                                                                                                                                                                                             gcgcTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGag                                                    <  1:93341/251‑1 (MQ=255)
                                                                                                                                                                                              cgcTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGGTAGACCAGATGCTGAAAGGGGGTAACGCCCCTGCGGCGAAGTGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGagc                                                   <  2:176876/251‑1 (MQ=255)
                                                                                                                                                                                                 tttGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGagca                                                  >  2:330186/1‑249 (MQ=255)
                                                                                                                                                                                                            tCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACAcgtcg                                      <  1:330186/250‑1 (MQ=255)
                                                                                                                                                                                                            tCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACAcgtc                                       <  1:246505/249‑1 (MQ=255)
                                                                                                                                                                                                             cTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACgtcgtc                                    >  2:316711/1‑251 (MQ=255)
                                                                                                                                                                                                                         aaaGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGCGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGCGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAggg                        >  2:232658/1‑251 (MQ=255)
                                                                                                                                                                                                                         aaaGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAACAGTGTTGGGTAACACCCGCAAAGAGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAggg                        >  2:224162/1‑251 (MQ=255)
                                                                                                                                                                                                                         aaaGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAggg                        >  1:78532/1‑251 (MQ=255)
                                                                                                                                                                                                                            gCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCgg                     >  1:88348/1‑250 (MQ=255)
                                                                                                                                                                                                                              cgcgCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGCAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGACGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGat                    >  1:58593/1‑250 (MQ=255)
                                                                                                                                                                                                                                  ccGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGAtattat               <  2:37905/251‑1 (MQ=255)
                                                                                                                                                                                                                                        aTTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCa         <  2:240381/251‑1 (MQ=255)
                                                                                                                                                                                                                                        aTTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCa         <  1:222110/251‑1 (MQ=255)
                                                                                                                                                                                                                                              cgcgAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTAGGGACCAAcgc     <  2:294391/249‑1 (MQ=255)
                                                                                                                                                                                                                                               tcgAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAAcgccgc  >  2:65791/2‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
CCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGC  >  NC_000913/4402767‑4403255

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