Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,525,969 T→A I61F (ATC→TTC)  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,9690TA98.4% 451.1 / ‑2.0 124I61F (ATC→TTC) nudEadenosine nucleotide hydrolase; Ap3A/Ap2A/ADP‑ribose/NADH hydrolase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base A (61/61);  total (62/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 142 total aligned reads displayed.

TTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGATCAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCATAAACACGCCGCACGCCATTGCTGAACTCCAG  >  NC_000913/3525872‑3526065
                                                                                                 |                                                                                                
  aGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGatcaggaacagg                                                                                              >  1:1120507/1‑100 (MQ=255)
  aGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGatcaggaacagg                                                                                              <  1:1853515/100‑1 (MQ=255)
      gCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCggatcaagatcggaagag                                                                                          >  1:2413454/1‑93 (MQ=255)
       cTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATc                                                                                         >  2:2123428/1‑100 (MQ=255)
       cTTCGTAGACGCTTTCACCCGGATAAAATAAAACTTACGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGAt                                                                                          <  2:283417/99‑1 (MQ=255)
          cGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTc                                                                                      >  1:2167258/1‑100 (MQ=255)
           gTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTcc                                                                                     <  1:2516354/100‑1 (MQ=255)
           gTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTcc                                                                                     <  1:2251519/100‑1 (MQ=255)
            tAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTcta                                                                                    <  1:2010985/100‑3 (MQ=255)
            tAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCca                                                                                    >  1:1019360/1‑100 (MQ=255)
             aGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCcac                                                                                   <  2:130828/100‑1 (MQ=255)
             aGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCcac                                                                                   <  1:2502840/100‑1 (MQ=255)
             aGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCcac                                                                                   <  2:55468/100‑1 (MQ=255)
                  cTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCg                                                                              >  2:426682/1‑100 (MQ=255)
                  cTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCg                                                                              <  2:665604/100‑1 (MQ=255)
                   tttCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCgg                                                                             <  2:856470/100‑1 (MQ=255)
                     tCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGca                                                                           <  1:491910/100‑1 (MQ=255)
                     tCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGca                                                                           <  1:746722/100‑1 (MQ=255)
                     tCACCCGGATCAATTAATCCTTTCGAAAAACCCAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGca                                                                           <  1:757596/100‑1 (MQ=255)
                      cACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGcac                                                                          <  2:443690/100‑1 (MQ=255)
                      cACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGcac                                                                          >  1:443692/1‑100 (MQ=255)
                       aCCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATCCGCGGATCAGGAACAGGTGATCGTCCACAATCGGcaca                                                                         >  2:1021285/1‑100 (MQ=255)
                        cccGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGcaca                                                                         >  2:736893/1‑99 (MQ=255)
                        cccGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACaa                                                                        >  1:894653/1‑100 (MQ=255)
                         ccGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAat                                                                       <  1:2546981/100‑1 (MQ=255)
                         ccGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAat                                                                       >  2:2146828/1‑100 (MQ=255)
                          cGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAatc                                                                      >  2:517514/1‑100 (MQ=255)
                             aTCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAatcatc                                                                   <  2:2409593/100‑1 (MQ=255)
                              tCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAAtcatca                                                                  >  1:320257/1‑100 (MQ=255)
                               cAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCAc                                                                 >  1:341680/1‑100 (MQ=255)
                                aaTTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACt                                                                >  2:24552/1‑100 (MQ=255)
                                aaTTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACt                                                                <  2:2554622/100‑1 (MQ=255)
                                aaTTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACt                                                                >  2:2948277/1‑100 (MQ=255)
                                aaTTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACt                                                                >  1:1214002/1‑100 (MQ=255)
                                 aTTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTg                                                               <  1:14870/100‑1 (MQ=255)
                                      tCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTcc                                                          <  2:2593419/100‑1 (MQ=255)
                                       ccTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTccc                                                         <  1:1259157/100‑1 (MQ=255)
                                        ntttCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCg                                                        >  2:882375/2‑100 (MQ=255)
                                        cTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGCTCGTCCACAATCGGCACAATCATCACTGCTTCCCg                                                        >  2:1723581/1‑100 (MQ=255)
                                           tCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCggtt                                                     >  2:2179435/1‑100 (MQ=255)
                                            cGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCggttg                                                    >  2:1516369/1‑100 (MQ=255)
                                             gAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCggttgg                                                   >  2:1376313/1‑100 (MQ=255)
                                             gAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCggttgg                                                   <  1:857599/100‑1 (MQ=255)
                                             gAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCggttgg                                                   <  1:2933294/100‑1 (MQ=255)
                                               aaaaCCTAATTCGTATGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCggttggtt                                                 <  1:2938153/100‑1 (MQ=255)
                                                aaaCCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGgttggttg                                                <  1:2853657/100‑1 (MQ=255)
                                                aaaCCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGgttggttg                                                <  1:1736103/100‑1 (MQ=255)
                                                 aaCCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGttggttgg                                               <  1:343248/100‑1 (MQ=255)
                                                    cTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACg                                            >  1:180091/1‑100 (MQ=255)
                                                     tAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGc                                           <  2:686362/100‑1 (MQ=255)
                                                     tAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGc                                           >  2:544779/1‑100 (MQ=255)
                                                     tAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGc                                           >  1:664193/1‑100 (MQ=255)
                                                       aTTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCAt                                         <  2:2083187/100‑1 (MQ=255)
                                                        ttCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCAtt                                        >  2:1974151/1‑100 (MQ=255)
                                                           gTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGt                                     <  1:48285/100‑1 (MQ=255)
                                                           gTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGt                                     >  1:1344910/1‑100 (MQ=255)
                                                           gTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGt                                     <  2:320255/100‑1 (MQ=255)
                                                            tAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGtt                                    >  1:2443602/1‑100 (MQ=255)
                                                             aGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTc                                   <  1:309334/100‑1 (MQ=255)
                                                              ggATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATGATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCa                                  >  1:2099430/1‑100 (MQ=255)
                                                              ggATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCa                                  >  2:257548/1‑100 (MQ=255)
                                                               gATTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCAt                                 >  2:2830398/1‑100 (MQ=255)
                                                                aTTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCATa                                <  2:1019355/100‑1 (MQ=255)
                                                                aTTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCATa                                <  2:1533365/100‑1 (MQ=255)
                                                                aTTCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCATa                                <  1:2948293/100‑1 (MQ=255)
                                                                 ttCAGTTCCCACTGCGTCTTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCATaa                               <  1:1021290/100‑1 (MQ=255)
                                                                 ttCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCATaa                               >  1:1712376/1‑100 (MQ=255)
                                                                 ttCAGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATAGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCATaa                               >  1:2044850/1‑100 (MQ=255)
                                                                    aGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCATAAaca                            >  2:265513/1‑100 (MQ=255)
                                                                    aGTTCCCACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCATAAaca                            >  1:2636591/1‑100 (MQ=255)
                                                                        cccACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCATAAACAcgcc                        >  1:2852746/1‑100 (MQ=255)
                                                                          cACTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCATAAACAcgccgc                      <  1:2546306/100‑1 (MQ=255)
                                                                           aCTGCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCATAAACACGCCGCa                     <  2:253486/100‑1 (MQ=255)
                                                                              gCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCATAAACACGCCGCACGc                  >  1:2241957/1‑100 (MQ=255)
                                                                              gCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCATAAACACGCCGCACGc                  >  1:2415959/1‑100 (MQ=255)
                                                                              gCGTATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCATAAACACGCCGCACGc                  >  1:93438/1‑100 (MQ=255)
                                                                                 tATTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCATAAACACGCCGCACGCCAt               >  2:1979735/1‑100 (MQ=255)
                                                                                  aTTCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCATAAACACGCCGCACGCCAtt              >  1:1003890/1‑100 (MQ=255)
                                                                                   ttCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCATAAACACGCCGCACGCCATTg             >  2:662088/1‑100 (MQ=255)
                                                                                   ttCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCATAAACACGCCGCACGCCATTg             <  1:2693433/100‑1 (MQ=255)
                                                                                    tCGCGGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCATAAACACGCCGCACGCCATTGc            <  1:874103/100‑1 (MQ=255)
                                                                                     cgcgGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCATAAACACGCCGCACGCCATTGCt           >  2:1866821/1‑100 (MQ=255)
                                                                                     cgcgGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCATAAACACGCCGCACGCCATTGCt           >  1:2668316/1‑100 (MQ=255)
                                                                                     cgcgGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCATAAACACGCCGCACGCCATTGCt           >  2:2971424/1‑100 (MQ=255)
                                                                                      gcgGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCATAAACACGCCGCACGCCATTGCTg          <  1:213596/100‑1 (MQ=255)
                                                                                      gcgGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCATAAACACGCCGCACGCCATTGCTg          <  1:1399651/100‑1 (MQ=255)
                                                                                      gcgGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCATAAACACGCCGCACGCCATTGCTg          <  2:331722/100‑1 (MQ=255)
                                                                                      gcgGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCATAAACACGCCGCACGCCATTGCTg          <  2:130579/100‑1 (MQ=255)
                                                                                       cgGATCAGGAACAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCATAAACACGCCGCACGCCATTGCTGa         <  2:663088/100‑1 (MQ=255)
                                                                                         gatcaggaacagGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCATAAACACGCCGCACGCCATTGCTGAAc       <  1:629144/100‑1 (MQ=255)
                                                                                         gatcaggaacagGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCATAAACACGCCGCACGCCATTGCTGAAc       <  1:544782/100‑1 (MQ=255)
                                                                                          atcaggaacaggTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCATAAACACGCCGCACGCCATTGCTGAACt      >  1:463454/1‑100 (MQ=255)
                                                                                            caggaacaggGGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGCTCTTAAACACGCCGCACGCCATTGCTGAACTcc    <  2:1962063/100‑1 (MQ=255)
                                                                                             aggaacaggTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCATAAACACGCCGCACGCCATTGCTGAACTCCa   >  1:932452/1‑100 (MQ=255)
                                                                                             aggaacaggTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCATAAACACGCCGCACGCCATTGCTGAACTCCa   <  1:1376321/100‑1 (MQ=255)
                                                                                             aggaacaggTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCATAAACACGCCGCACGCCATTGCTGAACTCCa   <  1:1052606/100‑1 (MQ=255)
                                                                                             aggaacaggTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCATAAACACGCCGCACGCCATTGCTGAACTCCa   <  2:2415944/100‑1 (MQ=255)
                                                                                              ggaacaggTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCATAAACACGCCGCACGCCATTGCTGAACTCCAg  >  1:1592589/1‑100 (MQ=255)
                                                                                              ggaacaggTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCATAAACACGCCGCACGCCATTGCTGAACTCCAg  >  1:1486118/1‑100 (MQ=255)
                                                                                              ggaacaggTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCATAAACACGCCGCACGCCATTGCTGAACTCCAg  >  1:612785/1‑100 (MQ=255)
                                                                                                 |                                                                                                
TTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGATCAGGTGATCGTCCACAATCGGCACAATCATCACTGCTTCCCGGTTGGTTGGACGCATTCGTTCATAAACACGCCGCACGCCATTGCTGAACTCCAG  >  NC_000913/3525872‑3526065

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