Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,139,790 T→A 6.6% A471A (GCT→GCA 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,7900TA 6.6% 260.9 / 5.9 87A471A (GCT→GCAflgKflagellar hook‑filament junction protein 1
Reads supporting (aligned to +/- strand):  ref base T (42/39);  new base A (4/2);  total (46/41)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CATGGATGTATTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACGGAAATCTGCAAC  >  NC_000913/1139643‑1139938
                                                                                                                                                   |                                                                                                                                                    
cATGGATGTATTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTa                                                                                                                                                    >  2:1178347/1‑150 (MQ=255)
  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                                                                                                      >  2:379137/1‑150 (MQ=255)
                                              gaagaaGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGaaa                                                                                                      >  1:166815/1‑150 (MQ=255)
                                              gaagaaGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGaaa                                                                                                      >  1:559430/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                                                                                              >  2:1598431/1‑150 (MQ=255)
                                                      tGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAg                                                                                              >  1:1117757/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:410487/1‑150 (MQ=255)
                                                                       aCCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAAtg                                                                             <  2:1322637/150‑1 (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:812901/150‑1 (MQ=255)
                                                                          gCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATgtgg                                                                          <  2:85465/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                                                     <  2:944002/150‑1 (MQ=255)
                                                                                               aTCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATcagc                                                     <  1:1123062/150‑1 (MQ=255)
                                                                                               aTCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATcagc                                                     <  1:221973/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)
                                                                                                                                                   |                                                                                                                                                    
CATGGATGTATTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACGGAAATCTGCAAC  >  NC_000913/1139643‑1139938

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: