Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,041,028 G→A 100% R15C (CGT→TGT)  sdhE ← flavinator of succinate dehydrogenase; antitoxin of CptAB toxin‑antitoxin pair

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,041,0280GA97.7% 266.8 / ‑2.9 85R15C (CGT→TGT) sdhEflavinator of succinate dehydrogenase; antitoxin of CptAB toxin‑antitoxin pair
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (34/49);  total (34/51)
Fisher's exact test for biased strand distribution p-value = 5.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

CAGACGAATAAAGATGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACGGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATCCAGTACCCGGCTATTGTAACGTGTTTTTCGC  >  NC_000913/3040933‑3041124
                                                                                               |                                                                                                
nagACGAATAAAGATGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCAcagcgg                                                                                              >  2:2408933/2‑100 (MQ=255)
cagACGAATAAAGATGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCAcagcgg                                                                                              <  1:2167086/100‑1 (MQ=255)
 agACGAATAAAGATGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACagcggc                                                                                             <  1:1919863/100‑1 (MQ=255)
 agACGAATAAAGATGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACagcggc                                                                                             <  2:137579/100‑1 (MQ=255)
  gACGAATAAAGATGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCa                                                                                            >  2:535718/1‑100 (MQ=255)
  gACGAATAAAGATGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCa                                                                                            <  1:1920410/100‑1 (MQ=255)
  gACGAATAAAGATGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCa                                                                                            >  1:1916649/1‑100 (MQ=255)
   aCGAATAAAGATGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCAt                                                                                           >  1:1236311/1‑100 (MQ=255)
    cGAATAAAGATGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATg                                                                                          >  1:2489951/1‑100 (MQ=255)
    cGAATAAAGATGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATg                                                                                          >  1:1668028/1‑100 (MQ=255)
      aaTAAAGATGCGTTTTTCGTCATCGCTTATGCTGTCGTACTCATGTTCGAAAAACGTCATGATTGTAATATCGAGTTCGCGCATACCACAGCGGCATGc                                                                                         >  1:2094245/1‑99 (MQ=255)
      aaTAAAGATGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGcc                                                                                        <  2:324416/100‑1 (MQ=255)
       aTAAAGATGCGTTTTTCGTCATCGCTTAAGCTGTCGTCCTCATGTTCGAAAACCGGCAGGTGTGAAATATCGAGTTCGCGCATACCACAGCGGCATGccc                                                                                       <  2:2205095/100‑1 (MQ=255)
       aTAAAGATGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACGGCGGCATGccc                                                                                       <  2:1266157/100‑1 (MQ=255)
       aTAAAGATGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGccc                                                                                       >  2:67498/1‑100 (MQ=255)
          aaGATGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTGATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCaat                                                                                    >  2:1370767/1‑100 (MQ=255)
          aaGATGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCaat                                                                                    <  1:652440/100‑1 (MQ=255)
          aaGATGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATCATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCACaa                                                                                     >  2:1642134/1‑99 (MQ=255)
           aGATGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCaatg                                                                                   <  1:1252651/100‑1 (MQ=255)
            gATGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCaatga                                                                                  <  1:827495/100‑1 (MQ=255)
              tGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCaatgaat                                                                                >  1:2509237/1‑100 (MQ=255)
               gCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCaatgaatg                                                                               <  2:1814242/100‑1 (MQ=255)
                  tttttCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCgg                                                                            <  2:1916635/100‑1 (MQ=255)
                   ttttCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCggg                                                                           >  1:78040/1‑100 (MQ=255)
                      tCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCtt                                                                        <  1:2462501/100‑1 (MQ=255)
                       cGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCttt                                                                       <  1:2196016/100‑1 (MQ=255)
                       cGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCttt                                                                       <  2:2489936/100‑1 (MQ=255)
                        gTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTttg                                                                      >  1:496778/1‑100 (MQ=255)
                         tCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTttgt                                                                     <  1:1490522/100‑1 (MQ=255)
                          cATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTttgtt                                                                    >  1:480926/1‑100 (MQ=255)
                          cATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTttgtt                                                                    >  2:1022502/1‑100 (MQ=255)
                          cATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGAGCATACCACAGCGGCATGCCCAATGAATGCGGGCTttgtt                                                                    >  1:521029/1‑100 (MQ=255)
                           aTCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTttgttg                                                                   <  2:333784/100‑1 (MQ=255)
                           aTCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTttgttg                                                                   <  2:89718/100‑1 (MQ=255)
                           aTCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATGGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTttgttg                                                                   <  2:1802617/100‑1 (MQ=255)
                            tCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTtgttgt                                                                  >  1:1552734/1‑100 (MQ=255)
                             cGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTgttgtt                                                                 <  2:1668016/100‑1 (MQ=255)
                              gCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTa                                                                >  1:208482/1‑100 (MQ=255)
                               cTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTaa                                                               <  2:2012422/100‑1 (MQ=255)
                               cTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTaa                                                               <  2:592444/100‑1 (MQ=255)
                                 tAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATg                                                             <  2:2294713/100‑1 (MQ=255)
                                 tAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATg                                                             <  2:1927980/100‑1 (MQ=255)
                                 tAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATg                                                             <  2:480921/100‑1 (MQ=255)
                                  aaGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCACAATGAATGCGGGCTTTGTTGTTAATg                                                             <  2:82750/99‑1 (MQ=255)
                                     cTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCa                                                         <  2:668090/100‑1 (MQ=255)
                                     cTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCa                                                         <  2:1848979/100‑1 (MQ=255)
                                       gTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCAtc                                                       >  1:392208/1‑100 (MQ=255)
                                           tACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCAtcttct                                                   <  1:135679/100‑1 (MQ=255)
                                            aCTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATcttctt                                                  <  2:878127/100‑1 (MQ=255)
                                             cTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATcttctt                                                  <  2:2108095/99‑1 (MQ=255)
                                             cTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCttcttc                                                 >  1:1181623/1‑100 (MQ=255)
                                             cTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCttcttc                                                 >  2:1116860/1‑100 (MQ=255)
                                              tCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATTCCCAATGAATGCGGGCTTTGTTGTTAATCTCCATCTTCTTcc                                                <  1:1735185/100‑1 (MQ=255)
                                              tCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTcc                                                <  2:1275877/100‑1 (MQ=255)
                                              tCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTcc                                                <  1:51569/100‑1 (MQ=255)
                                              tCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTcc                                                <  1:95679/100‑1 (MQ=255)
                                               cATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCt                                               >  2:1300749/1‑100 (MQ=255)
                                               cATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCt                                               >  2:1201856/1‑100 (MQ=255)
                                                  gTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGtc                                            >  1:2485148/1‑100 (MQ=255)
                                                   ttCGAACAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGtct                                           <  2:1775927/100‑1 (MQ=255)
                                                     cGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCa                                         >  1:183969/1‑100 (MQ=255)
                                                     cGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCa                                         <  2:563263/100‑1 (MQ=255)
                                                      gAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCAc                                        <  2:1433072/100‑1 (MQ=255)
                                                         aaaCGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGaa                                     <  2:877208/100‑1 (MQ=255)
                                                         aaaCGGCATGATTGAAATATCGAGTGCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGa                                      <  2:838769/99‑1 (MQ=255)
                                                             ggCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATc                                 <  1:1766051/100‑1 (MQ=255)
                                                              gCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATcc                                >  2:1276098/1‑100 (MQ=255)
                                                              gCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATcc                                <  2:110976/100‑1 (MQ=255)
                                                              gCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATcc                                <  1:844846/100‑1 (MQ=255)
                                                              gCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATcc                                <  1:1038331/100‑1 (MQ=255)
                                                              gCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATcc                                >  2:471686/1‑100 (MQ=255)
                                                                 tGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATCCa                               >  1:1178460/1‑98 (MQ=255)
                                                                      gAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATCCAGTACCCg                        >  1:1007036/1‑100 (MQ=255)
                                                                      gAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATCCAGTACCCg                        >  1:546029/1‑100 (MQ=255)
                                                                        aaTATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATCCAGTACCCGGc                      <  1:768920/100‑1 (MQ=255)
                                                                        aaTATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATCCAGTACCCGGc                      <  1:1082719/100‑1 (MQ=255)
                                                                         atatCGAGTTCGCGCATACCACGGCGGCATGCCCAATGAATGCGGGCTTTGTTATTAATGTCCATCTTCTTCCTGTCTCACGAAAATCCAGTACCCCGCt                     <  1:420752/100‑1 (MQ=255)
                                                                          tatCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATCCAgatcggaagag                    >  1:98247/1‑90 (MQ=255)
                                                                          tatCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATCCAGTACCCGGCTa                    <  2:2237724/100‑1 (MQ=255)
                                                                           atCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTCATGTCCATCTTCTTCCTGTCTCACGAAAATCCAGTACCCGGCTAt                   >  2:1489804/1‑100 (MQ=255)
                                                                           atCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCCCGAAAATCCAGTACCCGGCTAt                   >  2:2494148/1‑100 (MQ=255)
                                                                              gAGTTCGCGCATACTACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATCCAGTACCCGGCTATTg                 <  2:2004828/99‑1 (MQ=255)
                                                                              gAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATCCAGTACCCGGCTATTGt                <  1:1490777/100‑1 (MQ=255)
                                                                                  tCGCGCATACCACAGCGGCGTGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATCCAGTACCCGGCTATTGTAACg            <  1:1470167/100‑1 (MQ=255)
                                                                                   cgcgCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATCCAGTACCCGGCTATTGTAACgt           >  2:438041/1‑100 (MQ=255)
                                                                                     cgcATACCACAGCGGCATGCCCAATGAATGCGGGCGTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATCCAGTACCCGGCTATTGTAACgtgt         <  2:942600/100‑1 (MQ=255)
                                                                                           ccACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATCCAGTACCCGGCTATTGTAACGTGTTTTTCg   <  1:840521/100‑1 (MQ=255)
                                                                                           ccACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATCCAGTACCCGGCTATTGTAACGTGTTTTTCg   >  1:1520298/1‑100 (MQ=255)
                                                                                           ccACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATCCAGTACCCGGCTATTGTAACGTGTTTTTCg   <  1:1012113/100‑1 (MQ=255)
                                                                                            cacaGCGGCATGCCCAATGAATGCGGGCTTTTTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATCCAGTACCCGGCTATTGTAACGTGTTTTTCGc  <  1:2412811/100‑1 (MQ=255)
                                                                                               |                                                                                                
CAGACGAATAAAGATGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACGGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATCCAGTACCCGGCTATTGTAACGTGTTTTTCGC  >  NC_000913/3040933‑3041124

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: