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.

AGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTCGCTT  >  NC_000913/2856788‑2856986
                                                                                                    |                                                                                                   
 gaTACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGcc                                                                                                    >  2:1272455/1‑100 (MQ=255)
   tACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCaa                                                                                                  >  2:397851/1‑100 (MQ=255)
     cGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAAcc                                                                                                >  1:992016/1‑100 (MQ=255)
     cGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAAcc                                                                                                <  2:3760423/100‑1 (MQ=255)
     cGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGACAAcc                                                                                                >  1:2814451/1‑100 (MQ=255)
     cGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGC‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAAcc                                                                                                >  1:2356933/1‑100 (MQ=255)
       cTTGCCGGACATGCTGCCTCCACCTCATTAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGa                                                                                              <  2:1293015/100‑1 (MQ=255)
       cTTGCCGGACATGCTGCCTCCACCTCATTAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGa                                                                                              >  2:2028985/1‑100 (MQ=255)
         tGCCGGACATGCTGCCTCCACCTCATTAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATa                                                                                            <  2:3895569/100‑1 (MQ=255)
         tGCCGGACATGCTGCCTCCACCTCATTAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATa                                                                                            >  2:3288288/1‑100 (MQ=255)
          gCCGGACATGCTGCCTCCACCTCATTAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATAc                                                                                           >  2:1228722/1‑100 (MQ=255)
            cGGACATGCTGCCTCCACCTCATTAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACaa                                                                                         >  2:1360424/1‑100 (MQ=255)
               aCATGCTGCCTCCACCTCATTAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAAttt                                                                                      <  1:1335/100‑1 (MQ=255)
                cATGCTGCCTCCACCTCATTAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAAtttt                                                                                     >  2:3638123/1‑100 (MQ=255)
                  tgctgcCTCCACCTCATTAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGc                                                                                   >  1:563922/1‑100 (MQ=255)
                  tgctgcCTCCACCTCATTAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTg                                                                                    >  1:3105759/1‑99 (MQ=255)
                  tgctgcCTCCACCTCATTAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGc                                                                                   <  2:2111005/100‑1 (MQ=255)
                   gctgcCTCCACCTCATTAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACGCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCg                                                                                  <  1:2730076/100‑1 (MQ=255)
                    ctgcCTCCACCTCATTAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGt                                                                                 >  2:2960822/1‑100 (MQ=255)
                    ctgcCTCCACCTCATTAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGt                                                                                 <  2:3060556/100‑1 (MQ=255)
                    ctgcCTCCACCTCATTAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGATAGTCAACCGATACAATTTTGCGt                                                                                 >  2:1775203/1‑100 (MQ=255)
                      gcCTCCACCTCATTAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTAc                                                                               <  2:2005345/100‑1 (MQ=255)
                      gcCTCCACCTCATTAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTAc                                                                               >  1:4357027/1‑100 (MQ=255)
                      gcCTCCACCTCATTAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTAc                                                                               >  2:2787843/1‑100 (MQ=255)
                         tCCACCTCATTAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTg                                                                            <  2:4146808/100‑1 (MQ=255)
                         tCCACCTCATTAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTg                                                                            >  1:3382422/1‑100 (MQ=255)
                          ccACCTCATTAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTgt                                                                           <  1:38835/100‑2 (MQ=255)
                          ccACCTCATTAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGc                                                                           <  2:3369520/100‑1 (MQ=255)
                            aCCTCATTAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCtt                                                                         >  1:2287201/1‑100 (MQ=255)
                             ccTCATTAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTc                                                                        <  2:1364654/100‑1 (MQ=255)
                              cTCATTAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCa                                                                       >  2:4358455/1‑100 (MQ=255)
                              cTCATTAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAAGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCa                                                                       >  1:1269663/1‑100 (MQ=255)
                               tCATTAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCGAAAGAAGAAGGGTTAGTCAACAGATACAATTTTACGTACTTGCTTCAt                                                                      >  1:2604654/1‑100 (MQ=255)
                               agaTTAAGATGTAATATAAA‑TACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCAt                                                                      >  1:4292633/3‑100 (MQ=37)
                                gaTTAAGATGT‑AATATAAATACATCTTAATCTTAAAGGACACTATGACTCCAAAGAAGAAGGGTTAGCCAACAGATACAATTCTTCGAACTCGc                                                                           >  1:1534324/2‑94 (MQ=21)
                                cATTAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATa                                                                     <  2:2152964/100‑1 (MQ=255)
                                cATTAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATa                                                                     <  2:2703299/100‑1 (MQ=255)
                                cATTAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATa                                                                     <  2:2247303/100‑1 (MQ=255)
                                 aTTAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATaa                                                                    >  1:444915/1‑100 (MQ=255)
                                 aTTAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATaa                                                                    <  2:4253030/100‑1 (MQ=255)
                                  ttAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAg                                                                   <  1:2865300/100‑1 (MQ=255)
                                   tAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGc                                                                  <  2:3734097/100‑1 (MQ=255)
                                   tAAGATGT‑ATTTATATGACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAg                                                                   <  2:719710/99‑1 (MQ=255)
                                    aaGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCa                                                                 <  1:1669980/100‑1 (MQ=255)
                                       aTGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCa                                                              >  2:408245/1‑100 (MQ=255)
                                         gT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACg                                                            <  1:3514573/100‑1 (MQ=255)
                                         gT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACg                                                            <  2:2599269/100‑1 (MQ=255)
                                          t‑ATTTATATTACATCTTAATCTTAAAGGGCACTCTGACTCCAAAGAACTAGGGTTAGCCGACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGc                                                           <  1:139024/100‑1 (MQ=255)
                                          t‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGc                                                           >  1:1009716/1‑100 (MQ=255)
                                          t‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGc                                                           >  2:1008062/1‑100 (MQ=255)
                                          t‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGc                                                           <  2:3709606/100‑1 (MQ=255)
                                          t‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGc                                                           >  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)
                                                                                                    |                                                                                                   
AGATACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGT‑ATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTCGCTT  >  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: ‑