Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,856,836 T→A 6.5% intergenic (‑30/‑257) ygbA ← / → mutS protein YgbA/DNA mismatch repair protein MutS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,856,8360TA 6.5% 250.5 / 5.8 86intergenic (‑30/‑257)ygbA/mutSprotein YgbA/DNA mismatch repair protein MutS
Reads supporting (aligned to +/- strand):  ref base T (34/46);  new base A (4/2);  total (38/48)
Fisher's exact test for biased strand distribution p-value = 3.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

AACAAACAACGCCTCGTAATGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTC  >  NC_000913/2856690‑2856982
                                                                                                                                                  |                                                                                                                                                  
aacaaacaACGCCTCGTAATGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATAtt                                                                                                                                                 >  1:207773/1‑150 (MQ=255)
     acaacGCCTCGTAATGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACatc                                                                                                                                            <  1:1246321/150‑1 (MQ=255)
           ccTCGTAATGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACatcttaatc                                                                                                                                      >  1:1301192/1‑150 (MQ=255)
           ccTCGTAATGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACatcttaatc                                                                                                                                      <  2:585299/150‑1 (MQ=255)
           ccTCGTAATGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACatcttaatc                                                                                                                                      >  1:569053/1‑150 (MQ=255)
            cTCGTAATGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACatcttaatct                                                                                                                                     <  1:987459/150‑1 (MQ=255)
             tCGTAATGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACatcttaatctt                                                                                                                                    <  1:194577/150‑1 (MQ=255)
             tCGTAATGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTAATATAAATACatcttaatgag                                                                                                                                    <  1:495452/150‑4 (MQ=255)
              cGTAATGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGTTGCCTCCACCTCATTAAGATGTATTTATATTACatcttaatctta                                                                                                                                   <  2:589803/150‑1 (MQ=255)
              cGTAATGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACatcttaatctta                                                                                                                                   <  2:1166478/150‑1 (MQ=255)
                tAATGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTaaa                                                                                                                                 >  2:357985/1‑150 (MQ=255)
                 aaTGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAg                                                                                                                                <  1:294920/150‑1 (MQ=255)
                   tGCTCCGGCTCCGCTGACGCCTGCGGTCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAggg                                                                                                                              <  1:600066/150‑1 (MQ=255)
                   tGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAggg                                                                                                                              <  1:539921/150‑1 (MQ=255)
                    gctccggctccgCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGc                                                                                                                             <  2:410663/150‑1 (MQ=255)
                    gctccggctccgCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTCTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGc                                                                                                                             >  1:349131/1‑150 (MQ=255)
                          gctccgCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATg                                                                                                                       <  2:346454/150‑1 (MQ=255)
                          gctccgCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATg                                                                                                                       <  2:487775/150‑1 (MQ=255)
                             ccgCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACt                                                                                                                    >  2:1331413/1‑150 (MQ=255)
                               gCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTcc                                                                                                                  <  1:1229213/150‑1 (MQ=255)
                                  gACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCaaa                                                                                                               <  1:22319/150‑1 (MQ=255)
                                   aCGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAaag                                                                                                              <  2:446040/150‑1 (MQ=255)
                                    cGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAaaga                                                                                                             >  1:541153/1‑150 (MQ=255)
                                    cGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAaaga                                                                                                             <  1:74393/150‑1 (MQ=255)
                                    cGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAaaga                                                                                                             >  1:1328739/1‑150 (MQ=255)
                                     gCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAaagaa                                                                                                            >  2:569613/1‑150 (MQ=255)
                                       cTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTCCATCTTAATCTTAAAGGGCACTATGACTCCAaagaaga                                                                                                          <  1:416558/150‑1 (MQ=255)
                                           gggCATTTCTCTTGATAAAGATCGATCATTTTTTTATTCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAATATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAggg                                                                                                      <  1:994917/150‑1 (MQ=255)
                                            ggCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGt                                                                                                     >  1:1229198/1‑150 (MQ=255)
                                             gCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGtt                                                                                                    <  2:861941/150‑1 (MQ=255)
                                             gCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGATGCCTCCACATCATTAAGATGTATTTATATTACATATTAATATTAAAGGGCACTATGAATCAAAAGAAGAAGGGtt                                                                                                    >  1:159638/1‑150 (MQ=255)
                                              cATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTa                                                                                                   >  1:1042085/1‑150 (MQ=255)
                                              cATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTa                                                                                                   >  1:1037136/1‑150 (MQ=255)
                                                tttCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGc                                                                                                 <  2:122729/150‑1 (MQ=255)
                                                    gCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAAc                                                                                             <  2:1301192/150‑1 (MQ=255)
                                                     cTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAAcc                                                                                            <  1:1240020/150‑1 (MQ=255)
                                                       tGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGa                                                                                          >  2:922906/1‑150 (MQ=255)
                                                        gATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGAt                                                                                         >  2:713839/1‑150 (MQ=255)
                                                           aaaGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACa                                                                                      >  2:1112993/1‑150 (MQ=255)
                                                             aGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAAt                                                                                    <  1:225978/150‑1 (MQ=255)
                                                             aGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAAt                                                                                    <  2:849860/150‑1 (MQ=255)
                                                                        ttttttAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACtt                                                                         <  1:1029242/150‑1 (MQ=255)
                                                                          ttttAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGc                                                                       >  2:1423917/1‑150 (MQ=255)
                                                                            ttAATCGTCAGTTTTTCAGGTGAGTTTCGCTTGCGGGACATTCTGCCTGAACATGATTTAGATGTTTTTATATTACATCTTAATCTTAAAGGGCACTATGTCTCCATAGATGAAGGGTTAGCCAAGCGTTACAATTTTGCGTACTTGCtt                                                                     <  2:1140675/150‑1 (MQ=255)
                                                                            ttAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCtt                                                                     <  2:557464/150‑1 (MQ=255)
                                                                                 cGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATaa                                                                <  2:207773/150‑1 (MQ=255)
                                                                                 cGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATaa                                                                >  2:1293814/1‑150 (MQ=255)
                                                                                  gtaaGTTTTTCACGAGAGATAAGCTTGCCGGACATGCTGCCTCCAACTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTAATTGCTTAATAAg                                                               >  2:1253191/4‑150 (MQ=255)
                                                                                  gTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATaaa                                                               >  2:1315934/1‑149 (MQ=255)
                                                                                  gTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAg                                                               >  1:957526/1‑150 (MQ=255)
                                                                                  gTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAg                                                               >  1:1334972/1‑150 (MQ=255)
                                                                                  gTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAg                                                               <  1:1310294/150‑1 (MQ=255)
                                                                                    cAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCa                                                             <  1:1331413/150‑1 (MQ=255)
                                                                                       tttttCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCa                                                          <  1:1048617/150‑1 (MQ=255)
                                                                                            cACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCaa                                                     >  1:1103959/1‑150 (MQ=255)
                                                                                             aCGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCaaa                                                    >  1:1119036/1‑150 (MQ=255)
                                                                                              cGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCaaaa                                                   <  1:1422891/150‑1 (MQ=255)
                                                                                                         cTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACa                                        <  1:1293814/150‑1 (MQ=255)
                                                                                                            gCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCa                                     <  1:1017949/150‑1 (MQ=255)
                                                                                                                    tgctgcCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGCCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTACATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTccc                             >  1:1180703/1‑150 (MQ=255)
                                                                                                                     gctgcCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCg                            <  2:159638/150‑1 (MQ=255)
                                                                                                                          cTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAAcc                       <  1:713839/150‑1 (MQ=255)
                                                                                                                           tCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCa                      >  2:1318373/1‑150 (MQ=255)
                                                                                                                           accACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGACACTATGACTCCAAAGAAGAAGGGTTAGCCAACAGATACAATTTTGCGTAATTGCTTCATAAGCATAACGCAAAAGCTGCAAAACAGCATCTTTCCCGTCACCa                      >  1:903805/2‑150 (MQ=255)
                                                                                                                            ccACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTc                               >  2:797108/1‑140 (MQ=255)
                                                                                                                            ccACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTc                               <  1:797108/140‑1 (MQ=255)
                                                                                                                            ccACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAg                     >  1:205466/1‑150 (MQ=255)
                                                                                                                              aCCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCa                   <  2:1103959/150‑1 (MQ=255)
                                                                                                                               ccTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCAt                  >  2:1321210/1‑150 (MQ=255)
                                                                                                                                cTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATc                 <  1:729643/150‑1 (MQ=255)
                                                                                                                                aagaTTAAGATGTAATATAAATACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATc                 >  1:1027791/4‑150 (MQ=255)
                                                                                                                                 tCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCa                <  2:1201499/150‑1 (MQ=255)
                                                                                                                                 tCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCa                >  2:488469/1‑150 (MQ=255)
                                                                                                                                 agaTTAAGATGTAATATAAATACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCa                >  1:515072/3‑150 (MQ=255)
                                                                                                                                 agaTTAAGATGTAATATAAATACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCa                >  2:130528/3‑150 (MQ=255)
                                                                                                                                  cATTAAGATGTAATATAAATACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCaa               >  2:374641/1‑150 (MQ=255)
                                                                                                                                   aTTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAg              <  2:1229198/150‑1 (MQ=255)
                                                                                                                                   aTTAAGATGTAATATAAATACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAg              <  2:515072/150‑1 (MQ=255)
                                                                                                                                    ttAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGCCTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGa             >  2:1194283/1‑150 (MQ=255)
                                                                                                                                    ttAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGa             <  2:1334972/150‑1 (MQ=255)
                                                                                                                                    ttAAGATCTATTTATATTACATCTAAATCTTAAAGGGCACAATGACTCCAAAGAACAAGGGTTAGCGAACCGATACAATTTTGCGTACTAGCTTCATAAGCATCACGCAAAAGCTGCAAAACCGCATCTTTCCCGGAACCAGCATCAAGa             >  2:871682/1‑150 (MQ=255)
                                                                                                                                       aGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACt          <  1:1423917/150‑1 (MQ=255)
                                                                                                                                       aGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACt          <  1:1321210/150‑1 (MQ=255)
                                                                                                                                        taTGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTc         <  1:357985/149‑1 (MQ=255)
                                                                                                                                            tATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCg     >  1:220514/1‑150 (MQ=255)
                                                                                                                                               ttATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTc  <  2:349131/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AACAAACAACGCCTCGTAATGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTC  >  NC_000913/2856690‑2856982

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: