Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,397,925 Δ1 bp coding (454/594 nt) yhdE ← dTTP/UTP pyrophosphatase; m(5)UTP/m(5)CTP/pseudo‑UTP pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,397,9250T.99.2% 540.6 / ‑2.5 120coding (454/594 nt)yhdEdTTP/UTP pyrophosphatase; m(5)UTP/m(5)CTP/pseudo‑UTP pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (55/64);  total (56/64)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 120 total aligned reads displayed.

CTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACC  >  NC_000913/3397831‑3398019
                                                                                              |                                                                                              
cTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTa                                                                                          >  1:1924211/1‑100 (MQ=255)
cTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTa                                                                                          <  2:1869489/100‑1 (MQ=255)
cTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTa                                                                                          <  1:1846223/100‑1 (MQ=255)
cTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTa                                                                                          >  2:3027627/1‑100 (MQ=255)
cTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTa                                                                                          <  2:3124622/100‑1 (MQ=255)
 ttAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTAt                                                                                         <  1:487179/100‑1 (MQ=255)
  taataaTTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATg                                                                                        >  2:1609696/1‑100 (MQ=255)
  taataaTTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATg                                                                                        <  1:1303595/100‑1 (MQ=255)
   aataatTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGc                                                                                       >  1:237616/1‑100 (MQ=255)
    ataatTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCa                                                                                      >  2:2174351/1‑100 (MQ=255)
      aatTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCAcc                                                                                    >  2:1034892/1‑100 (MQ=255)
      aatTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCAcc                                                                                    <  1:2770514/100‑1 (MQ=255)
       atTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCt                                                                                   >  1:2902179/1‑100 (MQ=255)
        ttCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTg                                                                                  <  1:110530/100‑1 (MQ=255)
        ttCATACGTTTCAACCAGCGGTAAGCCGACTACCGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTg                                                                                  <  2:1614803/100‑1 (MQ=255)
         tCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTg                                                                                  >  2:1592480/1‑99 (MQ=255)
         tCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGc                                                                                 <  1:1089765/100‑1 (MQ=255)
          cATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCt                                                                                >  2:3012497/1‑100 (MQ=255)
           aTACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCtt                                                                               >  2:3157813/1‑100 (MQ=255)
            tACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCttt                                                                              >  1:1269581/1‑100 (MQ=255)
             aCGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTa                                                                             <  1:1336887/100‑1 (MQ=255)
             aCGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATCCACTTGCTTTa                                                                             >  2:1308422/1‑100 (MQ=255)
              cGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTa                                                                             <  2:2531076/99‑1 (MQ=255)
               gTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCGGCCCCTGAAACCGTATGCACCTGCTTTATc                                                                           <  2:2030712/100‑1 (MQ=255)
               gTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATc                                                                           >  1:2330939/1‑100 (MQ=255)
               gTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATc                                                                           <  2:2960757/100‑1 (MQ=255)
                   cAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGTCCCTGAAACCGTATGCACCTGCTTTATCTAAc                                                                       >  2:186654/1‑100 (MQ=255)
                    aaCCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACg                                                                      >  2:60884/1‑100 (MQ=255)
                     aCCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACgg                                                                     >  1:156580/1‑100 (MQ=255)
                     aCCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACgg                                                                     <  1:2547540/100‑1 (MQ=255)
                     aCCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACgg                                                                     >  1:101743/1‑100 (MQ=255)
                       cAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGtt                                                                   <  1:3171284/100‑1 (MQ=255)
                          cGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCa                                                                 <  1:3039878/99‑1 (MQ=255)
                           ggTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATc                                                               >  1:2068970/1‑100 (MQ=255)
                              aaGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCt                                                            <  1:772681/100‑1 (MQ=255)
                                gCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTgg                                                          >  1:738942/1‑100 (MQ=255)
                                gCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTgg                                                          >  2:758780/1‑100 (MQ=255)
                                gCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTgg                                                          >  1:932816/1‑100 (MQ=255)
                                gCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTg                                                           >  2:1804635/1‑99 (MQ=255)
                                  cGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCg                                                        <  1:399797/100‑1 (MQ=255)
                                   gACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGa                                                       <  2:237610/100‑1 (MQ=255)
                                       aCGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATa                                                   >  1:2086519/1‑100 (MQ=255)
                                        cGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAg                                                  >  1:2852215/1‑100 (MQ=255)
                                        cGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAg                                                  <  1:2832167/100‑1 (MQ=255)
                                         ggCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGc                                                 <  1:998309/100‑1 (MQ=255)
                                          gCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCGGGctatgtcgcc                                                <  2:232041/100‑10 (MQ=255)
                                          gCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCGGGctatgtcgcc                                                <  2:233940/100‑10 (MQ=255)
                                           cGTGATAGCTGCCATTTATCTTCCTGACAATACTGTCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAg                                                  <  2:2760590/97‑1 (MQ=255)
                                           cGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGccc                                               >  2:1605509/1‑100 (MQ=255)
                                            gTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCcg                                              >  1:2851765/1‑100 (MQ=255)
                                             tGATAACTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCACTGGCGACATAGCCcgc                                             >  2:3073406/1‑100 (MQ=255)
                                                tAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGAt                                          <  1:2911881/100‑1 (MQ=255)
                                                 aGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATg                                         >  2:704715/1‑100 (MQ=255)
                                                  gCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGt                                        >  2:2662077/1‑100 (MQ=255)
                                                  gCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGt                                        >  1:1646616/1‑100 (MQ=255)
                                                   cTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTc                                       >  2:453203/1‑100 (MQ=255)
                                                   cTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTc                                       <  2:722686/100‑1 (MQ=255)
                                                     gCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCtt                                     <  1:2988649/100‑1 (MQ=255)
                                                      ccATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTc                                    <  1:3012533/100‑1 (MQ=255)
                                                          ttATCTTCCTGACAAAACAGCCACCCAGTCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCt                                <  1:186659/100‑1 (MQ=255)
                                                          ttATCTTCCTGACAAAACAGCCACCCAACCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTc                                 <  2:156576/99‑1 (MQ=255)
                                                           tATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTg                               <  1:780120/100‑1 (MQ=255)
                                                            aTCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGt                              >  2:208495/1‑100 (MQ=255)
                                                             tCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGtt                             <  2:599456/100‑1 (MQ=255)
                                                             tCTTCCTGACAAAACAACCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGtt                             >  1:103363/1‑100 (MQ=255)
                                                              cTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTa                            <  1:985385/100‑1 (MQ=255)
                                                               ttCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTaa                           >  2:2929410/1‑100 (MQ=255)
                                                               ttCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTaa                           <  1:822098/100‑1 (MQ=255)
                                                                tCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAAc                          >  1:2189582/1‑100 (MQ=255)
                                                                 ccTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACg                         <  2:1269567/100‑1 (MQ=255)
                                                                  cTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGt                        >  1:1551713/1‑100 (MQ=255)
                                                                   tGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGtt                       >  2:338186/1‑100 (MQ=255)
                                                                   tGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGtt                       >  2:1156471/1‑100 (MQ=255)
                                                                    gACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTc                      <  1:1171996/100‑1 (MQ=255)
                                                                    gACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTc                      >  2:2612716/1‑100 (MQ=255)
                                                                      cAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTg                    >  1:1530652/1‑100 (MQ=255)
                                                                       aaaaCAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCt                     <  2:1163338/98‑1 (MQ=255)
                                                                        aaaCAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGATATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGaa                  <  1:758672/100‑1 (MQ=255)
                                                                          aCAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGaaa                 <  2:2750612/99‑1 (MQ=255)
                                                                          aCAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAg                <  2:2003360/100‑1 (MQ=255)
                                                                          aCAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCAGCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAg                <  2:1951179/100‑1 (MQ=255)
                                                                            aGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTc              <  1:1159984/100‑1 (MQ=255)
                                                                             gCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTca             <  1:1571861/100‑1 (MQ=255)
                                                                             gCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTca             <  2:425110/100‑1 (MQ=255)
                                                                              ccACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTcac            >  1:581252/1‑100 (MQ=255)
                                                                               cACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTcaca           <  1:2594362/100‑1 (MQ=255)
                                                                                aCCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACAt          <  1:2174377/100‑1 (MQ=255)
                                                                                 nccAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATc         >  1:325493/2‑100 (MQ=255)
                                                                                  ccAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCg        <  2:249803/100‑1 (MQ=255)
                                                                                   cAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCgg       >  1:1804762/1‑100 (MQ=255)
                                                                                   cAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCgg       <  2:1807643/100‑1 (MQ=255)
                                                                                    aGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGt      <  1:1034904/100‑1 (MQ=255)
                                                                                     gCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTg     <  2:101740/100‑1 (MQ=255)
                                                                                     gCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTg     <  2:1285974/100‑1 (MQ=255)
                                                                                      ccccTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTg     >  1:2837045/1‑100 (MQ=255)
                                                                                      ccccTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGa    <  2:2871313/100‑1 (MQ=255)
                                                                                      ccccTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGa    <  1:1864185/100‑1 (MQ=255)
                                                                                       cccTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGac   <  1:1403223/100‑1 (MQ=255)
                                                                                       cccTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGac   >  1:243159/1‑100 (MQ=255)
                                                                                        ccTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGacc  >  2:2063865/1‑100 (MQ=255)
                                                                                              |                                                                                              
CTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACC  >  NC_000913/3397831‑3398019

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