Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,397,925 Δ1 bp 79.9% 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.79.9% 344.0 / 118.4 134coding (454/594 nt)yhdEdTTP/UTP pyrophosphatase; m(5)UTP/m(5)CTP/pseudo‑UTP pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base T (12/15);  new base . (56/51);  total (68/66)
Fisher's exact test for biased strand distribution p-value = 5.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 135 total aligned reads displayed.

AATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACC  >  NC_000913/3397826‑3398019
                                                                                                   |                                                                                              
aaTTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAt                                                                                                <  2:1852354/100‑1 (MQ=255)
    aCTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCg                                                                                            >  2:2481219/1‑100 (MQ=255)
     cTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTa                                                                                          >  1:506241/1‑100 (MQ=255)
     cTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTa                                                                                          <  2:1229806/100‑1 (MQ=255)
       taataaTTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACGGTa                                                                                          >  2:898150/1‑99 (MQ=255)
       taataaTTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTAt                                                                                         >  2:1854506/1‑100 (MQ=255)
       taataaTTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATg                                                                                        >  2:2287741/1‑100 (MQ=255)
       taataaTTCATAAGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTAt                                                                                         >  2:898664/1‑100 (MQ=255)
          taatTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCAc                                                                                     >  1:2888298/1‑100 (MQ=255)
          taatTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCAc                                                                                     >  2:1805773/1‑100 (MQ=255)
           aatTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCAcc                                                                                    <  2:2632592/100‑1 (MQ=255)
           aatTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCAcc                                                                                    <  1:1805767/100‑1 (MQ=255)
            atTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCt                                                                                   <  2:2599395/100‑1 (MQ=255)
             ttCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCt                                                                                   <  2:1931912/100‑1 (MQ=255)
             ttCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTg                                                                                  <  1:3430/100‑1 (MQ=255)
             ttCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTg                                                                                  <  2:2441063/100‑1 (MQ=255)
             ttCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCATATGCACCTg                                                                                  <  2:2464704/100‑1 (MQ=255)
              tCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCATGACAAAACAGCCACCCAGCCCCTGAATTCCGTATGCAccgg                                                                                  <  2:1961552/100‑3 (MQ=255)
               cATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCt                                                                                <  1:1872501/100‑1 (MQ=255)
                aTACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCt                                                                                <  2:2802338/100‑1 (MQ=255)
                 tACGTTTCATCCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCACTTATCTTCCTGACAACACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCt                                                                                <  1:2380676/99‑1 (MQ=255)
                 tACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCtt                                                                               >  2:1680281/1‑100 (MQ=255)
                 tACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGACACCCAGCCCCTGAAACCGTATGCACCTGCttt                                                                              >  1:2008904/1‑100 (MQ=255)
                  aCGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTa                                                                             <  1:1794196/100‑1 (MQ=255)
                   cGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTAt                                                                            >  2:1789765/1‑100 (MQ=255)
                    gTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATc                                                                           >  2:128455/1‑100 (MQ=255)
                     tttCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTAt                                                                            <  2:2843159/99‑1 (MQ=255)
                      ttCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCt                                                                          >  1:9390/1‑100 (MQ=255)
                        cAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAAc                                                                       >  1:826367/1‑100 (MQ=255)
                        cAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAAc                                                                       <  2:192575/100‑1 (MQ=255)
                        cAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAAc                                                                       >  1:1543285/1‑100 (MQ=255)
                          aCCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACgg                                                                     <  2:2405890/100‑1 (MQ=255)
                            cAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGtt                                                                   <  1:1018781/100‑1 (MQ=255)
                               cGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCAt                                                                <  1:2663410/100‑1 (MQ=255)
                                ggTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATc                                                               >  1:2381880/1‑100 (MQ=255)
                                 gTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCg                                                              <  1:1066505/100‑1 (MQ=255)
                                   aaGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCt                                                            >  1:2316585/1‑100 (MQ=255)
                                    aGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCCCCTGCTTTATCTAACGGTTCATCGCTg                                                           >  1:992206/1‑100 (MQ=255)
                                    aGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTg                                                           >  1:1244964/1‑100 (MQ=255)
                                    aGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTg                                                           >  2:1887126/1‑100 (MQ=255)
                                    aGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTg                                                           <  1:2434416/100‑1 (MQ=255)
                                     gCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTgg                                                          >  1:689564/1‑100 (MQ=255)
                                     gCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTgg                                                          >  1:2317301/1‑100 (MQ=255)
                                      tcGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGc                                                         <  2:2486232/99‑1 (MQ=255)
                                       cGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGc                                                         >  2:2569447/1‑100 (MQ=255)
                                       cGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCg                                                        <  2:366768/100‑1 (MQ=255)
                                       cGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCg                                                        <  1:2102819/100‑1 (MQ=255)
                                        gACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCg                                                        <  1:179321/99‑1 (MQ=255)
                                         aCTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGAc                                                      <  1:1762913/100‑1 (MQ=255)
                                          cTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGAc                                                      <  2:515519/100‑1 (MQ=255)
                                          cTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACa                                                     >  2:1066985/1‑100 (MQ=255)
                                              ggCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGc                                                 >  2:429655/1‑100 (MQ=255)
                                               gCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGcc                                                >  1:1840165/1‑100 (MQ=255)
                                               gCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGcc                                                <  1:375205/100‑1 (MQ=255)
                                                  tGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCcgc                                             <  2:287947/100‑1 (MQ=255)
                                                   gATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCcgc                                             <  2:2156422/100‑1 (MQ=255)
                                                       gCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATg                                         >  2:1332870/1‑100 (MQ=255)
                                                        cTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTc                                       >  2:764717/1‑100 (MQ=255)
                                                        cTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTc                                       >  2:2372815/1‑100 (MQ=255)
                                                         tGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCt                                      >  2:793438/1‑100 (MQ=255)
                                                          gCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCtt                                     <  2:1894456/100‑1 (MQ=255)
                                                           ccATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCt                                      >  1:2230498/1‑99 (MQ=255)
                                                            cATTTATCTTCCTGACAAGACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCtt                                     <  1:2664121/98‑1 (MQ=255)
                                                            cATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCg                                   <  2:2523155/100‑1 (MQ=255)
                                                            cATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCg                                   <  2:2017434/100‑1 (MQ=255)
                                                            cATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCg                                   >  1:1455485/1‑100 (MQ=255)
                                                            cATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCg                                   <  1:910548/100‑1 (MQ=255)
                                                             aTTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTc                                    >  1:2312109/1‑99 (MQ=255)
                                                              tttATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTc                                 >  2:320144/1‑100 (MQ=255)
                                                              tttATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTc                                 <  1:879537/100‑1 (MQ=255)
                                                                 aTCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTg                               <  2:1136239/100‑1 (MQ=255)
                                                                   cTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTa                            <  1:762377/100‑1 (MQ=255)
                                                                   cTTCCTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTa                            >  1:2508520/1‑100 (MQ=255)
                                                                      ccTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACg                         <  2:1510272/100‑1 (MQ=255)
                                                                       cTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGt                        >  2:2186901/1‑100 (MQ=255)
                                                                       cTGACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGATTCATCGCTGGCGACATAGCCCGCGATGGCTTCGTCTGTTAAc                          >  1:673202/1‑98 (MQ=255)
                                                                         gACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTc                      >  1:190770/1‑100 (MQ=255)
                                                                         gACAAAACAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTc                      <  2:2140825/100‑1 (MQ=255)
                                                                             aaaCAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGaa                  <  2:573080/100‑1 (MQ=255)
                                                                              aaCAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGaaa                 >  1:1676511/1‑100 (MQ=255)
                                                                              aaCAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTAGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGaaa                 >  1:2039961/1‑100 (MQ=255)
                                                                               aCAGCCACCCAGCCCCTGAAAGCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAg                <  1:2497801/100‑1 (MQ=255)
                                                                               aCAGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAg                <  2:703699/100‑1 (MQ=255)
                                                                                cAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAg                <  1:1332868/100‑1 (MQ=255)
                                                                                 aGCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTc              <  2:2493344/100‑1 (MQ=255)
                                                                                  gCCACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTca             >  2:1029165/1‑100 (MQ=255)
                                                                                    cACCCAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTcaca           >  1:2217835/1‑100 (MQ=255)
                                                                                      cgcAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTAGTCTGTTAACGTTCTGAAAGttccat          >  1:17173/3‑94 (MQ=255)
                                                                                        cAGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCgg       >  2:1632325/1‑100 (MQ=255)
                                                                                         aGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCgg       <  2:1215454/100‑1 (MQ=255)
                                                                                         aGCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGt      >  2:2505138/1‑100 (MQ=255)
                                                                                          gCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGt      >  1:566960/1‑100 (MQ=255)
                                                                                          gCCCCTGAAGCCGTATGCGCCTGCGTTATCTATCGGTTCATCGCTGGCGACATAGCCCGAGATGACTTCGCCTGGTAACGTTCTGAAAGTCACATCGGTg     <  1:824189/100‑1 (MQ=255)
                                                                                          gCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTg     >  1:2573232/1‑100 (MQ=255)
                                                                                          gCCCCTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTg     >  2:772886/1‑100 (MQ=255)
                                                                                           ccccTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGa    >  2:2656023/1‑100 (MQ=255)
                                                                                            cccTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGac   <  2:1840173/100‑1 (MQ=255)
                                                                                             ccTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGacc  >  1:2593914/1‑100 (MQ=255)
                                                                                             ccTGAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGacc  <  2:1393707/100‑1 (MQ=255)
                                                                                               tGAATACCGTATGCACCTGCTTTATCTTACGGTTCATCGCTGCCGGCATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGacc  <  1:443897/99‑1 (MQ=255)
                                                                                                   |                                                                                              
AATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACC  >  NC_000913/3397826‑3398019

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: