Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,139,793 T→A 6.9% Y472* (TAT→TAA flgK → flagellar hook‑filament junction protein 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,139,7930TA 6.9% 265.6 / 8.0 87Y472* (TAT→TAAflgKflagellar hook‑filament junction protein 1
Reads supporting (aligned to +/- strand):  ref base T (41/39);  new base A (4/2);  total (46/41)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.10e-01

TGGATGTATTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACGGAAATCTGCAACGTTT  >  NC_000913/1139645‑1139942
                                                                                                                                                    |                                                                                                                                                     
tGGATGTATTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGTCCGTGCTGTATCTGCAAAGCAAGAGTAATACTGTGTTCGGTGCGAATTCCTTTAACGACGCTTATg                                                                                                                                                      <  2:895225/150‑1 (MQ=255)
tGGATGTATTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATg                                                                                                                                                      >  2:872447/1‑150 (MQ=255)
    tGTATTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTc                                                                                                                                                  <  1:300575/150‑1 (MQ=255)
      tATTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGt                                                                                                                                                >  2:1054156/1‑150 (MQ=255)
         tAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTag                                                                                                                                             <  1:1512119/150‑1 (MQ=255)
          aaTCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTagt                                                                                                                                            >  2:209490/1‑150 (MQ=255)
          aaTCACCGACGAAGCGAAAATAGCGAGGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACAGCAACGGTCAGGCCCTGCTGGATCTGGAAAGCAACAGTAGCCCGGTGGGCGGTGCGAAATCCTTTAACGGCGGTTAGGCCTCGTTagt                                                                                                                                            >  2:301203/1‑150 (MQ=255)
           aTCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTagtg                                                                                                                                           <  1:961118/150‑1 (MQ=255)
            tCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGAGGCATAAGCGTCGTTaaagg                                                                                                                                          <  2:904698/150‑5 (MQ=255)
             cACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTagtgag                                                                                                                                         <  2:1415437/150‑1 (MQ=255)
              aCCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTagtgagt                                                                                                                                        >  1:1429607/1‑150 (MQ=255)
               ccGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTagtgagtg                                                                                                                                       >  1:379420/1‑150 (MQ=255)
                  acgaAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGata                                                                                                                                    >  2:1186983/1‑150 (MQ=255)
                       gCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGt                                                                                                                               <  2:1077287/150‑1 (MQ=255)
                       gCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGt                                                                                                                               >  1:1538751/1‑150 (MQ=255)
                       gCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGt                                                                                                                               >  1:477131/1‑150 (MQ=255)
                         gAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGtaa                                                                                                                             >  2:221973/1‑150 (MQ=255)
                         gAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGtaa                                                                                                                             >  2:746344/1‑150 (MQ=255)
                         gAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGtaa                                                                                                                             >  1:515397/1‑150 (MQ=255)
                            aaTAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGtaataa                                                                                                                          >  1:944002/1‑150 (MQ=255)
                              tAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCATGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATaaaa                                                                                                                        >  2:1088913/1‑150 (MQ=255)
                                 cGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACcg                                                                                                                     >  2:1017560/1‑150 (MQ=255)
                                   aTGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACcgcg                                                                                                                   <  2:528749/150‑1 (MQ=255)
                                   aTGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACcgcg                                                                                                                   <  2:968377/150‑1 (MQ=255)
                                    tGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGa                                                                                                                  <  2:885428/150‑1 (MQ=255)
                                      gcgagcgaAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGGTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACg                                                                                                                <  2:1470633/150‑1 (MQ=255)
                                      gcgagcgaAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACg                                                                                                                <  1:1449491/150‑1 (MQ=255)
                                      gcgagcgaAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACg                                                                                                                <  2:1513558/150‑1 (MQ=255)
                                           cgaAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGaa                                                                                                           >  1:348098/1‑150 (MQ=255)
                                           cgaAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGaa                                                                                                           >  1:1596490/1‑150 (MQ=255)
                                            gaagaaGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGaaa                                                                                                          >  1:166815/1‑150 (MQ=255)
                                            gaagaaGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGaaa                                                                                                          >  1:559430/1‑150 (MQ=255)
                                            gaagaaGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGaaa                                                                                                          >  2:379137/1‑150 (MQ=255)
                                             aagaagATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGaaaa                                                                                                         <  1:1017560/150‑1 (MQ=255)
                                             aagaagATGCAGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGaaaa                                                                                                         >  2:151962/1‑150 (MQ=255)
                                               gaagaTGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAAcc                                                                                                       <  2:914635/150‑1 (MQ=255)
                                                 agaTGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAg                                                                                                     >  1:1044845/1‑150 (MQ=255)
                                                    tGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAg                                                                                                  >  1:1117757/1‑150 (MQ=255)
                                                    tGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAg                                                                                                  >  2:1598431/1‑150 (MQ=255)
                                                      cGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAgcg                                                                                                >  2:1018920/1‑150 (MQ=255)
                                                       gggTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGTGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAgcgc                                                                                               <  2:936841/150‑1 (MQ=255)
                                                                cGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAgg                                                                                      >  2:77182/1‑150 (MQ=255)
                                                                cGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAgg                                                                                      <  2:1044845/150‑1 (MQ=255)
                                                                 gATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGt                                                                                     <  1:1186983/150‑1 (MQ=255)
                                                                     aCCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAAtg                                                                                 >  1:1351782/1‑150 (MQ=255)
                                                                     aCCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAAtg                                                                                 >  1:133292/1‑150 (MQ=255)
                                                                     aCCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAAtg                                                                                 <  2:1322637/150‑1 (MQ=255)
                                                                     aCCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAAtg                                                                                 >  2:410487/1‑150 (MQ=255)
                                                                      ccGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTGATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAAtgt                                                                                >  1:300084/1‑150 (MQ=255)
                                                                      ccGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAAtgt                                                                                >  1:706583/1‑150 (MQ=255)
                                                                        gCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAAtg                                                                                 <  1:130299/147‑1 (MQ=255)
                                                                        gCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAAtg                                                                                 >  2:130299/1‑147 (MQ=255)
                                                                        gCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATgtgg                                                                              <  2:85465/150‑1 (MQ=255)
                                                                        gCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATgtgg                                                                              <  2:812901/150‑1 (MQ=255)
                                                                        gCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATgtgg                                                                              <  1:1178347/150‑1 (MQ=255)
                                                                              gTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGc                                                                        <  1:1244633/150‑1 (MQ=255)
                                                                                  ggCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCt                                                                    <  2:1098439/150‑1 (MQ=255)
                                                                                     ccTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTc                                                                 <  1:1286567/150‑1 (MQ=255)
                                                                                       tgctgGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCa                                                               <  2:311876/150‑1 (MQ=255)
                                                                                           ggATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATca                                                           >  1:674808/1‑150 (MQ=255)
                                                                                             aTCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATcagc                                                         <  1:221973/150‑1 (MQ=255)
                                                                                             aTCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATcagc                                                         <  2:944002/150‑1 (MQ=255)
                                                                                             aTCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATcagc                                                         <  1:1123062/150‑1 (MQ=255)
                                                                                                 gCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATcagcagca                                                     >  1:1537183/1‑150 (MQ=255)
                                                                                                  cAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATcagcagcag                                                    >  1:1463228/1‑150 (MQ=255)
                                                                                                   aaaGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGt                                                   >  1:442432/1‑150 (MQ=255)
                                                                                                    aaGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTc                                                  >  1:692208/1‑150 (MQ=255)
                                                                                                            gTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCg                                          <  1:1088913/150‑1 (MQ=255)
                                                                                                            gTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCg                                          <  2:477131/150‑1 (MQ=255)
                                                                                                             tAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCgg                                         <  1:209490/150‑1 (MQ=255)
                                                                                                               aaaCGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGgtg                                       <  2:559430/150‑1 (MQ=255)
                                                                                                               aaaCGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGgtg                                       <  2:1538751/150‑1 (MQ=255)
                                                                                                                 aCGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTc                                     <  2:515397/150‑1 (MQ=255)
                                                                                                                           gTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATg                           >  1:1226790/1‑150 (MQ=255)
                                                                                                                             gCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATgag                         <  1:872447/150‑1 (MQ=255)
                                                                                                                             gCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATgag                         <  2:157377/150‑1 (MQ=255)
                                                                                                                                    ccTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACg                  <  2:1429607/150‑1 (MQ=255)
                                                                                                                                     cactAACGAGGCATAAGCGTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACgg                 >  1:755223/4‑150 (MQ=255)
                                                                                                                                      tttAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACGGa                >  1:683997/1‑150 (MQ=255)
                                                                                                                                       ctAACGAGGCATAAGCGTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATc                                                            >  2:1233127/2‑105 (MQ=255)
                                                                                                                                       ctAACGAGGCATAAGCGTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATc                                                            <  1:1233127/104‑1 (MQ=255)
                                                                                                                                       ctAACGAGGCATAAGCGTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACGGaa               >  1:86534/2‑150 (MQ=255)
                                                                                                                                        tAACGAGGCATAAGCGTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCAAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACGGaaa              >  1:446100/1‑150 (MQ=255)
                                                                                                                                        tAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACGGaaa              >  2:833480/1‑150 (MQ=255)
                                                                                                                                             acgCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACGGAAATCTGc         <  1:151962/150‑1 (MQ=255)
                                                                                                                                                cTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACGGAAATCTGCAAc      >  2:822181/1‑150 (MQ=255)
                                                                                                                                                  tATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACGGAAATCTGCAACGt    >  1:3801/1‑150 (MQ=255)
                                                                                                                                                    tGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACGGAAATCTGCAACGttt  <  2:176859/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGGATGTATTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACGGAAATCTGCAACGTTT  >  NC_000913/1139645‑1139942

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