Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 383,897 C→T 15.5% M13I (ATG→ATA) ‡ yaiP ← putative family 2 glycosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913383,8970CT15.5% 265.8 / 35.9 110M13I (ATG→ATAyaiPputative family 2 glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (45/48);  new base T (6/11);  total (51/59)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.39e-02

Only 100 of 113 total aligned reads displayed.

CGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACC  >  NC_000913/383798‑383996
                                                                                                   |                                                                                                   
 gcCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTCTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACa                                                                                                    <  1:2317033/100‑1 (MQ=255)
     gACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGc                                                                                                <  1:522539/100‑1 (MQ=255)
       cacGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGTGTACTTAAAAACCATAGCACCATTGCt                                                                                               >  1:425370/1‑99 (MQ=255)
       cacGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTa                                                                                              >  2:962004/1‑100 (MQ=255)
        acGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTAt                                                                                             <  2:32223/100‑1 (MQ=255)
        acGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTAt                                                                                             >  2:1023079/1‑100 (MQ=255)
         cGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATg                                                                                            <  2:799594/100‑1 (MQ=255)
         cGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATg                                                                                            <  2:1386928/100‑1 (MQ=255)
         cGGGCCTTCGTTATACGCAGATATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTAt                                                                                             >  1:1373059/1‑99 (MQ=255)
          gggCCTTCGTTATACGCAGGTATAATGGCGTCAATAGAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATgg                                                                                           >  2:2052840/1‑100 (MQ=255)
          gggCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAATGTTGCTATgg                                                                                           <  1:245332/100‑1 (MQ=255)
          gggCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATgg                                                                                           <  1:105071/100‑1 (MQ=255)
          gggCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATgg                                                                                           <  2:690796/100‑1 (MQ=255)
           ggCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAATGTTGCTATggt                                                                                          <  2:564404/100‑2 (MQ=255)
           ggCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATg                                                                                            >  1:2470690/1‑98 (MQ=255)
           ggCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGa                                                                                          <  1:2125118/100‑1 (MQ=255)
            gCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAATGTTGCTATggtt                                                                                         <  1:899157/100‑3 (MQ=255)
            gCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAATGTTGCTATggtt                                                                                         <  2:2501227/100‑3 (MQ=255)
            gCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGAc                                                                                         <  1:2008164/100‑1 (MQ=255)
             ccTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAATGTTGCTATggttt                                                                                        <  2:887618/100‑4 (MQ=255)
             ccTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACa                                                                                        <  1:381832/100‑1 (MQ=255)
             ccTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACa                                                                                        <  2:671901/100‑1 (MQ=255)
              cTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAATGTTGCTATggtttt                                                                                       <  2:1126697/100‑5 (MQ=255)
               ttCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAATGTTGCTATggttttt                                                                                      <  1:2853619/100‑6 (MQ=255)
               ttCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAATGTTGCTATggttttt                                                                                      <  1:2235072/100‑6 (MQ=255)
               ttCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATa                                                                                      <  1:463924/100‑1 (MQ=255)
                 cGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACa                                                                                    <  1:2052829/100‑1 (MQ=255)
                  gTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAg                                                                                   >  1:2307711/1‑100 (MQ=255)
                  gTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAg                                                                                   <  1:1195405/100‑1 (MQ=255)
                   ttATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAATGTTGCTATggtttttaagt                                                                                  <  2:1037914/100‑10 (MQ=38)
                       aCGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATaaa                                                                              <  2:927941/100‑1 (MQ=255)
                       aCGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGAGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATaaa                                                                              <  1:52011/100‑1 (MQ=255)
                        cGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAat                                                                             >  2:746107/1‑100 (MQ=255)
                        cGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAat                                                                             >  2:2011373/1‑100 (MQ=255)
                         gCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAata                                                                            <  2:2072919/100‑1 (MQ=255)
                         gCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGATTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAata                                                                            <  1:1005989/100‑1 (MQ=255)
                          cAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAatat                                                                           >  2:2060585/1‑100 (MQ=255)
                          cAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAatat                                                                           >  2:728473/1‑100 (MQ=255)
                            ggTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATcc                                                                         <  1:1805560/100‑1 (MQ=255)
                              tataATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAg                                                                       <  2:1625187/100‑1 (MQ=255)
                              tataATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTCACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGCCATACAGATAAATATcc                                                                         >  1:333676/1‑98 (MQ=255)
                              tataATGGCGTCAATACAGCCTTTCTTTTGACTAGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAg                                                                       >  1:2789205/1‑100 (MQ=255)
                                  aTGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGttt                                                                   >  2:2624686/1‑100 (MQ=255)
                                  aTGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGttt                                                                   >  1:847731/1‑100 (MQ=255)
                                   tGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGtttt                                                                  <  1:2345352/100‑1 (MQ=255)
                                    ggCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTc                                                                 >  1:2179236/1‑100 (MQ=255)
                                     gCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCa                                                                <  1:30991/100‑1 (MQ=255)
                                      cGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCAt                                                               >  1:1326298/1‑100 (MQ=255)
                                       gTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCAt                                                               <  2:2755858/99‑1 (MQ=255)
                                       gTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATg                                                              <  2:1129350/100‑1 (MQ=255)
                                        tCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGc                                                             <  1:80327/100‑1 (MQ=255)
                                          aaTACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCaa                                                           <  1:114368/100‑1 (MQ=255)
                                           aTACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAAt                                                          <  1:2024188/100‑1 (MQ=255)
                                             aCAGGCTTTCTTTTGACTGGGTTTCCGTCTTAGCGTACTTAAAAACAATAGCAACATTGCTATGAACATACAGATAAATAACCAGGTTTTCATGCAAttt                                                        <  1:2325786/100‑1 (MQ=255)
                                              cAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAAttt                                                        >  1:79273/1‑99 (MQ=255)
                                               aGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGc                                                      <  1:2331054/100‑1 (MQ=255)
                                               aGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGc                                                      <  2:1656671/100‑1 (MQ=255)
                                                gCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCt                                                     >  2:2066687/1‑100 (MQ=255)
                                                  ctttctttTGACTGGGTTTACGTCTTATCGTACTTAAAAGCCATAGCAACATTGCTATGGTCATACAGATAAATATCCAGGTTTTCATGCAAGTTGCt                                                     >  1:2186085/1‑98 (MQ=255)
                                                  ctttctttTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCa                                                   <  2:707123/100‑1 (MQ=255)
                                                  ctttctttTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCa                                                   >  2:2582385/1‑100 (MQ=255)
                                                   tttctttTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCaa                                                  >  2:408913/1‑100 (MQ=255)
                                                   tttctttTCACTGGTTTTATGTCTTAGCTTACTTTAATACCATAGCAACATCGCTATGGACATACCGATAAACATCCAGGTTTTCATGCAATTTGCTc                                                    <  1:2666152/98‑1 (MQ=255)
                                                      ctttTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGa                                               >  1:49606/1‑100 (MQ=255)
                                                      ctttTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACAGTGCTATAGACATACAGATGAATATCCAGGTTTTCATGCAATTTGCTCaaa                                                 <  1:1058238/98‑1 (MQ=255)
                                                        tttGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAAt                                             <  1:2863018/100‑1 (MQ=255)
                                                          tGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCa                                           >  2:1293435/1‑100 (MQ=255)
                                                           gACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCAt                                          >  2:133681/1‑100 (MQ=255)
                                                             cTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCAttt                                        <  2:1656283/100‑1 (MQ=255)
                                                              tGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCAtttt                                       >  2:2036530/1‑100 (MQ=255)
                                                               gggTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCGAAGGATCATTTTa                                      >  1:2919494/1‑100 (MQ=255)
                                                                 gTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATg                                    <  2:1431508/100‑1 (MQ=255)
                                                                 gTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATg                                    >  2:1583907/1‑100 (MQ=255)
                                                                    tACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCCCAAAGAATCATTTTATGAAt                                 <  1:1785437/100‑1 (MQ=255)
                                                                     aCGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAAtt                                <  1:2864764/100‑1 (MQ=255)
                                                                        tCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACa                             <  2:1200055/100‑1 (MQ=255)
                                                                          ttAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACaaa                           >  2:336368/1‑100 (MQ=255)
                                                                            aGCGTACTTAAAATCCATAGCAACATTTCTATTTATATACAGATAAATATCCAGGTCTTCATGCAATTTCCTCAAAGAATCATTTTATGAATTACAAAGc                         <  1:1918440/100‑1 (MQ=255)
                                                                            aGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGc                         >  2:2183517/1‑100 (MQ=255)
                                                                            aGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGc                         >  1:2532600/1‑100 (MQ=255)
                                                                            aGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGc                         <  2:2006505/100‑1 (MQ=255)
                                                                             gCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGcc                        >  1:2447736/1‑100 (MQ=255)
                                                                             gCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGcc                        <  2:1014387/100‑1 (MQ=255)
                                                                              cGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCt                       <  2:2508518/100‑1 (MQ=255)
                                                                                tACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTc                     >  2:1597252/1‑100 (MQ=255)
                                                                                  cTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCAc                   >  1:1651766/1‑100 (MQ=255)
                                                                                  cTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATAATTTTATGAATTACAAAGCCtt                      >  1:1646934/1‑97 (MQ=255)
                                                                                    tAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCAccc                 >  2:2613456/1‑100 (MQ=255)
                                                                                     aaaaaCCATAGCACCATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCa                <  2:1267309/100‑1 (MQ=255)
                                                                                       tgtccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGa              >  1:777940/4‑100 (MQ=255)
                                                                                       ataCCATAGCAACATTGCTATGGACATACAGATAAATATCCAAGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGa              >  1:1278028/3‑100 (MQ=255)
                                                                                        gtccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGAt             >  1:1623105/3‑100 (MQ=255)
                                                                                        gtccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGAt             >  1:1907502/3‑100 (MQ=255)
                                                                                         tccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATc            >  2:932334/2‑100 (MQ=255)
                                                                                          ccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCg           >  2:2258370/1‑100 (MQ=255)
                                                                                             tAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTa        >  1:1798390/1‑100 (MQ=255)
                                                                                              aGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTAt       <  2:49607/100‑1 (MQ=255)
                                                                                              aGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTAt       >  2:1418966/1‑100 (MQ=255)
                                                                                                   cATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGAcc  <  1:2148531/100‑1 (MQ=255)
                                                                                                   cATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGAcc  <  2:1181641/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACC  >  NC_000913/383798‑383996

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: