Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,856,887 C→T 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,8870CT98.9% 309.0 / ‑4.2 90intergenic (‑81/‑206)ygbA/mutSuncharacterized protein/methyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (51/38);  total (52/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

ACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGT  >  NC_000913/2856792‑2856980
                                                                                               |                                                                                             
aCGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAAc                                                                                           <  2:1585793/100‑1 (MQ=255)
 cGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAAcc                                                                                          <  2:937557/100‑1 (MQ=255)
 cGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAAcc                                                                                          <  1:2471765/100‑1 (MQ=255)
  gCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCg                                                                                         <  2:1959812/100‑1 (MQ=255)
   cTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGa                                                                                        >  1:2371266/1‑100 (MQ=255)
   cTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGa                                                                                        <  2:2134732/100‑1 (MQ=255)
     tGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATa                                                                                      <  2:820113/100‑1 (MQ=255)
      gCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATAc                                                                                     <  1:2681365/100‑1 (MQ=255)
             aTGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTg                                                                              <  1:2708706/100‑1 (MQ=255)
             aTGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTg                                                                              >  2:2125494/1‑100 (MQ=255)
              tgctgcCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGc                                                                             >  1:1735812/1‑100 (MQ=255)
               gctgcCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCg                                                                            <  1:914017/100‑1 (MQ=255)
               gctgcCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCg                                                                            <  2:2223607/100‑1 (MQ=255)
               gctgcCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCg                                                                            <  2:2154876/100‑1 (MQ=255)
                ctgcCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGt                                                                           <  2:2168624/100‑1 (MQ=255)
                  nccTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTAc                                                                         >  2:1016143/2‑100 (MQ=255)
                   ccTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACt                                                                        >  1:501278/1‑100 (MQ=255)
                      ccACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGc                                                                     >  2:2191584/1‑100 (MQ=255)
                      ccACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGc                                                                     <  2:1836092/100‑1 (MQ=255)
                       cACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCt                                                                    <  2:530722/100‑1 (MQ=255)
                       cACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCt                                                                    >  2:1672094/1‑100 (MQ=255)
                       cACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCt                                                                    <  1:2372051/100‑1 (MQ=255)
                        aCCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCtt                                                                   >  2:2661190/1‑100 (MQ=255)
                         ccTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTc                                                                  <  2:841273/100‑1 (MQ=255)
                          cTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCa                                                                 >  1:2685263/1‑100 (MQ=255)
                          cTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCa                                                                 >  1:483566/1‑100 (MQ=255)
                            cATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTTTGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATa                                                               <  2:2306464/100‑1 (MQ=255)
                             aTTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATaa                                                              <  2:1847391/100‑1 (MQ=255)
                              ttAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAGGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAg                                                             <  1:2364104/100‑1 (MQ=255)
                               tAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATTCAATTTTGCGTACTTGCTTCATAAGc                                                            <  1:1680249/100‑1 (MQ=255)
                                ccgatctATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCa                                                           <  1:63420/94‑1 (MQ=255)
                                aaGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCa                                                           <  1:1096985/100‑1 (MQ=255)
                                aaGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCa                                                           >  2:678959/1‑100 (MQ=255)
                                aaGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCAAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTACTTCATAAGCa                                                           <  2:375417/100‑1 (MQ=255)
                                  gATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATc                                                         <  2:474547/100‑1 (MQ=255)
                                  gATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATc                                                         >  2:2339141/1‑100 (MQ=255)
                                    tGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCAc                                                       <  1:1350350/100‑1 (MQ=255)
                                     gTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACg                                                      >  2:1786281/1‑100 (MQ=255)
                                      tATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGc                                                     >  2:1257085/1‑100 (MQ=255)
                                      tATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGc                                                     >  1:1983713/1‑100 (MQ=255)
                                       aTTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTTAACTGATACAATTTTGCGTACTTGCTTCATAAGCATCAcact                                                    >  2:2382821/1‑97 (MQ=255)
                                       aTTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCaagatcgg                                                    >  2:63419/1‑93 (MQ=255)
                                        tttATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCaa                                                   >  1:2234665/1‑100 (MQ=255)
                                           atatTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAg                                                >  2:1913480/1‑100 (MQ=255)
                                            tatTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGc                                               >  1:2238957/1‑100 (MQ=255)
                                             atTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCt                                              >  2:1523729/1‑100 (MQ=255)
                                             atTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCt                                              >  1:728286/1‑100 (MQ=255)
                                              ttACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTg                                             >  1:707001/1‑100 (MQ=255)
                                               tACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGc                                            >  1:340733/1‑100 (MQ=255)
                                                  atcttaatcttaAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCaaa                                         >  1:2026237/1‑100 (MQ=255)
                                                   tcttaatcttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCaaaa                                        >  1:1120298/1‑100 (MQ=255)
                                                   tcttaatcttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCaaaa                                        >  2:2252101/1‑100 (MQ=255)
                                                   tcttaatcttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCaaaa                                        >  1:1691028/1‑100 (MQ=255)
                                                    cttaatcttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCaaaa                                        >  2:44662/1‑99 (MQ=255)
                                                    cttaatcttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAAc                                       >  1:950754/1‑100 (MQ=255)
                                                    cttaatcttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAAc                                       <  1:786320/100‑1 (MQ=255)
                                                       aatcttaaAGGGCACTATGACTCCAAAGAAAAAAGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCAGCAAAACAGc                                    >  2:1811814/1‑100 (MQ=255)
                                                        ttaataaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCa                                   >  1:441794/5‑100 (MQ=255)
                                                         tcttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCAt                                  >  1:2603419/1‑100 (MQ=255)
                                                          cttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATc                                 >  2:454439/1‑100 (MQ=255)
                                                          cttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATc                                 >  2:258185/1‑100 (MQ=255)
                                                          cttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATc                                 >  2:2233918/1‑100 (MQ=255)
                                                           ttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATAGAATTTTGCGTACTTGCTTCATGAGCATCACGCAAAAGCTGCAAAACAGCAt                                  >  1:549641/1‑98 (MQ=255)
                                                            naaaGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCtt                               >  2:325999/2‑100 (MQ=255)
                                                              aaGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTc                             <  1:334697/100‑1 (MQ=255)
                                                               aGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTcc                            <  1:44663/100‑1 (MQ=255)
                                                                 ggCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCg                          <  2:2685243/100‑1 (MQ=255)
                                                                 ggCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCg                          <  1:1000823/100‑1 (MQ=255)
                                                                 ggCACTATGAATCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCg                          <  1:2318175/100‑1 (MQ=255)
                                                                 cgatcTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCg                          <  1:2040414/96‑1 (MQ=255)
                                                                   cACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGa                        >  2:232013/1‑100 (MQ=255)
                                                                   cACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGa                        >  2:256703/1‑100 (MQ=255)
                                                                    aCTATGACGGCAAAGCAGAAGGGTTACTCACCAGATACAGTTTTGCGAACTTTCTTCATAAGCATCACGCAAAATCTGCAAGACATCATCtttggcg                          >  2:675900/1‑93 (MQ=255)
                                                                     cTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAAc                      >  1:1639370/1‑100 (MQ=255)
                                                                      tATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAAcc                     <  2:340730/100‑1 (MQ=255)
                                                                       aTGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCgagatcg                    >  2:2040397/1‑94 (MQ=255)
                                                                        tGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAg                   >  1:1985659/1‑100 (MQ=255)
                                                                         gACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGTATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGc                  >  2:1934454/1‑100 (MQ=255)
                                                                           cTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCAt                >  1:1375957/1‑100 (MQ=255)
                                                                            tCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATc               <  2:2371247/100‑1 (MQ=255)
                                                                               aaaGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAg            >  1:682966/1‑100 (MQ=255)
                                                                                aagaagaagGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGa           <  1:2006083/100‑1 (MQ=255)
                                                                                aagaagaagGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGa           >  2:1452592/1‑100 (MQ=255)
                                                                                    agaagGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTc       <  2:2212252/100‑1 (MQ=255)
                                                                                    agaagGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATAACGCAAAAGCTGCAAACCAGCATCTTTCCCGGAACCAGCATCAAGAACt        >  1:227042/1‑99 (MQ=255)
                                                                                        gggTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCg   >  2:2001865/1‑100 (MQ=255)
                                                                                        gggTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCg   >  1:1112811/1‑100 (MQ=255)
                                                                                        gggTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCg   <  2:298108/100‑1 (MQ=255)
                                                                                         ggTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGt  <  1:933619/100‑1 (MQ=255)
                                                                                         ggTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGt  >  1:354601/1‑100 (MQ=255)
                                                                                               |                                                                                             
ACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGT  >  NC_000913/2856792‑2856980

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