Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,304,649 C→T 15.8% intergenic (+256/+459) eco → / ← mqo ecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,304,6490CT15.8% 194.6 / 23.3 87intergenic (+256/+459)eco/mqoecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base C (43/30);  new base T (10/4);  total (53/34)
Fisher's exact test for biased strand distribution p-value = 5.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01

ACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGC  >  NC_000913/2304550‑2304714
                                                                                                   |                                                                 
tctccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTc                                                                   <  1:285444/97‑1 (MQ=17)
 cgccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCg                                                                  <  1:814261/100‑1 (MQ=33)
 cgccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCg                                                                  <  2:1907974/100‑1 (MQ=33)
  gccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTgg                                                                 >  1:1102505/1‑100 (MQ=18)
  gccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCgg                                                                 <  1:1519778/100‑1 (MQ=33)
  gccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCgg                                                                 <  1:441319/100‑1 (MQ=32)
  gccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCgg                                                                 <  1:1660484/100‑1 (MQ=32)
  gccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCgg                                                                 <  2:2776164/100‑1 (MQ=33)
  gccgcATCCGACAGCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGCCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCgg                                                                 <  1:2074008/100‑1 (MQ=12)
   ccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCgg                                                                 >  2:1862013/1‑99 (MQ=33)
   ccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGa                                                                <  2:2899793/100‑1 (MQ=33)
   ccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTAGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCAGa                                                                >  2:285446/1‑100 (MQ=12)
    cgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAt                                                               <  1:327332/100‑1 (MQ=32)
    cgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAt                                                               <  2:733362/100‑1 (MQ=33)
     gcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTGCACCCCGTAGGTCGGAt                                                               <  1:2728224/99‑1 (MQ=17)
     gcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATa                                                              >  2:173910/1‑100 (MQ=33)
     gcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATa                                                              <  2:2343013/100‑1 (MQ=32)
      cATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATaa                                                             >  2:2335797/1‑100 (MQ=33)
        tCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAgg                                                           >  2:1194274/1‑100 (MQ=33)
         ccGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGc                                                          <  1:153774/100‑1 (MQ=14)
         ccGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGc                                                          <  2:367302/100‑1 (MQ=33)
          cGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGc                                                          <  1:2118152/99‑1 (MQ=14)
          cGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAgg                                                           <  1:761981/98‑1 (MQ=32)
          cGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCg                                                         >  2:2617518/1‑100 (MQ=33)
           gtcATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCgg                                                        <  1:1768263/98‑1 (MQ=17)
           gACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCgg                                                        <  1:2617507/100‑1 (MQ=17)
           gACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCgg                                                        <  2:1677809/100‑1 (MQ=33)
           gACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCgg                                                        <  2:709898/100‑1 (MQ=33)
           gACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCgg                                                        >  2:62897/1‑100 (MQ=33)
            aCATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGt                                                       <  1:1577269/100‑1 (MQ=33)
             cctCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGtt                                                      >  2:1957775/3‑100 (MQ=18)
             cATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGtt                                                      <  1:2515178/100‑1 (MQ=14)
             cATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGtt                                                      <  1:2903002/100‑1 (MQ=32)
             cATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGtt                                                      <  1:1077659/100‑1 (MQ=32)
              aTCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTa                                                     <  1:2280221/100‑1 (MQ=32)
               tCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTCGGTCGGATAAGGCGGTTAc                                                    >  1:719990/1‑100 (MQ=18)
               tCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAc                                                    >  2:1873085/1‑100 (MQ=33)
                ccAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGTTTAcg                                                   >  1:543312/1‑100 (MQ=18)
                 cAACGCCCGGGCCGGTTGCCTGATGGGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAcgc                                                  <  2:209348/100‑1 (MQ=12)
                  aaCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAcgcc                                                 >  2:2581918/1‑100 (MQ=32)
                  aaCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAcgcc                                                 >  2:1410449/1‑100 (MQ=32)
                  aaCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAcgcc                                                 <  1:1879289/100‑1 (MQ=33)
                   aCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAcgccg                                                <  2:2510474/100‑1 (MQ=33)
                   aCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAcgccg                                                <  2:1661892/100‑1 (MQ=33)
                   aCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCCCCACCCCGTAGGTCGGATAAGGCGGTTAcgccg                                                <  1:749432/100‑1 (MQ=18)
                    cGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAcgccg                                                <  1:1080365/99‑1 (MQ=32)
                     gCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCa                                              >  2:581164/1‑100 (MQ=31)
                      cccGAGCCGGTTGCCTGATGAGACGCTGGCGCGTCTTATCAGCACTCCACAGCTCTGAAGTGCTCCACCCCGTATGTCGTATAAGGCGGCTAcgc                                                  >  1:1707760/1‑95 (MQ=11)
                            gtgtaggcctgataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGa                                         >  1:2389907‑M2/32‑98 (MQ=255)
                               taggcctgataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATc                                     >  2:1792941‑M2/29‑99 (MQ=255)
                                aggcctgataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTa                                   >  2:2403336‑M2/28‑100 (MQ=255)
                                aggcctgataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTa                                   >  1:1684657‑M2/28‑100 (MQ=255)
                                aggcctgataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGCCGGATAAGGCGGTTACGCCGCATCCGACATCTa                                   >  2:1152171‑M2/28‑100 (MQ=255)
                                aggcctgataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACatcca                                   >  1:1239331‑M2/28‑98 (MQ=255)
                                 ggcctgataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGCTAAAGAGGTTACGCCGCAGCCGACATCCaa                                  >  1:771617‑M2/27‑100 (MQ=255)
                                 ggcctgataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTaa                                  >  1:2057169‑M2/27‑100 (MQ=255)
                                 ggcctgataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTaa                                  >  1:2597534‑M2/27‑100 (MQ=255)
                                 ggcctgataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTaa                                  >  2:2037153‑M2/27‑100 (MQ=255)
                                 ggcctgataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTa                                   >  2:2158202‑M2/27‑99 (MQ=255)
                                 ggcctgataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTACTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTaa                                  >  1:745855‑M2/27‑100 (MQ=255)
                                  gcctgataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAAc                                 >  2:1904760‑M2/26‑100 (MQ=255)
                                  gcctgataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAAc                                 >  1:585965‑M2/26‑100 (MQ=255)
                                  gcctgataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAAc                                 >  1:35803‑M2/26‑100 (MQ=255)
                                  gcctgataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAAc                                 >  2:1430948‑M2/26‑100 (MQ=255)
                                  gcctgataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAAc                                 >  2:2396334‑M2/26‑100 (MQ=255)
                                  gcctgataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAAc                                 >  2:705454‑M2/26‑100 (MQ=255)
                                  gcctgataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAAc                                 >  1:2581567‑M2/26‑100 (MQ=255)
                                   cctgataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACg                                >  1:2312960‑M2/25‑100 (MQ=255)
                                   cctgataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAAc                                 >  2:1463371‑M2/25‑99 (MQ=255)
                                   cctgataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACg                                >  1:425803‑M2/25‑100 (MQ=255)
                                   cctgataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACg                                >  1:1807162‑M2/25‑100 (MQ=255)
                                      gataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGccc                             >  2:1670050‑M2/22‑100 (MQ=255)
                                         aagacgcgccagcgtcttATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAg                          >  2:1525038‑M2/19‑100 (MQ=255)
                                               cgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGtt                    >  2:2836010‑M2/13‑100 (MQ=255)
                                                          cATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGa         >  1:354828‑M2/2‑100 (MQ=255)
                                                            gCAGTCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCGGATAAGGCGGTCACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGAc        <  1:2598962‑M2/98‑1 (MQ=255)
                                                            TCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACg       >  2:2537166‑M2/1‑100 (MQ=255)
                                                            TCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACg       <  2:440057‑M2/100‑1 (MQ=255)
                                                            TCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCAAGCCGGTTGCCTGATGCGACg       >  1:2394810‑M2/1‑100 (MQ=255)
                                                             CAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGc      >  1:472084‑M2/1‑100 (MQ=255)
                                                             CAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGc      >  1:2223682‑M2/1‑100 (MQ=255)
                                                              AGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGATGCCTGATGCGACg       >  1:1895729‑M2/1‑98 (MQ=255)
                                                               GGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTg    >  2:1498103‑M2/1‑100 (MQ=255)
                                                               GGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTg    >  2:1497097‑M2/1‑100 (MQ=255)
                                                               GGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTg    >  2:1925356‑M2/1‑100 (MQ=255)
                                                               GGCCTACACCGCAGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTg    >  2:1913302‑M2/1‑100 (MQ=255)
                                                                GCCTAACCCGCCGTGCAGTGCTTCTCCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTgg   <  1:726311‑M2/100‑1 (MQ=255)
                                                                 CCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGc  >  1:1897700‑M2/1‑100 (MQ=255)
                                                                                                   |                                                                 
ACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGC  >  NC_000913/2304550‑2304714

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: