Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,525,601 Δ1 bp 100% coding (549/561 nt) nudE ← adenosine nucleotide hydrolase; Ap3A/Ap2A/ADP‑ribose/NADH hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,525,5980C.79.7% 341.4 / 113.8 128coding (552/561 nt)nudEadenosine nucleotide hydrolase; Ap3A/Ap2A/ADP‑ribose/NADH hydrolase
Reads supporting (aligned to +/- strand):  ref base C (14/12);  new base . (45/57);  total (59/69)
Fisher's exact test for biased strand distribution p-value = 3.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 138 total aligned reads displayed.

TATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCCATCATATGCG  >  NC_000913/3525500‑3525692
                                                                                                  |                                                                                              
tatccggcTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGcc                                                                                               >  1:311368‑M2/9‑100 (MQ=255)
 atccggcTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGccc                                                                                              <  2:1487173‑M2/93‑1 (MQ=255)
 atccggcTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGccc                                                                                              <  2:262788‑M2/93‑1 (MQ=255)
   ccggcTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTcgccctg                                                                                            >  2:2007135‑M2/6‑98 (MQ=255)
    cggcTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACtcgccctgc                                                                                           >  1:2828597‑M2/5‑97 (MQ=255)
    cggcTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACtcgccctgc                                                                                           >  1:1407082‑M2/5‑97 (MQ=255)
     ggcTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACActcgccctgcc                                                                                          >  2:2636554‑M2/4‑96 (MQ=255)
        tAACAAAATCGTGAAAACTCAATAAATTGCAGAAACAACATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCtt                                                                                      >  1:1900378/1‑100 (MQ=255)
          aCAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAATCATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTc                                                                                    <  2:2428673/100‑1 (MQ=255)
          aCAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCttt                                                                                     >  2:1278201/1‑100 (MQ=255)
          aCAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTc                                                                                    <  1:1659535/100‑1 (MQ=255)
          aCAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTc                                                                                    <  2:2291087/100‑1 (MQ=255)
           cAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTc                                                                                    >  2:120618/1‑100 (MQ=255)
           cAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCa                                                                                   <  1:1504816/100‑1 (MQ=255)
           cAAAATCGTGAAAACTCAATAAATTGCAGAAACACCATAAGCCTGATAAACATTGTACATCAGGCAAACTTCACCCATTTACAGTCGCCCCTGCCCTTTc                                                                                    >  2:2792329/1‑100 (MQ=255)
            aaaaTCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCaa                                                                                  <  2:1948680/100‑1 (MQ=255)
             aaaTCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAAc                                                                                 <  1:2396555/100‑1 (MQ=255)
             aaaTCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAAc                                                                                 <  1:587857/100‑1 (MQ=255)
               aTCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCa                                                                               <  2:2838891/100‑1 (MQ=255)
                 cGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCAtt                                                                             >  2:2077541/1‑100 (MQ=255)
                  gTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTc                                                                            >  2:1027342/1‑100 (MQ=255)
                   tGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTc                                                                            >  1:1177046/1‑100 (MQ=255)
                   tGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTc                                                                            >  1:1222944/1‑100 (MQ=255)
                    gAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTcgc                                                                          <  1:1433841/100‑1 (MQ=255)
                    gAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTcgc                                                                          <  2:1438639/100‑1 (MQ=255)
                     aaaaCTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTcgcg                                                                         >  2:157768/1‑100 (MQ=255)
                        aCTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCa                                                                       >  2:1585093/1‑100 (MQ=255)
                           cAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAg                                                                   >  2:186505/1‑100 (MQ=255)
                           cAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAg                                                                   <  1:1100067/100‑1 (MQ=255)
                            aaTAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGa                                                                    >  1:654361/1‑99 (MQ=255)
                            aaTAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAgg                                                                  >  2:400876/1‑100 (MQ=255)
                             aTAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAgg                                                                  >  2:559806/1‑100 (MQ=255)
                             aTAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAgg                                                                  >  2:515721/1‑100 (MQ=255)
                              tAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGaa                                                                >  1:418501/1‑100 (MQ=255)
                                 aTTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATATACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACa                                                              >  1:1864860/1‑100 (MQ=255)
                                   tGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAgcg                                                           <  1:2381425/100‑1 (MQ=255)
                                   tGCACAAACCCCATAAGCATGATAAACAGTGCGCAGCAGCCAAACTTCACGCATTCACACTCGCCCCGCCCTTTCCACCAGTCGCGCACGATAAACAgc                                                            <  1:1813299/99‑1 (MQ=255)
                                    gCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAgcg                                                           >  1:369675/1‑100 (MQ=255)
                                    gCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTAAAATCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAgcg                                                           >  1:2626363/1‑99 (MQ=255)
                                     cAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCa                                                         <  1:515381/100‑1 (MQ=255)
                                       gAAACCCCATAAGCCTGATAAACATTGTGCATCAGACAAACTTCACGCATTTACAATCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCAc                                                        >  1:1004726/1‑99 (MQ=255)
                                            cccATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACAt                                                  >  1:231276/1‑100 (MQ=255)
                                                 aaGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTgcgc                                             <  1:1027340/100‑1 (MQ=255)
                                                    ccTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCt                                          <  1:2227283/100‑1 (MQ=255)
                                                    ccTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCt                                          <  2:1944856/100‑1 (MQ=255)
                                                    ccTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCt                                          >  1:1527495/1‑100 (MQ=255)
                                                     cTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCt                                          <  1:515720/100‑1 (MQ=255)
                                                     cTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCtt                                         >  2:1898197/1‑100 (MQ=255)
                                                      tGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCtt                                         >  1:1324570/1‑100 (MQ=255)
                                                       gATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCa                                       <  1:2049959/100‑1 (MQ=255)
                                                        aTAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCAt                                      <  1:73623/100‑1 (MQ=255)
                                                         tAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCAtt                                     >  2:285608/1‑100 (MQ=255)
                                                         tAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCAtt                                     <  2:11272/100‑1 (MQ=255)
                                                         tAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCAtt                                     >  2:239831/1‑100 (MQ=255)
                                                          aaaCATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCGTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTg                                    <  1:1000868/100‑1 (MQ=255)
                                                           aaCATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTg                                    <  2:78124/100‑1 (MQ=255)
                                                           aaCATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGa                                   >  1:1325432/1‑100 (MQ=255)
                                                           aaCATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGa                                   <  1:2823175/100‑1 (MQ=255)
                                                           aaCATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGa                                   <  2:915058/100‑1 (MQ=255)
                                                            aCATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGaa                                  >  1:2513988/1‑100 (MQ=255)
                                                            aCATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGaa                                  >  2:1925301/1‑100 (MQ=255)
                                                            aCATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGaa                                  >  2:1433669/1‑100 (MQ=255)
                                                             cATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAg                                 <  2:318626/100‑1 (MQ=255)
                                                             cATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAg                                 <  2:889465/100‑1 (MQ=255)
                                                              aTTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAg                                 <  1:1620458/100‑1 (MQ=255)
                                                                    cATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAggg                           <  2:689328/100‑1 (MQ=255)
                                                                    cATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGCAGACAggg                           >  1:69446/1‑99 (MQ=255)
                                                                    cATCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGt                          >  1:34370/1‑100 (MQ=255)
                                                                    cATCAGGCAAACTTCACGCATTTACACGCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGt                          <  2:391103/100‑1 (MQ=255)
                                                                     aTCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGt                          <  1:2805096/100‑1 (MQ=255)
                                                                     aTCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAggg                           <  1:714980/98‑1 (MQ=255)
                                                                     aTCAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAggg                           >  1:748413/1‑98 (MQ=255)
                                                                       cAGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCtt                       >  1:1939510/1‑100 (MQ=255)
                                                                        aGGCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTc                      >  2:402919/1‑100 (MQ=255)
                                                                         ggCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCCTTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTcc                     <  1:1234546/100‑1 (MQ=255)
                                                                         ggCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTcc                     >  2:1286902/1‑100 (MQ=255)
                                                                         ggCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTcc                     <  1:1286900/100‑1 (MQ=255)
                                                                         ggCAAACTTCACGCATTTACACTCGCCCTGCCCTTACAACCATTCGCGCACGAGGAACAGCGCACTGACATTGAGAGCTTCATTGAAGTCAGGGTCTTcc                     >  2:1733599/1‑100 (MQ=255)
                                                                          gCAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCa                    >  1:898241/1‑100 (MQ=255)
                                                                           cAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAg                   <  2:2529214/100‑1 (MQ=255)
                                                                           cAAACTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAg                   <  2:2265596/100‑1 (MQ=255)
                                                                            aaaCTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAg                   <  1:2737067/100‑1 (MQ=255)
                                                                             aaCTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCa                 <  2:1135635/100‑1 (MQ=255)
                                                                             aaCTTCACGCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCa                 <  1:740815/100‑1 (MQ=255)
                                                                              aCTTCACGCATTTACAGTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCaa                <  1:1701245/100‑1 (MQ=255)
                                                                                 tCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAAt              <  2:460350/100‑1 (MQ=255)
                                                                                    cGCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCcat          >  1:2409818/1‑100 (MQ=255)
                                                                                    cGCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCcat          <  1:586905/100‑1 (MQ=255)
                                                                                    cGCATTTACACTCGCCCTGCACTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCcat          >  1:1778860/1‑100 (MQ=255)
                                                                                     gCATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCcatn         <  1:859674/100‑2 (MQ=255)
                                                                                      cATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCcatc         <  1:453708/100‑1 (MQ=255)
                                                                                      cATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCcatca        <  1:1093753/100‑1 (MQ=255)
                                                                                      cATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCcatca        <  1:2315351/100‑1 (MQ=255)
                                                                                      cATTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCcatca        >  1:2234581/1‑100 (MQ=255)
                                                                                       aTTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCcatca        <  2:762431/100‑1 (MQ=255)
                                                                                       aTTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCcatcat       <  1:2778804/100‑1 (MQ=255)
                                                                                       aTTTACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGAGTCTTCCAGCAAATCcatcat       <  1:2509997/100‑1 (MQ=255)
                                                                                        tttACACTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCCATCata      <  1:1413831/100‑1 (MQ=255)
                                                                                           acacTCGCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCCATCATATgc   >  2:2400788/1‑100 (MQ=255)
                                                                                             acTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCCATCATATgcg  <  1:341765/100‑1 (MQ=255)
                                                                                                  |                                                                                              
TATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCATTTACACTCGCCCCTGCCCTTTCAACCATTCGCGCACGAGGAACAGCGCACTGACATTGCGCGCTTCATTGAAGTCAGGGTCTTCCAGCAAATCCATCATATGCG  >  NC_000913/3525500‑3525692

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: