Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,856,887 C→T 32.3% 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,8870CT32.3% 206.6 / 140.3 168intergenic (‑81/‑206)ygbA/mutSuncharacterized protein/methyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base C (51/62);  new base T (29/25);  total (81/87)
Fisher's exact test for biased strand distribution p-value = 3.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.29e-01

Only 100 of 172 total aligned reads displayed.

AGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTCGCTT  >  NC_000913/2856788‑2856986
                                                                                                    |                                                                                                   
 gaTACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGcc                                                                                                    >  2:1272455/1‑100 (MQ=255)
   tACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCaa                                                                                                  >  2:397851/1‑100 (MQ=255)
     cGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAAcc                                                                                                >  1:992016/1‑100 (MQ=255)
     cGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAAcc                                                                                                <  2:3760423/100‑1 (MQ=255)
     cGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGACAAcc                                                                                                >  1:2814451/1‑100 (MQ=255)
     cGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGCATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAAcc                                                                                                >  1:2356933/1‑100 (MQ=255)
       cTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGa                                                                                              <  2:1293015/100‑1 (MQ=255)
       cTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGa                                                                                              >  2:2028985/1‑100 (MQ=255)
         tGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATa                                                                                            <  2:3895569/100‑1 (MQ=255)
         tGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATa                                                                                            >  2:3288288/1‑100 (MQ=255)
          gCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATAc                                                                                           >  2:1228722/1‑100 (MQ=255)
            cGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACaa                                                                                         >  2:1360424/1‑100 (MQ=255)
               aCATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAAttt                                                                                      <  1:1335/100‑1 (MQ=255)
                cATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAAtttt                                                                                     >  2:3638123/1‑100 (MQ=255)
                  tgctgcCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGc                                                                                   >  1:563922/1‑100 (MQ=255)
                  tgctgcCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTg                                                                                    >  1:3105759/1‑99 (MQ=255)
                  tgctgcCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGc                                                                                   <  2:2111005/100‑1 (MQ=255)
                   gctgcCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACGCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCg                                                                                  <  1:2730076/100‑1 (MQ=255)
                    ctgcCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGt                                                                                 >  2:2960822/1‑100 (MQ=255)
                    ctgcCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGt                                                                                 <  2:3060556/100‑1 (MQ=255)
                    ctgcCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGATAGTCAACCGATACAATTTTGCGt                                                                                 >  2:1775203/1‑100 (MQ=255)
                      gcCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTAc                                                                               <  2:2005345/100‑1 (MQ=255)
                      gcCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTAc                                                                               >  1:4357027/1‑100 (MQ=255)
                      gcCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTAc                                                                               >  2:2787843/1‑100 (MQ=255)
                         tCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTg                                                                            <  2:4146808/100‑1 (MQ=255)
                         tCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTg                                                                            >  1:3382422/1‑100 (MQ=255)
                          ccACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTgt                                                                           <  1:38835/100‑2 (MQ=255)
                          ccACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGc                                                                           <  2:3369520/100‑1 (MQ=255)
                            aCCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCtt                                                                         >  1:2287201/1‑100 (MQ=255)
                             ccTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTc                                                                        <  2:1364654/100‑1 (MQ=255)
                              cTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCa                                                                       >  2:4358455/1‑100 (MQ=255)
                              cTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAAGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCa                                                                       >  1:1269663/1‑100 (MQ=255)
                               tCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCGAAAGAAGAAGGGTTAGTCAACAGATACAATTTTACGTACTTGCTTCAt                                                                      >  1:2604654/1‑100 (MQ=255)
                               agaTTAAGATGTAATATAAATACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCAt                                                                      >  1:4292633/3‑100 (MQ=37)
                                gaTTAAGATGTAATATAAATACATCTTAATCTTAAAGGACACTATGACTCCAAAGAAGAAGGGTTAGCCAACAGATACAATTCTTCGAACTCGc                                                                           >  1:1534324/2‑94 (MQ=21)
                                cATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATa                                                                     <  2:2152964/100‑1 (MQ=255)
                                cATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATa                                                                     <  2:2703299/100‑1 (MQ=255)
                                cATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATa                                                                     <  2:2247303/100‑1 (MQ=255)
                                 aTTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATaa                                                                    >  1:444915/1‑100 (MQ=255)
                                 aTTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATaa                                                                    <  2:4253030/100‑1 (MQ=255)
                                  ttAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAg                                                                   <  1:2865300/100‑1 (MQ=255)
                                   tAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGc                                                                  <  2:3734097/100‑1 (MQ=255)
                                   tAAGATGTATTTATATGACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAg                                                                   <  2:719710/99‑1 (MQ=255)
                                    aaGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCa                                                                 <  1:1669980/100‑1 (MQ=255)
                                       aTGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCa                                                              >  2:408245/1‑100 (MQ=255)
                                         gTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACg                                                            <  1:3514573/100‑1 (MQ=255)
                                         gTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACg                                                            <  2:2599269/100‑1 (MQ=255)
                                          tATTTATATTACATCTTAATCTTAAAGGGCACTCTGACTCCAAAGAACTAGGGTTAGCCGACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGc                                                           <  1:139024/100‑1 (MQ=255)
                                          tATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGc                                                           >  1:1009716/1‑100 (MQ=255)
                                          tATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGc                                                           >  2:1008062/1‑100 (MQ=255)
                                          tATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGc                                                           <  2:3709606/100‑1 (MQ=255)
                                          tATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGc                                                           >  2:2050949/1‑100 (MQ=255)
                                            aTTTATATTACATCGTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCa                                                          <  1:1442295/100‑1 (MQ=255)
                                             tttATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCaa                                                         >  1:1114241/1‑100 (MQ=255)
                                              ttATATTACCTCTTACTCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCaaa                                                        <  1:37787/100‑1 (MQ=255)
                                              ttATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCaaa                                                        <  2:1210041/100‑1 (MQ=255)
                                               tataTTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCaaaa                                                       >  1:430499/1‑100 (MQ=255)
                                                atatTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAg                                                      <  1:3174148/100‑1 (MQ=255)
                                                 tatTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGTGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGc                                                     <  1:1544871/100‑1 (MQ=255)
                                                 tatTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGc                                                     <  2:664118/100‑1 (MQ=255)
                                                 tatTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGc                                                     <  1:2634472/100‑1 (MQ=255)
                                                   ttACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTg                                                   <  2:910915/100‑1 (MQ=255)
                                                   ttACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTg                                                   <  2:804445/100‑1 (MQ=255)
                                                      cATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCaa                                                <  1:3852432/100‑1 (MQ=255)
                                                       atcttaatcttaAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGATTCATAAGCATCACGCAAAAGGTGCaaa                                               >  2:3625720/1‑100 (MQ=255)
                                                        tcttaatcttaaAGGGCACTATGACTCCAAAGAAGGAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCaaaa                                              >  2:1090184/1‑100 (MQ=255)
                                                        tcttaatcttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCaaaa                                              >  1:4055331/1‑100 (MQ=255)
                                                         cttaatcttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAAc                                             <  2:182455/100‑1 (MQ=255)
                                                           taatcttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAg                                           <  2:3922663/100‑1 (MQ=255)
                                                           taatcttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAg                                           <  2:4301478/100‑1 (MQ=255)
                                                           taatcttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAg                                           <  1:1350037/100‑1 (MQ=255)
                                                            aatcttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGc                                          >  1:1924183/1‑100 (MQ=255)
                                                            aatcttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGc                                          >  2:3131053/1‑100 (MQ=255)
                                                               cttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATc                                       >  1:7341/1‑100 (MQ=255)
                                                                   aaGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTc                                   <  2:444916/100‑1 (MQ=255)
                                                                   aaGGGCACTATGACTCCAAAGAAGAAGGGTTAACCAACCGATACAATTTTGTGTACTTGCTTCATAAGCATCACGCAAAAGTTGCAAAACAGCATCTTTc                                   >  1:3045611/1‑100 (MQ=255)
                                                                      ggCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCg                                >  1:3601951/1‑100 (MQ=255)
                                                                      ggCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCg                                <  1:1869151/100‑1 (MQ=255)
                                                                             tGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAg                         <  1:2627273/100‑1 (MQ=255)
                                                                               aCTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCa                       <  2:2702525/100‑1 (MQ=255)
                                                                               aCTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCa                       <  2:3057527/100‑1 (MQ=255)
                                                                                cTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCAt                      >  2:3519293/1‑100 (MQ=255)
                                                                                cTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCAt                      >  1:1360819/1‑100 (MQ=255)
                                                                                 tCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATc                     >  1:3790605/1‑100 (MQ=255)
                                                                                  ccAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAgagatc                    >  1:2069174/1‑95 (MQ=255)
                                                                                  ccAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCa                    >  1:3127390/1‑100 (MQ=255)
                                                                                  ccAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCa                    >  1:1273680/1‑100 (MQ=255)
                                                                                  ccAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTATCCCGGAACCAGCATCa                    >  1:823539/1‑100 (MQ=255)
                                                                                  ccAAAAAAGAAGGGTTAGCCAACATATAGAATTTTGCGTACATGGTTAATAAGTATCACGAAAAAGCTTCAAAACAGCATCTTTAACGTAACCAga                        >  1:3645553/1‑95 (MQ=255)
                                                                                   cAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCaa                   >  1:4289574/1‑100 (MQ=255)
                                                                                     aagaagaagGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGa                 <  1:2857530/100‑1 (MQ=255)
                                                                                     aagaagaagGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCTGGAACCAGCATCAAGa                 <  2:1962290/100‑1 (MQ=255)
                                                                                     aagaagaagGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGa                 <  1:1165564/100‑1 (MQ=255)
                                                                                         agaagGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTc             <  1:2312731/100‑1 (MQ=255)
                                                                                         agaagGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTc             >  2:1762558/1‑100 (MQ=255)
                                                                                           aagGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGc           <  1:1987180/100‑1 (MQ=255)
                                                                                            agGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCACCAAGAACTCGcc          >  1:3942699/1‑100 (MQ=255)
                                                                                              ggTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGt        >  2:2192720/1‑100 (MQ=255)
                                                                                                ttAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTc      <  1:3807542/100‑1 (MQ=255)
                                                                                                 tAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTCg     >  2:3207772/1‑100 (MQ=255)
                                                                                                    |                                                                                                   
AGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTCGCTT  >  NC_000913/2856788‑2856986

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