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
                                                                                              |                                                                                              
cTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTa                                                                                          >  1:1924211/1‑100 (MQ=255)
cTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTa                                                                                          <  2:1869489/100‑1 (MQ=255)
cTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTa                                                                                          <  1:1846223/100‑1 (MQ=255)
cTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTa                                                                                          >  2:3027627/1‑100 (MQ=255)
cTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTa                                                                                          <  2:3124622/100‑1 (MQ=255)
 ttAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTAt                                                                                         <  1:487179/100‑1 (MQ=255)
  taataaTTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATg                                                                                        >  2:1609696/1‑100 (MQ=255)
  taataaTTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATg                                                                                        <  1:1303595/100‑1 (MQ=255)
   aataatTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGc                                                                                       >  1:237616/1‑100 (MQ=255)
    ataatTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCa                                                                                      >  2:2174351/1‑100 (MQ=255)
      aatTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCAcc                                                                                    >  2:1034892/1‑100 (MQ=255)
      aatTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCAcc                                                                                    <  1:2770514/100‑1 (MQ=255)
       atTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCt                                                                                   >  1:2902179/1‑100 (MQ=255)
        ttCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTg                                                                                  <  1:110530/100‑1 (MQ=255)
        ttCATACGTTTCAACCAGCGGTAAGCCGACTACCGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTg                                                                                  <  2:1614803/100‑1 (MQ=255)
         tCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTg                                                                                  >  2:1592480/1‑99 (MQ=255)
         tCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGc                                                                                 <  1:1089765/100‑1 (MQ=255)
          cATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCt                                                                                >  2:3012497/1‑100 (MQ=255)
           aTACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCtt                                                                               >  2:3157813/1‑100 (MQ=255)
            tACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCttt                                                                              >  1:1269581/1‑100 (MQ=255)
             aCGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTa                                                                             <  1:1336887/100‑1 (MQ=255)
             aCGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATCCACTTGCTTTa                                                                             >  2:1308422/1‑100 (MQ=255)
              cGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTa                                                                             <  2:2531076/99‑1 (MQ=255)
               gTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCGGCCCCTGAA‑ACCGTATGCACCTGCTTTATc                                                                           <  2:2030712/100‑1 (MQ=255)
               gTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATc                                                                           >  1:2330939/1‑100 (MQ=255)
               gTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATc                                                                           <  2:2960757/100‑1 (MQ=255)
                   cAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGTCCCTGAA‑ACCGTATGCACCTGCTTTATCTAAc                                                                       >  2:186654/1‑100 (MQ=255)
                    aaCCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACg                                                                      >  2:60884/1‑100 (MQ=255)
                     aCCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACgg                                                                     >  1:156580/1‑100 (MQ=255)
                     aCCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACgg                                                                     <  1:2547540/100‑1 (MQ=255)
                     aCCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACgg                                                                     >  1:101743/1‑100 (MQ=255)
                       cAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGtt                                                                   <  1:3171284/100‑1 (MQ=255)
                          cGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCa                                                                 <  1:3039878/99‑1 (MQ=255)
                           ggTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATc                                                               >  1:2068970/1‑100 (MQ=255)
                              aaGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCt                                                            <  1:772681/100‑1 (MQ=255)
                                gCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTgg                                                          >  1:738942/1‑100 (MQ=255)
                                gCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTgg                                                          >  2:758780/1‑100 (MQ=255)
                                gCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTgg                                                          >  1:932816/1‑100 (MQ=255)
                                gCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTg                                                           >  2:1804635/1‑99 (MQ=255)
                                  cGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCg                                                        <  1:399797/100‑1 (MQ=255)
                                   gACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGa                                                       <  2:237610/100‑1 (MQ=255)
                                       aCGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATa                                                   >  1:2086519/1‑100 (MQ=255)
                                        cGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAg                                                  >  1:2852215/1‑100 (MQ=255)
                                        cGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAg                                                  <  1:2832167/100‑1 (MQ=255)
                                         ggCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGc                                                 <  1:998309/100‑1 (MQ=255)
                                          gCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCGGGctatgtcgcc                                                <  2:232041/100‑10 (MQ=255)
                                          gCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCGGGctatgtcgcc                                                <  2:233940/100‑10 (MQ=255)
                                           cGTGATAGCTGCCATTTATCTTCCTGACAATACTGTCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAg                                                  <  2:2760590/97‑1 (MQ=255)
                                           cGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGccc                                               >  2:1605509/1‑100 (MQ=255)
                                            gTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCcg                                              >  1:2851765/1‑100 (MQ=255)
                                             tGATAACTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCACTGGCGACATAGCCcgc                                             >  2:3073406/1‑100 (MQ=255)
                                                tAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGAt                                          <  1:2911881/100‑1 (MQ=255)
                                                 aGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATg                                         >  2:704715/1‑100 (MQ=255)
                                                  gCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGt                                        >  2:2662077/1‑100 (MQ=255)
                                                  gCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGt                                        >  1:1646616/1‑100 (MQ=255)
                                                   cTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTc                                       >  2:453203/1‑100 (MQ=255)
                                                   cTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTc                                       <  2:722686/100‑1 (MQ=255)
                                                     gCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCtt                                     <  1:2988649/100‑1 (MQ=255)
                                                      ccATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTc                                    <  1:3012533/100‑1 (MQ=255)
                                                          ttATCTTCCTGACAAAACAGCCACCCAGTCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCt                                <  1:186659/100‑1 (MQ=255)
                                                          ttATCTTCCTGACAAAACAGCCACCCAACCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTc                                 <  2:156576/99‑1 (MQ=255)
                                                           tATCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTg                               <  1:780120/100‑1 (MQ=255)
                                                            aTCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGt                              >  2:208495/1‑100 (MQ=255)
                                                             tCTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGtt                             <  2:599456/100‑1 (MQ=255)
                                                             tCTTCCTGACAAAACAACCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGtt                             >  1:103363/1‑100 (MQ=255)
                                                              cTTCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTa                            <  1:985385/100‑1 (MQ=255)
                                                               ttCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTaa                           >  2:2929410/1‑100 (MQ=255)
                                                               ttCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTaa                           <  1:822098/100‑1 (MQ=255)
                                                                tCCTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAAc                          >  1:2189582/1‑100 (MQ=255)
                                                                 ccTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACg                         <  2:1269567/100‑1 (MQ=255)
                                                                  cTGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGt                        >  1:1551713/1‑100 (MQ=255)
                                                                   tGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGtt                       >  2:338186/1‑100 (MQ=255)
                                                                   tGACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGtt                       >  2:1156471/1‑100 (MQ=255)
                                                                    gACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTc                      <  1:1171996/100‑1 (MQ=255)
                                                                    gACAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTc                      >  2:2612716/1‑100 (MQ=255)
                                                                      cAAAACAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTg                    >  1:1530652/1‑100 (MQ=255)
                                                                       aaaaCAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCt                     <  2:1163338/98‑1 (MQ=255)
                                                                        aaaCAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGATATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGaa                  <  1:758672/100‑1 (MQ=255)
                                                                          aCAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGaaa                 <  2:2750612/99‑1 (MQ=255)
                                                                          aCAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAg                <  2:2003360/100‑1 (MQ=255)
                                                                          aCAGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCAGCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAg                <  2:1951179/100‑1 (MQ=255)
                                                                            aGCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTc              <  1:1159984/100‑1 (MQ=255)
                                                                             gCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTca             <  1:1571861/100‑1 (MQ=255)
                                                                             gCCACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTca             <  2:425110/100‑1 (MQ=255)
                                                                              ccACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTcac            >  1:581252/1‑100 (MQ=255)
                                                                               cACCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTcaca           <  1:2594362/100‑1 (MQ=255)
                                                                                aCCCAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACAt          <  1:2174377/100‑1 (MQ=255)
                                                                                 nccAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATc         >  1:325493/2‑100 (MQ=255)
                                                                                  ccAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCg        <  2:249803/100‑1 (MQ=255)
                                                                                   cAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCgg       >  1:1804762/1‑100 (MQ=255)
                                                                                   cAGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCgg       <  2:1807643/100‑1 (MQ=255)
                                                                                    aGCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGt      <  1:1034904/100‑1 (MQ=255)
                                                                                     gCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTg     <  2:101740/100‑1 (MQ=255)
                                                                                     gCCCCTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTg     <  2:1285974/100‑1 (MQ=255)
                                                                                      ccccTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTg     >  1:2837045/1‑100 (MQ=255)
                                                                                      ccccTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGa    <  2:2871313/100‑1 (MQ=255)
                                                                                      ccccTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGa    <  1:1864185/100‑1 (MQ=255)
                                                                                       cccTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGac   <  1:1403223/100‑1 (MQ=255)
                                                                                       cccTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGac   >  1:243159/1‑100 (MQ=255)
                                                                                        ccTGAA‑ACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGacc  >  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: ‑