Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,134,118 A→C Q106P (CAG→CCG)  flgF → flagellar component of cell‑proximal portion of basal‑body rod

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,134,1180AC100.0% 248.3 / NA 69Q106P (CAG→CCG) flgFflagellar component of cell‑proximal portion of basal‑body rod
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (37/32);  total (37/32)

TGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCAGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAG  >  NC_000913/1133885‑1134346
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                     
tGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGGATACGCGGAATGGCAGCATTCCgg                                                                                                                                                                                                                                     >  1:199257/1‑236 (MQ=255)
                ttgcgcgagTAGGTGAATGCTTGACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCGCGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCa                                                                                                                                                                                                                     <  2:643967/233‑1 (MQ=255)
                   cgcgcgCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCCGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGGATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACt                                                                                                                                                                                                                  >  1:254047/1‑236 (MQ=255)
                   cgcgcgCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACt                                                                                                                                                                                                                  <  1:570051/236‑1 (MQ=255)
                     cgcgcAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGa                                                                                                                                                                                                                <  2:65548/236‑1 (MQ=255)
                     cgcgcAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTGCTCTGCCCACGCGCACGGTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGa                                                                                                                                                                                                                <  1:676967/236‑1 (MQ=255)
                                 tGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCGCGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGc                                                                                                                                                                                                    <  1:267046/236‑1 (MQ=255)
                                    tttACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCCCGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATc                                                                                                                                                                                                 <  1:645414/236‑1 (MQ=255)
                                    tttACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACGCCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATc                                                                                                                                                                                                 <  2:157489/236‑1 (MQ=255)
                                    tttACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATc                                                                                                                                                                                                 >  1:107159/1‑236 (MQ=255)
                                    tttACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATc                                                                                                                                                                                                 <  1:463927/236‑1 (MQ=255)
                                    tttACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATc                                                                                                                                                                                                 >  1:368891/1‑236 (MQ=255)
                                             ggTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGGATAAGCGGAAGTGCAGCCTTCCCGTTGATCCCACCGGGCAACAGACAATCCAGGGGCGACCggcggtgc                                                                                                                                                                                        >  2:393179/1‑230 (MQ=255)
                                               ttcaaGTGGAAGGGCGTGCCCGGCCCAAGCGCAAGTTGGTCACGGCGTCAACGCCGGGCGCAGAGATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAgg                                                                                                                                                                                       <  2:687846/231‑1 (MQ=255)
                                                  cAGTGGAAGGGCTTTCTCCGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAAGGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAgggg                                                                                                                                                                                                     <  1:335198/218‑1 (MQ=255)
                                                         aGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGc                                                                                                                                                                            >  1:107826/1‑236 (MQ=255)
                                                              tttCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGGTGATCCCACCGGGCAACTGACAATTCAGGGGCCTACGGTGAAAGGCGAGGCTGGGCCAAtt                                                                                                                                                                       >  2:67648/1‑236 (MQ=255)
                                                                 ctctGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGc                                                                                                                                                                     <  1:81096/235‑1 (MQ=255)
                                                                      cccACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCGGCATTCCGGTTGATCCCACCGGGCAACTGACAATTACGGGGCATACGGTGATAGGCGAGGCTGGGCGAATTGCTGTGcc                                                                                                                                                               >  2:30272/1‑236 (MQ=255)
                                                                        cACGCCCACGTTGGTCCGGCCGCACACGGCGGGCGCAGATATGACGCACGGCAAAGTGGATTACCCCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCgg                                                                                                                                                             <  2:634437/236‑1 (MQ=255)
                                                                          cgcgCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTTCCGGaa                                                                                                                                                           >  2:143489/1‑236 (MQ=255)
                                                                          cgcgCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCAGGTGATAGGCTAGGCTGGGCCAAATGCTGTGCCGGac                                                                                                                                                           >  2:276182/1‑235 (MQ=255)
                                                                           gcgcACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGGGCCGGAAg                                                                                                                                                          >  1:472456/1‑236 (MQ=255)
                                                                            cgcACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAgg                                                                                                                                                                                                       <  2:671582/190‑1 (MQ=255)
                                                                            cgcACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAgg                                                                                                                                                                                                       >  1:671582/1‑190 (MQ=255)
                                                                                tgTTGGTCACTGCGTCAACGCCGGGCGCAGATATGACGCCAGGCAAAATGGATTACGCCGCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCg                                                                                                                                                     <  2:107159/235‑1 (MQ=255)
                                                                                         cGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTg                                                                                                                                                                                            <  1:672244/188‑1 (MQ=255)
                                                                                         cGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTg                                                                                                                                                                                            >  2:672244/1‑188 (MQ=255)
                                                                                               cAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATtgct                                                                                                                                      >  1:724164/1‑236 (MQ=255)
                                                                                                 tttgCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATtgctg                                                                                                                                     >  1:463907/4‑236 (MQ=255)
                                                                                                    gCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCg                                                                                                                                  <  2:689126/236‑1 (MQ=255)
                                                                                                      cGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGAt                                                                                                                                <  2:472456/236‑1 (MQ=255)
                                                                                                       gggCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGc                                                                                                                                                                                      <  2:252959/181‑1 (MQ=255)
                                                                                                       gggCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGc                                                                                                                                                                                      >  1:252959/1‑181 (MQ=255)
                                                                                                        ggCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCCGGGGCATCCGGTGATAGGCGAGGCCGGGCCACTTGCTGTTCAGGTAGGGGCGGGAATCACTATTGCTGCCGAttg                                                                                                                              >  2:595216/1‑234 (MQ=255)
                                                                                                             aGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCGGGGCCAATTGCTGTGCCGGAATGGGGGGAAATCACCACTGATGCCGATGGcaca                                                                                                                          >  2:542129/1‑235 (MQ=255)
                                                                                                                   gACGCCCGGCAAAGTGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAAtctc                                                                                                                    <  1:321062/234‑1 (MQ=255)
                                                                                                                   gACGCCCGGCAAAAGGGATTACCCCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCg                                                                                                                   <  1:30272/235‑1 (MQ=255)
                                                                                                                    aCGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGgc                                                                                                                 >  1:113198/1‑236 (MQ=255)
                                                                                                                       cccgcaaaaaTGGATTCCCCCTCGCGCCCCGTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGgcg                                                                                                                <  2:199257/228‑1 (MQ=255)
                                                                                                                         cGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTc                                                                                                             >  1:606319/1‑235 (MQ=255)
                                                                                                                          ggCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCGGTGCCGGCAGGGGCGGAAATCACTATTGCTGCCGATGGCTCAAGTCTGGCGCTCaa                                                                                                           >  2:50787/1‑236 (MQ=255)
                                                                                                                           gCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAAt                                                                                                          <  1:669129/236‑1 (MQ=255)
                                                                                                                                      tACACCTCGCGCCGGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGAt                                                                                                 <  2:254047/234‑1 (MQ=255)
                                                                                                                                             cgcgCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGGAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCCCAATCCGGCCGATCCgggcaaa                                                                                        >  2:472158/1‑231 (MQ=255)
                                                                                                                                              gcgcCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGGGGTAGGCGAGGCTGGGCCAATTGCTGTGCCGGACGGGGCGGAAATCACTATTGCTGCCGCTGGCACAAATTCGGGCGTCAATCCCGGGGATCCTGCa                                                                                           >  2:6177/1‑232 (MQ=255)
                                                                                                                                               cgcCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCCGTGCCGGCAGGGGCGGAAATCACCATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCCGGCCATACGGCCAATAc                                                                                      >  2:686793/1‑236 (MQ=255)
                                                                                                                                                              gcgTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAgg                                                                       >  1:53487/1‑236 (MQ=255)
                                                                                                                                                              gcgTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAgg                                                                       <  2:113198/236‑1 (MQ=255)
                                                                                                                                                              gcgTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAgg                                                                       <  2:463907/236‑1 (MQ=255)
                                                                                                                                                                     agcagGAAGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAgg                                                                       <  1:311984/229‑1 (MQ=255)
                                                                                                                                                                       cagGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGgcg                                                                                                                <  1:448918/186‑1 (MQ=255)
                                                                                                                                                                            tggctggctggcCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGGGCACGTAGGGCGGCTTaaa                                                             >  2:311984/1‑232 (MQ=255)
                                                                                                                                                                              gctggctggcCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGGTTGCACCACCTCGGCGCTCAATc                                                                                                         >  2:448918/1‑186 (MQ=255)
                                                                                                                                                                                 ggctggcCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATcc                                                                                                                                                                                                <  2:333065/97‑1 (MQ=255)
                                                                                                                                                                                 ggctggcCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATcc                                                                                                                                                                                                >  1:333065/1‑97 (MQ=255)
                                                                                                                                                                                    tggcCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGcc                                                 >  1:713694/1‑236 (MQ=255)
                                                                                                                                                                                        cGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACgg                                             >  1:405790/1‑236 (MQ=255)
                                                                                                                                                                                         gTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGGCTTAAAATGGTGAAAAGCACgtg                                            >  2:586413/1‑234 (MQ=255)
                                                                                                                                                                                          tGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACggg                                            <  1:143489/235‑1 (MQ=255)
                                                                                                                                                                                          tGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGc                                           >  1:648070/1‑236 (MQ=255)
                                                                                                                                                                                            cAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAg                                         >  1:284045/1‑236 (MQ=255)
                                                                                                                                                                                                     gACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCa                                <  1:388952/236‑1 (MQ=255)
                                                                                                                                                                                                           aGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAgaccgg                          >  1:14903/1‑231 (MQ=255)
                                                                                                                                                                                                                  ggTATACGCGTAATGGAAGCATTCCGGTTGATCCCGCCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTgacga                    <  1:256543/235‑1 (MQ=255)
                                                                                                                                                                                                                            tAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTc         >  1:44694/1‑236 (MQ=255)
                                                                                                                                                                                                                             aaTGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACAGGCAGCGAAGGTCAGCCCGGGGGTGACGGCAttttcat        >  2:473835/1‑233 (MQ=255)
                                                                                                                                                                                                                              aTGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGt       >  2:339589/1‑236 (MQ=255)
                                                                                                                                                                                                                                   agcaTTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAg  <  2:648070/236‑1 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                     
TGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCAGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAG  >  NC_000913/1133885‑1134346

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