Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,856,887 C→T 100% 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,8870CT90.2% 341.1 / 18.0 126intergenic (‑81/‑206)ygbA/mutSuncharacterized protein/methyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base C (5/8);  new base T (55/58);  total (60/66)
Fisher's exact test for biased strand distribution p-value = 5.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 126 total aligned reads displayed.

ACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTC  >  NC_000913/2856792‑2856982
                                                                                               |                                                                                               
aCGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAAc                                                                                             >  1:2441270/1‑100 (MQ=255)
  gCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCg                                                                                           <  1:501955/100‑1 (MQ=255)
         ggACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAAt                                                                                    >  1:2399923/1‑100 (MQ=255)
         ggACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAAt                                                                                    >  1:2400351/1‑100 (MQ=255)
            cATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAAtttt                                                                                 <  1:440084/100‑1 (MQ=255)
            cATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAAtttt                                                                                 >  2:10913/1‑100 (MQ=255)
             aTGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTg                                                                                <  2:68820/100‑1 (MQ=255)
             aTGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTg                                                                                <  2:1630676/100‑1 (MQ=255)
              tgctgcCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGTGTTAGTCAACCGATACAATTTTGc                                                                               >  2:1263844/1‑100 (MQ=255)
               gctgcCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCg                                                                              >  2:77402/1‑100 (MQ=255)
               gctgcCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCg                                                                              <  1:28838/100‑1 (MQ=255)
                  gcCTCCACCTCATTAAGCTGTATTTATATTACATCTTAATCTCAACGGGCACTATAACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTAc                                                                           >  2:2326281/1‑100 (MQ=255)
                  gcCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTan                                                                           <  1:1101982/100‑2 (MQ=255)
                  gcCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTAc                                                                           >  2:1847187/1‑100 (MQ=255)
                  gcCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTAc                                                                           <  1:1406610/100‑1 (MQ=255)
                   ccTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACt                                                                          >  2:2185034/1‑100 (MQ=255)
                    cTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACtt                                                                         >  2:1953229/1‑100 (MQ=255)
                    cTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACtt                                                                         >  2:1229446/1‑100 (MQ=255)
                     tCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTg                                                                        >  1:1863953/1‑100 (MQ=255)
                     tCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTg                                                                        >  1:1298514/1‑100 (MQ=255)
                      ccACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGc                                                                       >  2:1443539/1‑100 (MQ=255)
                       cACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCt                                                                      <  2:781976/100‑1 (MQ=255)
                        aCCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCtt                                                                     <  1:2153808/100‑1 (MQ=255)
                         ccTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTc                                                                    <  1:1221203/100‑1 (MQ=255)
                           tCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCAt                                                                  >  1:691208/1‑100 (MQ=255)
                            cATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATa                                                                 >  1:876042/1‑100 (MQ=255)
                            cATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATa                                                                 >  1:1077142/1‑100 (MQ=255)
                            cATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATa                                                                 <  2:857392/100‑1 (MQ=255)
                            cATTAAGATGTATTTATATTACATCTTAATCTCAAAGAGGAATATGATTTAAAAGAAGAAGGGTTAGTTAACTGATAAAATTTTGCATACTTGCGTtata                                                                 >  2:1408408/1‑100 (MQ=255)
                             aTTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATaa                                                                <  2:1365635/100‑1 (MQ=255)
                             aTTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATaa                                                                <  2:765443/100‑1 (MQ=255)
                              ttAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAg                                                               <  1:1631096/100‑1 (MQ=255)
                              ttAAGATGTATTTATAATACATATTAATCTGAAAGGGAAATATGACTCCAAAGTAGAAGGGTCAATAAACGAATACAATTATGCATAGTTGCtgacta                                                                 >  2:1550129/1‑93 (MQ=255)
                               tAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGc                                                              <  1:1915879/100‑1 (MQ=255)
                                aaGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCa                                                             <  1:1373180/100‑1 (MQ=255)
                                 aGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCAt                                                            <  1:2376302/100‑1 (MQ=255)
                                 aGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCAt                                                            >  1:252926/1‑100 (MQ=255)
                                  gATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATc                                                           <  1:199952/100‑1 (MQ=255)
                                   aTGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTACAAAGAATAAGTGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAgaacc                                                           >  2:2208834/1‑95 (MQ=255)
                                    tGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCGTAAGCATCAc                                                         >  1:486952/1‑100 (MQ=255)
                                      tATTTATATTACATCTTAATCTTAAAGGGCACTATTACTCCAAAGACGAAGATTTAATAAACCGAAACAATTTTTCGTACTTGATTCATAAGCAACACGc                                                       >  2:856435/1‑100 (MQ=255)
                                      tATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGTTTAGTCAACCGCTACAATTTTACGTACTTGTATCATAAGCATCACGc                                                       >  2:856353/1‑100 (MQ=255)
                                      tATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGc                                                       >  2:1597549/1‑100 (MQ=255)
                                       aTTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGAGTACTTGCTTCATAAGCATCACGCa                                                      <  2:424771/100‑1 (MQ=255)
                                         ttATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCaaa                                                    <  1:1597558/100‑1 (MQ=255)
                                         ttATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCaa                                                     >  2:1652789/1‑99 (MQ=255)
                                          tataTTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCACCACGCaaaa                                                   >  1:1434480/1‑100 (MQ=255)
                                             atTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGACGAAGGGTTAATCAACCCCTACAATTTTGAGTACTTGCTTCATAAGACTCACGCaaa                                                    >  2:856760/1‑96 (MQ=255)
                                             atTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCt                                                <  2:2209847/100‑1 (MQ=255)
                                             atTACATCTTAATCTTAAAGGGCACTATGACTACAAAGAAGAAGGGTTAGTCAAACGATACAATTTTGCGTACTTGCGTCATAAGCATCACGCAAAATCt                                                >  2:1287066/1‑100 (MQ=255)
                                                aCATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCa                                             <  2:412728/100‑1 (MQ=255)
                                                aCATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCa                                             <  2:92907/100‑1 (MQ=255)
                                                  atcttaatcttaAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCaaa                                           >  1:1884216/1‑100 (MQ=255)
                                                    cttaatcttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAAc                                         >  1:442109/1‑100 (MQ=255)
                                                     ttaatcttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACa                                        >  2:2297767/1‑100 (MQ=255)
                                                      taatcttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAACACAg                                       >  2:404337/1‑100 (MQ=255)
                                                      taatcttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAg                                       >  1:769262/1‑100 (MQ=255)
                                                       aatcttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAAc                                         <  2:2523307/97‑1 (MQ=255)
                                                        atcttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCa                                     >  1:1241252/1‑100 (MQ=255)
                                                         tcttaaAGCGCACTATGACTCCAAATAAGAAGCGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCAt                                    <  2:2465362/100‑1 (MQ=255)
                                                          cttaaAGGGCACTATTACTCCAAAGAAGAAGGGATAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCTAAACAGCATc                                   >  2:1715728/1‑100 (MQ=255)
                                                          cttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATc                                   <  1:702817/100‑1 (MQ=255)
                                                          cttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGCTACACTTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGc                                      >  2:2066736/1‑97 (MQ=255)
                                                           ttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCt                                  >  1:1594470/1‑100 (MQ=255)
                                                            taaAGGGCAGTATGACTCCAAGTAAGAAGGGTTAGCGAACCGATACAATTTTGCGTACTGGCTTCAGAAGCATCACGCAAAAGCTGCAAAACAGCa                                     <  2:517363/96‑1 (MQ=255)
                                                              aaGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTc                               <  1:1017936/100‑1 (MQ=255)
                                                               aGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTcc                              <  2:627035/100‑1 (MQ=255)
                                                               aGGGCACTATGAACCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTcc                              <  1:2319276/100‑1 (MQ=255)
                                                                gggCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTccc                             >  1:1973041/1‑100 (MQ=255)
                                                                 ggCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCg                            <  1:898694/100‑1 (MQ=255)
                                                                 ggCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCg                            >  2:1106004/1‑100 (MQ=255)
                                                                    aCTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGaa                         >  2:901899/1‑100 (MQ=255)
                                                                     cTATGAGTCCAAAGAAGAAGGGTTAGTCAACCGATACAAGTTTGCTTACTTGCTTCATAAGCATCACGCAAAATCTGCAAAACAGCATCTTTCCCGGAAc                        <  1:1517778/100‑1 (MQ=255)
                                                                     cTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAAc                        >  2:1628601/1‑100 (MQ=255)
                                                                      tATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAAcc                       <  2:1653502/100‑1 (MQ=255)
                                                                      tATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAAcc                       <  1:411810/100‑1 (MQ=255)
                                                                       aTGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCa                      <  2:1973031/100‑1 (MQ=255)
                                                                        tgccTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAg                     <  2:2410680/97‑1 (MQ=255)
                                                                        tGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAg                     >  2:2455613/1‑100 (MQ=255)
                                                                          aCTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCa                   >  1:2303803/1‑100 (MQ=255)
                                                                            tCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATc                 <  2:1982589/100‑1 (MQ=255)
                                                                            tCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATc                 <  1:249380/100‑1 (MQ=255)
                                                                            tCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATc                 <  2:1269397/100‑1 (MQ=255)
                                                                             ccAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCa                <  1:206969/100‑1 (MQ=255)
                                                                             ccAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCa                <  1:2458199/100‑1 (MQ=255)
                                                                             ccAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCa                <  2:2301468/100‑1 (MQ=255)
                                                                              cAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCaa               >  2:1938519/1‑100 (MQ=255)
                                                                               aaaGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAg              <  1:1737272/100‑1 (MQ=255)
                                                                                aagaagaagGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAg              <  2:1049253/99‑1 (MQ=255)
                                                                                aagaagaagGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGa             <  1:916953/100‑1 (MQ=255)
                                                                                aagaagaagGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGa             <  2:590174/100‑1 (MQ=255)
                                                                                 agaagaagGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGaa            >  1:438759/1‑100 (MQ=255)
                                                                                  gaagaagGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGaa            <  2:807677/99‑1 (MQ=255)
                                                                                     gaagGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCg        <  1:2519210/100‑1 (MQ=255)
                                                                                      aagGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGc       <  2:946518/100‑1 (MQ=255)
                                                                                       agGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGcc      <  2:1338520/100‑1 (MQ=255)
                                                                                       agGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGcc      <  1:1436194/100‑1 (MQ=255)
                                                                                         ggTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGt    >  1:1225252/1‑100 (MQ=255)
                                                                                          gTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGtt   >  1:1085619/1‑100 (MQ=255)
                                                                                           ttAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTc  <  2:653260/100‑1 (MQ=255)
                                                                                               |                                                                                               
ACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTC  >  NC_000913/2856792‑2856982

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: