Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,856,887 C→T 81.4% intergenic (‑81/‑206) ygbA ← / → mutS uncharacterized protein/methyl‑directed mismatch repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,856,8870CT81.4% 209.8 / 39.8 100intergenic (‑81/‑206)ygbA/mutSuncharacterized protein/methyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base C (11/8);  new base T (40/41);  total (51/49)
Fisher's exact test for biased strand distribution p-value = 6.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTCGCT  >  NC_000913/2856791‑2856985
                                                                                                 |                                                                                                  
tACGCTTGCCGTACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGTTTAGTCa                                                                                                  <  1:665392/99‑1 (MQ=255)
 aCGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAAc                                                                                                >  2:1701172/1‑100 (MQ=255)
  cGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAAcc                                                                                               >  2:2813326/1‑100 (MQ=255)
   gCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTCTTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGTTTATTCAACCg                                                                                              >  1:561801/1‑100 (MQ=255)
   gCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAAcc                                                                                               >  1:2448680/1‑99 (MQ=255)
   gCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCg                                                                                              <  2:2287773/100‑1 (MQ=255)
     ttGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGAt                                                                                            <  2:191793/100‑1 (MQ=255)
      tGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATAGTACATCTTAATCGTAAAGGGCACTATGAGTCCAAAGAAGAAGGGTGAGCCAACCGATa                                                                                           >  1:168780/1‑100 (MQ=255)
           gACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAAtt                                                                                      >  2:1505365/1‑100 (MQ=255)
            aCATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAAttt                                                                                     <  2:247234/100‑1 (MQ=255)
              aTGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTg                                                                                   <  2:1131581/100‑1 (MQ=255)
               tgctgcCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGc                                                                                  >  2:2675270/1‑100 (MQ=255)
               tgctgcCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGc                                                                                  <  1:1001514/100‑1 (MQ=255)
                gctgcCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCg                                                                                 >  1:432279/1‑100 (MQ=255)
                gctgcCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCg                                                                                 <  1:1541262/100‑1 (MQ=255)
                     tctcACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACtt                                                                            <  2:218694/97‑1 (MQ=255)
                      tCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTg                                                                           <  2:297385/100‑1 (MQ=255)
                       ccACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGc                                                                          >  2:138041/1‑100 (MQ=255)
                        cACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACttaga                                                                         >  1:218693/1‑97 (MQ=255)
                        cACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCt                                                                         <  1:1083328/100‑1 (MQ=255)
                          ccTCATTAAGATGTATTTATATTACATCTTAATCTTATAGGGCACTGTGACTCAAACGAAGAAGGACTAGTAAACCGATACAATTTTGCGTACTTGATTc                                                                       >  1:2209358/1‑100 (MQ=255)
                          ccTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTc                                                                       <  1:1360730/100‑1 (MQ=255)
                           cTCCTTAAGAGGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTc                                                                       <  2:1599134/99‑1 (MQ=255)
                           cTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCa                                                                      <  1:1631526/100‑1 (MQ=255)
                             gaTTAAGATGTAATATAAATACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATa                                                                    >  1:2443524/2‑100 (MQ=37)
                               ttAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGTTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAg                                                                  >  1:2259158/1‑100 (MQ=255)
                               ttAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAg                                                                  >  1:118739/1‑100 (MQ=255)
                                 aaGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCa                                                                >  1:1475030/1‑100 (MQ=255)
                                 aaGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCa                                                                >  2:2604185/1‑100 (MQ=255)
                                 aaGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATGTTGCGTACTTGCTTCATAAGCa                                                                <  1:1693630/100‑1 (MQ=255)
                                 aaGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCa                                                                <  2:919706/100‑1 (MQ=255)
                                   gATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATc                                                              >  1:589503/1‑100 (MQ=255)
                                    aTGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCa                                                             <  2:1622579/100‑1 (MQ=255)
                                      gTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACg                                                           <  2:2697394/100‑1 (MQ=255)
                                            tataTTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCaaaa                                                      <  2:2774882/100‑1 (MQ=255)
                                            tataTTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCaaaa                                                      >  1:974717/1‑100 (MQ=255)
                                            tataTTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCaaaa                                                      <  1:2580174/100‑1 (MQ=255)
                                                 tACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGc                                                 >  2:966760/1‑100 (MQ=255)
                                                  aCATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCa                                                <  1:1785695/100‑1 (MQ=255)
                                                     tcttaatcttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCaaaa                                             >  1:1087688/1‑100 (MQ=255)
                                                      cttaatcttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAAc                                            <  1:447983/100‑1 (MQ=255)
                                                      cttaatcttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAAc                                            >  2:1860824/1‑100 (MQ=255)
                                                       ttaatcttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACa                                           >  2:2140827/1‑100 (MQ=255)
                                                         aatcttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGc                                         <  1:2813308/100‑1 (MQ=255)
                                                         aatcttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGc                                         <  2:1729035/100‑1 (MQ=255)
                                                         aatcttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAGAACAGc                                         >  2:387770/1‑100 (MQ=255)
                                                         aatcttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGc                                         >  1:535911/1‑100 (MQ=255)
                                                          atcttaaAGGGCACTATGCCTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCa                                        <  1:2671760/100‑1 (MQ=255)
                                                          atcttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCa                                        >  1:431310/1‑100 (MQ=255)
                                                          atcttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCa                                        <  1:2818538/100‑1 (MQ=255)
                                                          atcttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCa                                        >  2:2439763/1‑100 (MQ=255)
                                                          atcttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCa                                        >  2:964813/1‑100 (MQ=255)
                                                            cttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCCTCACGCAAAAGCTGCAAAACAGCATc                                      >  1:2180265/1‑100 (MQ=255)
                                                            cttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCa                                        >  1:396176/1‑98 (MQ=255)
                                                            cttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATc                                      >  1:664259/1‑100 (MQ=255)
                                                            cttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTCCGTACTTGCTTCATACACATCACGCAAAAGCTGCAAAACAGCATc                                      >  1:982144/1‑100 (MQ=255)
                                                            cttaaAGGGCAATATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCa                                        <  1:1982594/98‑1 (MQ=255)
                                                              taaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCtt                                    >  2:1760903/1‑100 (MQ=255)
                                                                 aGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTcc                                 >  1:978050/1‑100 (MQ=255)
                                                                  gggCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTccc                                <  2:2810739/100‑1 (MQ=255)
                                                                  gggCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTccc                                >  2:495272/1‑100 (MQ=255)
                                                                  gggCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTccc                                <  1:2242154/100‑1 (MQ=255)
                                                                  gggCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTccc                                >  2:271722/1‑100 (MQ=255)
                                                                   ggCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCg                               <  2:1241842/100‑1 (MQ=255)
                                                                   ggCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCg                               <  2:2551475/100‑1 (MQ=255)
                                                                    gCACTATGACTCCAAAGAAGCAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTccc                                <  1:2164358/98‑1 (MQ=255)
                                                                    gCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCgg                              >  1:138214/1‑100 (MQ=255)
                                                                    gCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCg                               >  1:925368/1‑99 (MQ=255)
                                                                    gCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTCGCTTCATAAGCATCACGCAACAGCTGCAAAACAGCATCTTTCCCgg                              >  1:650122/1‑100 (MQ=255)
                                                                       cTAGGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCAGCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAAc                           <  1:2302544/100‑1 (MQ=255)
                                                                         aTGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCa                         <  1:2599295/100‑1 (MQ=255)
                                                                         aTGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCa                         >  2:2170039/1‑100 (MQ=255)
                                                                            aCTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCa                      <  1:101053/100‑1 (MQ=255)
                                                                             cTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCAt                     <  1:2028121/100‑1 (MQ=255)
                                                                             cTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCAt                     >  2:2410377/1‑100 (MQ=255)
                                                                              tCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATc                    <  1:1230501/100‑1 (MQ=255)
                                                                               cgagggaagaagGTTGAGTCAACCTATACAATCTGGAGGACCTACTGCATAAGCATCCCTCAAAAGCTGCAAAACACCATATTTCCCTGAACCa                         <  1:2541834/89‑1 (MQ=255)
                                                                               ccAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCa                   <  1:1028131/100‑1 (MQ=255)
                                                                               ccAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCa                   <  2:432281/100‑1 (MQ=255)
                                                                                cAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCaa                  >  2:2283225/1‑100 (MQ=255)
                                                                                cAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCaa                  >  1:77411/1‑100 (MQ=255)
                                                                                cAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGAATCaa                  >  1:2843897/1‑100 (MQ=255)
                                                                                cAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCaa                  >  1:112911/1‑100 (MQ=255)
                                                                                 aaaGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAg                 <  1:2283210/100‑1 (MQ=255)
                                                                                  aagaagaagGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGa                <  2:773591/100‑1 (MQ=255)
                                                                                  aagaagaagGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGa                <  1:176681/100‑1 (MQ=255)
                                                                                  aagaagaagGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTCTCCCGGAACCAGCATCAAGa                <  2:1844564/100‑1 (MQ=255)
                                                                                    gaagaagGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAAc              <  2:2744492/100‑1 (MQ=255)
                                                                                      agaagGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTc            >  2:1652590/1‑100 (MQ=255)
                                                                                      agaagGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTc            <  1:1666624/100‑1 (MQ=255)
                                                                                      agaagGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCAAGCAAAAACTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTc            >  2:1266952/1‑100 (MQ=255)
                                                                                       gaagGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCg           >  2:2613632/1‑100 (MQ=255)
                                                                                        aagGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGc          <  2:1925892/100‑1 (MQ=255)
                                                                                        aagGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGc          <  2:1164430/100‑1 (MQ=255)
                                                                                         agGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGcc         <  2:1344869/100‑1 (MQ=255)
                                                                                          gggTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCg        >  1:2614126/1‑100 (MQ=255)
                                                                                          gggTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCg        >  2:2250298/1‑100 (MQ=255)
                                                                                            gTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGtt      >  2:1632495/1‑100 (MQ=255)
                                                                                            gTTAGCCAACCGATACAATTTTGCGTACTTGCCTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACGCGCCGtt      <  1:140955/100‑1 (MQ=255)
                                                                                                gCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTCGct  >  2:2818658/1‑100 (MQ=255)
                                                                                                 |                                                                                                  
TACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTCGCT  >  NC_000913/2856791‑2856985

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: