Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 889,028 T→C 17.6% L264L (CTA→CTG)  ybjL ← putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913889,0280TC17.6% 233.8 / 39.2 107L264L (CTA→CTG) ybjLputative transporter
Reads supporting (aligned to +/- strand):  ref base T (38/50);  new base C (11/8);  total (49/58)
Fisher's exact test for biased strand distribution p-value = 3.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

AGGTCACGATCGAAAACTTCTTTACCGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGCCGGTTTGTCGATAAATACCCAGTTCACGCAGATTT  >  NC_000913/888930‑889127
                                                                                                  |                                                                                                   
aGGTCACGATCGAAAACTTCTTTACCGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTa                                                                                                    <  2:104464/100‑1 (MQ=255)
  gTCACGATCGAAAACTTCTTTACCGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGCAGc                                                                                                  >  1:1998008/1‑100 (MQ=255)
  gTCACGATCGAAAACTTCTTTACCGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGCAGc                                                                                                  <  1:584691/100‑1 (MQ=255)
   tCACGATCGAAAACTTCTTTACCGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGc                                                                                                  >  2:768114/1‑99 (MQ=255)
   tCACGATCGAAAACTTCTTTACCGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGc                                                                                                  >  2:1874124/1‑99 (MQ=255)
    cACGATCGAAGACTTCTTTACCGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGc                                                                                                  <  1:9644/98‑1 (MQ=255)
    cACGATCGAAAACTTCTTTACCGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCAc                                                                                                <  2:1470223/100‑1 (MQ=255)
     aCGATCGAAAACTTCTTTACCGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACg                                                                                               >  2:1804710/1‑100 (MQ=255)
      cGATCGAAAACTTCTTTACCGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACgg                                                                                              <  2:2437068/100‑1 (MQ=255)
      cGATCGAAAACTTCTTTACCGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGCAGCACgg                                                                                              <  2:788633/100‑1 (MQ=255)
        aTCGAAAACTTCTTTACCGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCa                                                                                            <  2:145133/100‑1 (MQ=255)
        aTCGAAAACTTCTTTACCGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGCAGCACGGCa                                                                                            <  2:1077818/100‑1 (MQ=255)
        aTCGAAAACTTCTTTACCGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGCGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCa                                                                                            <  2:691865/100‑1 (MQ=255)
         tCGAAAACTTCTTTACCGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCAt                                                                                           <  2:1538104/100‑1 (MQ=255)
         tCGAAAACTTCTTTACCGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCAt                                                                                           >  1:674850/1‑100 (MQ=255)
          cGAAAACTTCTTTACCGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATc                                                                                          <  2:1886559/100‑1 (MQ=255)
          cGAAAACTTCTTTACCGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATc                                                                                          >  2:750332/1‑100 (MQ=255)
          cGAAAACTTCTTTACCGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATc                                                                                          >  2:615544/1‑100 (MQ=255)
           gAAAACTTCTTTACCGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCa                                                                                         >  2:216932/1‑100 (MQ=255)
           gAAAACTTCTTTACCGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCa                                                                                         <  2:2396985/100‑1 (MQ=255)
              aaCTTCTTTACCGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCg                                                                                      >  2:2353184/1‑100 (MQ=255)
                  tcttTACCGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATTACCGTCTg                                                                                  <  2:1408179/100‑1 (MQ=255)
                  tcttTACCGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTg                                                                                  <  1:795431/100‑1 (MQ=255)
                  tcttTACCGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTg                                                                                  <  2:149231/100‑1 (MQ=255)
                     ttACCGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGAt                                                                               <  1:75965/100‑1 (MQ=255)
                      tACCGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGAtt                                                                              >  1:996183/1‑100 (MQ=255)
                       aCCGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGAttt                                                                             >  1:558311/1‑100 (MQ=255)
                        ccGTTGCGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCCGCCAACGCTATTTCATCGCGCAATTGTAGCACGGCATCAGCGtgtg                                                                                  <  1:797158/94‑1 (MQ=255)
                        ccGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTg                                                                            <  1:2206688/100‑1 (MQ=255)
                        ccGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTg                                                                            <  1:1299109/100‑1 (MQ=255)
                        ccGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTg                                                                            <  2:2545721/100‑1 (MQ=255)
                        ccGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGCAGCACGGCATCACCGTCTGGATTTg                                                                            >  1:1204061/1‑100 (MQ=255)
                         cGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGc                                                                           <  2:703306/100‑1 (MQ=255)
                         cGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGc                                                                           <  2:677358/100‑1 (MQ=255)
                            tACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAg                                                                        <  2:2523888/100‑1 (MQ=255)
                            tACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGCAGCACGGCATCACCGTCTGGATTTGCCAg                                                                        <  1:1142968/100‑1 (MQ=255)
                              cGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGAGAGCCTACCATCGCTATTTCATCGCGCAGTTGTAGCACGGCATCACTGTCTGGATTTg                                                                            <  1:1233378/94‑1 (MQ=255)
                               ggAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAAt                                                                     <  1:1493648/100‑1 (MQ=255)
                                gAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATc                                                                    <  1:216930/100‑1 (MQ=255)
                                 aaGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGCAGCACGGCATCACCGTCTGGATTTGCCAGAATcc                                                                   >  2:2382680/1‑100 (MQ=255)
                                 aaGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGCAGCACGGCATCACCGTCTGGATTTGCCAGAATcc                                                                   >  2:1159774/1‑100 (MQ=255)
                                  aGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCGACCAACGCTATTTCAGCGCCCATTTGCAGGACGGCATCACCGTCTGGATTTGCCAGAATccc                                                                  >  1:995685/1‑100 (MQ=255)
                                    cTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGt                                                                >  1:137373/1‑100 (MQ=255)
                                    cTGGGATCGAGTCGGGCATGGGCGTCGGGATACCCTACCAACGCTATTTCATCGCCCATTTGCAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGt                                                                >  1:1954651/1‑100 (MQ=255)
                                     tGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGtt                                                               >  2:1524335/1‑100 (MQ=255)
                                      gggATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGCAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTa                                                              >  1:2303640/1‑100 (MQ=255)
                                      gggATCGAGGCGGGCAAGGGGGTCGGGATAGCCTACTAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTa                                                              >  1:1236350/1‑100 (MQ=255)
                                       ggATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTAc                                                             >  1:1028999/1‑100 (MQ=255)
                                       ggATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGCAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTAc                                                             >  1:2440169/1‑100 (MQ=255)
                                        gATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACg                                                            >  2:1731691/1‑100 (MQ=255)
                                        gATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACg                                                            >  1:925655/1‑100 (MQ=255)
                                         aTCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGt                                                           >  1:2078356/1‑100 (MQ=255)
                                         aTCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGCAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGt                                                           >  2:1028415/1‑100 (MQ=255)
                                                 gggCATGGGCGTCGGGAGAGCCTACGAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGGTACGTCGAAt                                                      <  1:197233/97‑1 (MQ=255)
                                                  ggCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCTATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGt                                                  <  1:1812043/100‑1 (MQ=255)
                                                  ggCATGGGCGGCGTGATAGCCTAGCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTATGGATTTGCTAGAATCCCGTTACGTCGAATACGt                                                  <  1:594502/100‑1 (MQ=255)
                                                   gCATGGGCGGCGGGGTAGCCTCCCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGGATCCCGTTACGTCGAATACGtt                                                 <  1:755335/100‑1 (MQ=255)
                                                    cATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTc                                                <  2:674853/100‑1 (MQ=255)
                                                    cATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTc                                                >  1:648625/1‑100 (MQ=255)
                                                    cATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTc                                                <  2:733711/100‑1 (MQ=255)
                                                     aTGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCa                                               <  1:1917914/100‑1 (MQ=255)
                                                        ggCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATg                                            <  1:842170/100‑1 (MQ=255)
                                                         gCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGt                                           <  1:1804700/100‑1 (MQ=255)
                                                          cGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTa                                          <  2:1236354/100‑1 (MQ=255)
                                                          cGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTa                                          >  1:791603/1‑100 (MQ=255)
                                                          cGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTa                                          >  1:2375616/1‑100 (MQ=255)
                                                               ggATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTagcagc                                     >  2:2118292/1‑100 (MQ=255)
                                                                   aGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGCCGGt                                 <  2:190035/100‑1 (MQ=255)
                                                                      ctcccAACGTTTTTGCATCGCCCATTTGCGGGACGGGGTCACCGTCTGGATTTGCCATAATCCAGTTAGGTCGAATACGTTCCATGTAGCAGCCGGt                                 <  1:2255116/94‑1 (MQ=255)
                                                                      cTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGCCGGTTTg                              >  1:774695/1‑100 (MQ=255)
                                                                      cTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAAGCCCGTTACGGCGAATACGTTCAATGTAGCAGCCGGTTTg                              <  2:1606569/100‑1 (MQ=255)
                                                                        aCCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGCCGGTTTGTc                            >  1:1408701/1‑100 (MQ=255)
                                                                         ccAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATTTAGCAGCCGGTTTGTc                            >  1:893816/1‑99 (MQ=255)
                                                                         ccAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGCCGGTTTGTCg                           <  1:236655/100‑1 (MQ=255)
                                                                          cAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTCCGTCGAATACGTTCAATGTAGCAGCCGGTTTGTCg                           >  1:320182/1‑99 (MQ=255)
                                                                           aaCGCTATTTCATCGCCCATTTGCAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGCCGGTTTGTCGAt                         <  2:2034253/100‑1 (MQ=255)
                                                                            aCGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGCCGGTTTGTCGATa                        <  2:517332/100‑1 (MQ=255)
                                                                            aCGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGCCGGTTTGTCGATa                        >  1:2462706/1‑100 (MQ=255)
                                                                            aCGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGCCGGTTTGTCGATa                        <  1:2353169/100‑1 (MQ=255)
                                                                             cactactccATCGCCCATGTGCAGCACGGCCTCATCGTCTGGATTTGCCAGAATCCCGTTTCGGCGAATACGTTCAATGTAGGATCCGGTTTGTCg                           <  1:1542389/88‑1 (MQ=255)
                                                                             cGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGCCGGTTTGTCGATaa                       <  1:421790/100‑1 (MQ=255)
                                                                             cGCTATTTCATCGCCCATTTGCAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGCCGGTTTGTCGATaa                       >  2:659022/1‑100 (MQ=255)
                                                                                 aTTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGCCGGTTTGTCGATAAATAc                   >  1:819362/1‑100 (MQ=255)
                                                                                 aTTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGCCGGTTTGTCGATAAATAc                   <  1:1003718/100‑1 (MQ=255)
                                                                                   ttCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGCCGGTTTGTCGATAAATAccc                 >  2:683445/1‑100 (MQ=255)
                                                                                   ttCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGCCGGTTTGTCGATAAATAccc                 <  2:2140679/100‑1 (MQ=255)
                                                                                    tCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGCCGGTTTGTCGATAAATACCCa                >  1:1042254/1‑100 (MQ=255)
                                                                                    tCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGCCGGTTTGTCGATAAATACCCa                <  1:2212731/100‑1 (MQ=255)
                                                                                    tCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGCCGGTTTGTCGATAAATACCCa                <  2:2248377/100‑1 (MQ=255)
                                                                                     cATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGCCGGTTTGTCGATAAATACCCAg               >  1:907352/1‑100 (MQ=255)
                                                                                     cATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGCCGGTTTGTCGATAAATACCCAg               >  1:2120970/1‑100 (MQ=255)
                                                                                      aTCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGCCGGTTTGTCGATAAATACCCAGt              >  1:1008056/1‑100 (MQ=255)
                                                                                      aTCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGCCGGTTTGTCGATAAATACCCAGt              >  1:2491630/1‑100 (MQ=255)
                                                                                      aTCGCCCATTTGCAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGCCGGTTTGTCGATAAATACCCAGt              >  2:900119/1‑100 (MQ=255)
                                                                                         gCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGCCGGTTTGTCGATAAATACCCAGTTCa           >  1:77985/1‑100 (MQ=255)
                                                                                          cccATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGCCGGTTTGTCGATAAATACCCAGTTCAc          >  2:2312534/1‑100 (MQ=255)
                                                                                          cccATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGCCGGTTTGTCGATAAATACCCAGTTCAc          >  2:137013/1‑100 (MQ=255)
                                                                                           ccATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGCCGGTTTGTCGATAAATACCCAGTTCACg         <  2:1304183/100‑1 (MQ=255)
                                                                                                tGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGCCGGTTTGTCGATAAATACCCAGTTCACGCAGAt    >  1:14887/1‑100 (MQ=255)
                                                                                                  tAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGCCGGTTTGTCGATAAATACCCAGTTCACGCAGAttt  >  1:842591/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
AGGTCACGATCGAAAACTTCTTTACCGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGCCGGTTTGTCGATAAATACCCAGTTCACGCAGATTT  >  NC_000913/888930‑889127

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