Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 383,898 A→G 15.5% M13T (ATG→ACG) ‡ yaiP ← putative family 2 glycosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913383,8980AG15.5% 246.3 / 38.3 110M13T (ATG→ACG) yaiPputative family 2 glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (44/49);  new base G (6/11);  total (50/60)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.93e-01

Only 100 of 112 total aligned reads displayed.

GCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACC  >  NC_000913/383799‑383996
                                                                                                   |                                                                                                  
gcCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTCTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACa                                                                                                    <  1:2317033/100‑1 (MQ=255)
    gACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGc                                                                                                <  1:522539/100‑1 (MQ=255)
      caggggCCGTCGTTATACGCAGGCAGAATGGCCTCAATACAGCCTTCCTATAGACTGGCTTTACGTCTCAGCCTACCTCAACACCATAGCAACa                                                                                                    >  1:1433342/4‑94 (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)
        caggCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATg                                                                                            >  2:1967725/3‑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)
          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)
           gCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAATGTTGCTATggt                                                                                          <  1:2457048/99‑2 (MQ=255)
            ccTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAATGTTGCTATggttt                                                                                        <  2:887618/100‑4 (MQ=255)
            ccTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACa                                                                                        <  1:2281547/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:1195405/100‑1 (MQ=255)
                  ttATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAATGTTGCTATggtttttaagt                                                                                  <  2:1037914/100‑10 (MQ=38)
                  ttATAAGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAATGTTGCTATggtttttaag                                                                                   <  1:1981831/99‑9 (MQ=38)
                    ataCGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATa                                                                                <  1:305018/100‑1 (MQ=255)
                      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: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)
                              ataATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAgg                                                                      <  2:1378928/100‑1 (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)
                                         aaTACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCaa                                                           <  1:114368/100‑1 (MQ=255)
                                          aTACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAAt                                                          <  1:2024188/100‑1 (MQ=255)
                                           tACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAAtt                                                         >  2:2791503/1‑100 (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:2582385/1‑100 (MQ=255)
                                                 ctttctttTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCa                                                   <  2:707123/100‑1 (MQ=255)
                                                  tttctttTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCaa                                                  >  2:408913/1‑100 (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:1583907/1‑100 (MQ=255)
                                                                gTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATg                                    <  2:1431508/100‑1 (MQ=255)
                                                                   tACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCCCAAAGAATCATTTTATGAAt                                 <  1:1785437/100‑1 (MQ=255)
                                                                    aCGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAAtt                                <  1:2864764/100‑1 (MQ=255)
                                                                         ttAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACaaa                           >  2:336368/1‑100 (MQ=255)
                                                                          tAGCGTACTTAAAAACCATAGCAACATTGCTATGGTAATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAg                          >  2:103594/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                        <  2:1014387/100‑1 (MQ=255)
                                                                            gCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGcc                        >  1:2447736/1‑100 (MQ=255)
                                                                             cGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCt                       <  2:2508518/100‑1 (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)
                                                                                      aaaCCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGa              >  2:2123039/1‑100 (MQ=255)
                                                                                       gtccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGAt             >  1:1907502/3‑100 (MQ=255)
                                                                                        tccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATc            >  2:932334/2‑100 (MQ=255)
                                                                                         ccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCg           >  2:2258370/1‑100 (MQ=255)
                                                                                         ccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCg           >  1:660274/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  <  2:1181641/100‑1 (MQ=255)
                                                                                                  cATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGAcc  <  1:2148531/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
GCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACC  >  NC_000913/383799‑383996

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: