Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 383,897 C→T 15.9% 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.9% 454.8 / 64.1 183M13I (ATG→ATA) yaiPputative family 2 glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (79/74);  new base T (19/10);  total (99/84)
Fisher's exact test for biased strand distribution p-value = 2.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.04e-01

Only 100 of 187 total aligned reads displayed.

CGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACC  >  NC_000913/383798‑383996
                                                                                                   |                                                                                                   
cgcCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAAc                                                                                                     >  2:969637/1‑100 (MQ=255)
cgcCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAAc                                                                                                     >  1:2743339/1‑100 (MQ=255)
 gcCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACa                                                                                                    >  2:1067322/1‑100 (MQ=255)
   cAGACACGTGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACAtt                                                                                                  <  1:1195885/100‑1 (MQ=255)
    aGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTg                                                                                                 <  1:3537833/100‑1 (MQ=255)
    aGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTg                                                                                                 <  2:4156460/100‑1 (MQ=255)
    aGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTg                                                                                                 <  2:3318862/100‑1 (MQ=255)
     gACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGc                                                                                                >  1:1557334/1‑100 (MQ=255)
        acGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTAt                                                                                             <  2:2596622/100‑1 (MQ=255)
         cGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATg                                                                                            >  1:1287290/1‑100 (MQ=255)
         cGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATg                                                                                            <  2:730614/100‑1 (MQ=255)
          gggCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAATGTTGCTATgg                                                                                           <  2:3588848/100‑1 (MQ=255)
            gCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGTTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTAtn                                                                                            <  1:1589028/97‑2 (MQ=255)
            gCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGAc                                                                                         <  2:3423229/100‑1 (MQ=255)
             ccTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAATGTTGCTATggttt                                                                                        <  2:1558644/100‑4 (MQ=255)
             ccTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACa                                                                                        <  2:883385/100‑1 (MQ=255)
              cTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACAt                                                                                       >  2:3692692/1‑100 (MQ=255)
                tCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATAc                                                                                     <  1:3398516/100‑1 (MQ=255)
                  gTTATACGCAGTTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAATGTTGCTATggttttt                                                                                      <  1:3985654/97‑6 (MQ=38)
                  gTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAg                                                                                   >  2:281933/1‑100 (MQ=255)
                     ataCGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATa                                                                                <  2:663061/100‑1 (MQ=255)
                       aCGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATcaa                                                                              >  2:1914998/1‑100 (MQ=255)
                       aCGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATaaa                                                                              <  1:3202009/100‑1 (MQ=255)
                       aCGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATaaa                                                                              >  2:4056731/1‑100 (MQ=255)
                          cAGGTATAATGTCATCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAatat                                                                           <  1:3494310/100‑1 (MQ=255)
                          cAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAatat                                                                           >  2:2657368/1‑100 (MQ=255)
                            ggTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATcc                                                                         >  2:251962/1‑100 (MQ=255)
                              tataATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAg                                                                       >  1:3043965/1‑100 (MQ=255)
                               ataATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAgg                                                                      >  2:3162147/1‑100 (MQ=255)
                                taATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGt                                                                     >  1:2047578/1‑100 (MQ=255)
                                  aTGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGttt                                                                   <  1:1809573/100‑1 (MQ=255)
                                  aTGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAAAATTGCTATGGACATACAGATAAATATCCAGGttt                                                                   >  1:2320578/1‑100 (MQ=255)
                                    ggCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTc                                                                 >  2:3417017/1‑100 (MQ=255)
                                     gCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCa                                                                <  1:4252700/100‑1 (MQ=255)
                                      cGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCAt                                                               <  1:3342861/100‑1 (MQ=255)
                                      cGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCAt                                                               <  1:3814194/100‑1 (MQ=255)
                                      cGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCAt                                                               <  2:1846386/100‑1 (MQ=255)
                                         cAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCa                                                            <  2:2695098/100‑1 (MQ=255)
                                         cAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCa                                                            <  2:4260176/100‑1 (MQ=255)
                                         cAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCa                                                            >  2:4212282/1‑100 (MQ=255)
                                          aaTACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCaa                                                           <  2:286826/100‑1 (MQ=255)
                                          aaTACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCa                                                            <  1:629846/99‑1 (MQ=255)
                                           aTACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGAACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAAt                                                          >  2:932762/1‑100 (MQ=255)
                                            tACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAAtt                                                         >  2:4054842/1‑100 (MQ=255)
                                            tACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAAtt                                                         <  1:3589093/100‑1 (MQ=255)
                                              cAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTg                                                       >  2:3304955/1‑100 (MQ=255)
                                               aGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGc                                                      <  1:908352/100‑1 (MQ=255)
                                               aGCCTTTCTATTTACTCGGTTTACGTCTTAGCGTACTTAATAACCATAGCAACATTTCTATGGACATACAGATAAATATTCAGGTTTTCATGCAAttttc                                                      <  1:596792/100‑3 (MQ=255)
                                                gCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCt                                                     >  2:308908/1‑100 (MQ=255)
                                                 ccTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTc                                                    <  2:233921/100‑1 (MQ=255)
                                                 ccTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTc                                                    >  2:1100013/1‑100 (MQ=255)
                                                 ccTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTc                                                    >  2:573044/1‑100 (MQ=255)
                                                  ctttctttTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCa                                                   <  2:800470/100‑1 (MQ=255)
                                                  ctttctttTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCa                                                   >  1:2694854/1‑100 (MQ=255)
                                                     tctttTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAg                                                <  2:3464944/100‑1 (MQ=255)
                                                     tctttTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAg                                                <  1:3274747/100‑1 (MQ=255)
                                                        tttGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAAt                                             >  1:3899533/1‑100 (MQ=255)
                                                        tttGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAAt                                             >  1:2197614/1‑100 (MQ=255)
                                                         ttGACTGGGTTTGCATATTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAAt                                             <  1:1928890/99‑1 (MQ=255)
                                                          tGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGACTCa                                           >  1:1632160/1‑100 (MQ=255)
                                                          tGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCa                                           >  2:3894249/1‑100 (MQ=255)
                                                           gACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCAt                                          >  2:3384250/1‑100 (MQ=255)
                                                           gACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCAt                                          >  1:2777776/1‑100 (MQ=255)
                                                           gACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCAt                                          >  1:3023394/1‑100 (MQ=255)
                                                              tGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCAtttt                                       <  2:2847640/100‑1 (MQ=255)
                                                                ggTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTAt                                     >  2:3674208/1‑100 (MQ=255)
                                                                 gTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATg                                    <  1:4113225/100‑1 (MQ=255)
                                                                   ttACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGaa                                  <  1:1520082/100‑1 (MQ=255)
                                                                    tACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAAt                                 >  1:553841/1‑100 (MQ=255)
                                                                       gTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTAc                              <  2:167037/100‑1 (MQ=255)
                                                                         cTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACaa                            <  1:899550/100‑1 (MQ=255)
                                                                            aGCGTACTTAAAAACAATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGc                         <  1:4014720/100‑1 (MQ=255)
                                                                             gCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGcc                        >  2:1093158/1‑100 (MQ=255)
                                                                              cGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCt                       >  2:1300197/1‑100 (MQ=255)
                                                                                 aCTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCa                    >  2:1315612/1‑100 (MQ=255)
                                                                                  cTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCAc                   <  2:2445300/100‑1 (MQ=255)
                                                                                    tAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCAccc                 >  2:1948860/1‑100 (MQ=255)
                                                                                    tAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCAccc                 >  2:2136892/1‑100 (MQ=255)
                                                                                     aaaaaCCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCa                >  2:3638999/1‑100 (MQ=255)
                                                                                      aaaaCCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAg               <  2:4318757/100‑1 (MQ=255)
                                                                                       tgtccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGa              >  1:1655721/4‑100 (MQ=255)
                                                                                       tgtccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGa              >  1:716906/4‑100 (MQ=255)
                                                                                       tgtccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGa              >  1:2896387/4‑100 (MQ=255)
                                                                                        gtccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGAt             >  1:2634533/3‑100 (MQ=255)
                                                                                        gtccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGAt             >  2:525577/3‑100 (MQ=255)
                                                                                        gtccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGAt             >  1:1945576/3‑100 (MQ=255)
                                                                                         tccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATc            >  2:778183/2‑100 (MQ=255)
                                                                                         tccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATc            >  1:1367971/2‑100 (MQ=255)
                                                                                         aCCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATc            >  1:551849/1‑100 (MQ=255)
                                                                                          ccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGAt             >  1:371671/1‑98 (MQ=255)
                                                                                          ccATAGCAACATTGCTATGGACATACAGATAAATATCCAGGCTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGAt             >  1:1447060/1‑98 (MQ=255)
                                                                                           cATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGc          >  1:2665333/1‑100 (MQ=255)
                                                                                           cATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGc          >  1:461061/1‑100 (MQ=255)
                                                                                            aTAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCt         >  1:2465432/1‑100 (MQ=255)
                                                                                            aTAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCt         <  1:1980141/100‑1 (MQ=255)
                                                                                               gCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTAtt      >  2:66158/1‑100 (MQ=255)
                                                                                                cAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTg     <  1:2578117/100‑1 (MQ=255)
                                                                                                cAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTg     <  1:677587/100‑1 (MQ=255)
                                                                                                 aaCATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGa    <  1:729603/100‑1 (MQ=255)
                                                                                                   cATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGAcc  <  2:2266669/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: ‑