Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 383,898 A→G 15.9% 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.9% 422.4 / 68.7 182M13T (ATG→ACG) yaiPputative family 2 glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (77/76);  new base G (19/10);  total (96/86)
Fisher's exact test for biased strand distribution p-value = 1.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.22e-01

Only 100 of 189 total aligned reads displayed.

GCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCT  >  NC_000913/383799‑383997
                                                                                                   |                                                                                                   
gcCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACa                                                                                                     >  2:1067322/1‑100 (MQ=255)
gcCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACa                                                                                                     >  1:2673082/1‑100 (MQ=255)
 ccAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACAt                                                                                                    >  1:1972117/1‑100 (MQ=255)
  cAGACACGTGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACAtt                                                                                                   <  1:1195885/100‑1 (MQ=255)
  cAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACAtt                                                                                                   >  2:737576/1‑100 (MQ=255)
  cAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACAtt                                                                                                   >  1:1571558/1‑100 (MQ=255)
   aGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTg                                                                                                  <  2:3318862/100‑1 (MQ=255)
   aGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTg                                                                                                  <  1:3537833/100‑1 (MQ=255)
    gACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGc                                                                                                 <  1:692882/100‑1 (MQ=255)
       acGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTAt                                                                                              <  2:2596622/100‑1 (MQ=255)
        cGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATg                                                                                             >  1:1287290/1‑100 (MQ=255)
         gggCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAATGTTGCTATgg                                                                                            <  2:3588848/100‑1 (MQ=255)
           gCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAATGTTGCTATggtt                                                                                          <  2:911384/100‑3 (MQ=255)
           gCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAATGTTGCTATggtt                                                                                          <  2:836969/100‑3 (MQ=255)
           gCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGAc                                                                                          <  2:3423229/100‑1 (MQ=255)
           gCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGAc                                                                                          <  1:2258474/100‑1 (MQ=255)
            ccTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAATGTTGCTATggttt                                                                                         <  2:1558644/100‑4 (MQ=255)
             cTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACAt                                                                                        >  2:3692692/1‑100 (MQ=255)
               tCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAATGTTGCTATggttttta                                                                                      <  1:1998790/100‑7 (MQ=39)
                 gTTATACGCAGTTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAATGTTGCTATggttttt                                                                                       <  1:3985654/97‑6 (MQ=38)
                 gTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAg                                                                                    >  2:281933/1‑100 (MQ=255)
                    ataCGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATa                                                                                 <  1:2719868/100‑1 (MQ=255)
                    ataCGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATa                                                                                 <  2:3126886/100‑1 (MQ=255)
                      aCGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATcaa                                                                               >  2:1914998/1‑100 (MQ=255)
                       cGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGATTGAGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATaaa                                                                               <  1:4015364/99‑1 (MQ=255)
                         cAGGTATAATTGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTTCTTAAAAACTATAGCAACATTGCTATGGACATACAGATAAatat                                                                            >  1:3390531/1‑100 (MQ=255)
                         cAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAatat                                                                            >  2:2657368/1‑100 (MQ=255)
                           ggTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATcc                                                                          >  2:251962/1‑100 (MQ=255)
                           ggTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATcc                                                                          >  2:4215906/1‑100 (MQ=255)
                           ggTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATcc                                                                          >  2:2515220/1‑100 (MQ=255)
                              ataATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAgg                                                                       <  1:3137302/100‑1 (MQ=255)
                                 aTGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGttt                                                                    >  2:79197/1‑100 (MQ=255)
                                 aTGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGttt                                                                    <  1:1809573/100‑1 (MQ=255)
                                    gCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCa                                                                 <  1:4252700/100‑1 (MQ=255)
                                     cGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCAt                                                                <  2:1846386/100‑1 (MQ=255)
                                     cGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCAt                                                                <  1:3814194/100‑1 (MQ=255)
                                        cAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCa                                                             <  2:2695098/100‑1 (MQ=255)
                                         aaTACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCa                                                             <  1:629846/99‑1 (MQ=255)
                                           tACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAAtt                                                          <  1:3589093/100‑1 (MQ=255)
                                           tACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAAtt                                                          >  2:4054842/1‑100 (MQ=255)
                                           tACAGCCTTACTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGATATGGACATACAGATAAATATACAGGTTTTCATGCAAtt                                                          >  1:193873/1‑100 (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:3658429/1‑100 (MQ=255)
                                                ccTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTc                                                     >  2:573044/1‑100 (MQ=255)
                                                ccTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTc                                                     >  2:1100013/1‑100 (MQ=255)
                                                   ttctttTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCaaa                                                  >  2:719214/1‑100 (MQ=255)
                                                    tctttTGACTGGTTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAg                                                 <  1:2658686/100‑1 (MQ=255)
                                                    tctttTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAg                                                 <  2:3464944/100‑1 (MQ=255)
                                                      ttttGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGaa                                               >  2:956021/1‑100 (MQ=255)
                                                       tttGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAAt                                              >  1:3899533/1‑100 (MQ=255)
                                                         tGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCa                                            >  2:3894249/1‑100 (MQ=255)
                                                          gACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCAt                                           >  1:2777776/1‑100 (MQ=255)
                                                          gACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCAt                                           >  1:3231279/1‑100 (MQ=255)
                                                           aCTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCAtt                                          <  2:1557339/100‑1 (MQ=255)
                                                            cTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCAttt                                         >  2:4060616/1‑100 (MQ=255)
                                                             tGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCAtttt                                        <  1:1849019/100‑1 (MQ=255)
                                                             tGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCAtttt                                        <  2:2847640/100‑1 (MQ=255)
                                                               ggTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTAt                                      >  2:3674208/1‑100 (MQ=255)
                                                                gTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTGCATGCAATTTGCTCAAAGAATCATTTTATg                                     <  2:337261/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)
                                                                         ttAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACaaa                            >  2:3209873/1‑100 (MQ=255)
                                                                            gCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGcc                         >  2:1093158/1‑100 (MQ=255)
                                                                             cGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCt                        >  2:1300197/1‑100 (MQ=255)
                                                                               tACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTc                      <  1:2126177/100‑1 (MQ=255)
                                                                                 cTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCAc                    <  2:4282355/100‑1 (MQ=255)
                                                                                   tAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCAccc                  <  2:3082375/100‑1 (MQ=255)
                                                                                   tAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCAccc                  >  2:1948860/1‑100 (MQ=255)
                                                                                   tAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCAccc                  >  1:4012598/1‑100 (MQ=255)
                                                                                   tAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCAccc                  >  2:2136892/1‑100 (MQ=255)
                                                                                   tAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCAccc                  <  1:1609135/100‑1 (MQ=255)
                                                                                    aaaaaCCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCa                 >  2:3638999/1‑100 (MQ=255)
                                                                                      tgtccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGa               >  1:2896387/4‑100 (MQ=255)
                                                                                      tgtccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGa               >  1:1655721/4‑100 (MQ=255)
                                                                                       gtccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGAt              >  2:3394976/3‑100 (MQ=255)
                                                                                       gtccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGAt              >  2:525577/3‑100 (MQ=255)
                                                                                       gtccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGAt              >  1:2634533/3‑100 (MQ=255)
                                                                                       gtccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAACCTTCACACAGAt              >  1:3456181/3‑100 (MQ=255)
                                                                                        tccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATc             >  2:52376/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)
                                                                                         ccATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCg            >  2:4200032/1‑100 (MQ=255)
                                                                                         ccATAGCAACATTGCTATGGACATACAGATAAATATCCAGGCTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGAt              >  1:1447060/1‑98 (MQ=255)
                                                                                          cATAGCAATGTTGCTATGGACATACAGATTAATATCCAGGTTTTCATGCAATTTGCTCAAAGACTCATTTTATGAATTACAAAGCCTTCACCCAGATCGc           >  2:2414511/1‑100 (MQ=255)
                                                                                          cATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGc           >  1:461061/1‑100 (MQ=255)
                                                                                          cATAGCAACATTGCTATGGACATACAGATAAATATCCAGGCTTTCATTCAATTTGCTCAAAGAATCATTTTCTGAATTACAAACCCTCCACCCAGATCGc           >  1:3275473/1‑100 (MQ=255)
                                                                                           aTAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCt          >  1:2465432/1‑100 (MQ=255)
                                                                                              gCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTAtt       >  1:3038913/1‑100 (MQ=255)
                                                                                              gCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTAtt       >  2:66158/1‑100 (MQ=255)
                                                                                               cAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTg      <  1:677587/100‑1 (MQ=255)
                                                                                               cAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTg      <  1:2578117/100‑1 (MQ=255)
                                                                                                aaCATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGa     <  1:729603/100‑1 (MQ=255)
                                                                                                  taTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGAcc   <  1:3342416/99‑1 (MQ=255)
                                                                                                  cATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGAcc   <  2:2266669/100‑1 (MQ=255)
                                                                                                   aTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCt  <  1:1551090/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCT  >  NC_000913/383799‑383997

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: