Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,856,832 T→A 6.9% intergenic (‑26/‑261) 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,8320TA 6.9% 264.0 / 6.4 85intergenic (‑26/‑261)ygbA/mutSprotein YgbA/DNA mismatch repair protein MutS
Reads supporting (aligned to +/- strand):  ref base T (34/45);  new base A (4/2);  total (38/47)
Fisher's exact test for biased strand distribution p-value = 4.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

GCGCGTAAACAAACAACGCCTCGTAATGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCG  >  NC_000913/2856683‑2856979
                                                                                                                                                     |                                                                                                                                                   
gcgcGTAAACAAACAACGCCTCGTAATGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTAt                                                                                                                                                     >  1:487775/1‑150 (MQ=255)
 cgcgTAAACAAACAACGCCTCGTAATGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTAtt                                                                                                                                                    <  2:1031131/150‑1 (MQ=255)
       aacaaacaACGCCTCGTAATGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATAtt                                                                                                                                              >  1:207773/1‑150 (MQ=255)
            acaacGCCTCGTAATGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACatc                                                                                                                                         <  1:1246321/150‑1 (MQ=255)
                  ccTCGTAATGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACatcttaatc                                                                                                                                   >  1:569053/1‑150 (MQ=255)
                  ccTCGTAATGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACatcttaatc                                                                                                                                   <  2:585299/150‑1 (MQ=255)
                  ccTCGTAATGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACatcttaatc                                                                                                                                   >  1:1301192/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:74393/150‑1 (MQ=255)
                                           cGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAaaga                                                                                                          >  1:541153/1‑150 (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:1037136/1‑150 (MQ=255)
                                                     cATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTa                                                                                                >  1:1042085/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:1310294/150‑1 (MQ=255)
                                                                                         gTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAg                                                            >  1:1334972/1‑150 (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                            <  1:797108/140‑1 (MQ=255)
                                                                                                                                   ccACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTc                            >  2:797108/1‑140 (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             >  2:130528/3‑150 (MQ=255)
                                                                                                                                        agaTTAAGATGTAATATAAATACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCa             >  1:515072/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)
                                                                                                                                                     |                                                                                                                                                   
GCGCGTAAACAAACAACGCCTCGTAATGCTCCGGCTCCGCTGACGCCTGCGGGCATTTCGCTTGATAAAGATCGATCATTTTTTTAATCGTCAGTTTTTCACGAGAGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCG  >  NC_000913/2856683‑2856979

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: