Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,859,191 T→C 52.1% D309G (GAT→GGT)  yidK ← putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,859,1910TC52.1% 42.6 / 282.8 212D309G (GAT→GGT) yidKputative transporter
Reads supporting (aligned to +/- strand):  ref base T (50/51);  new base C (50/59);  total (102/110)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 214 total aligned reads displayed.

ATCACCGCACCAAATAACACCGCGCCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTC  >  NC_000913/3859092‑3859288
                                                                                                   |                                                                                                 
 tCACCGCACCAAATAACACCGCGCCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATc                                                                                                  <  2:2561618/100‑1 (MQ=255)
  cACCGCACCAAATAACACCGCGCCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATcc                                                                                                 >  1:130856/1‑100 (MQ=255)
   aCCGCACCAAATAACACCGCGCCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCt                                                                                                <  1:3234142/100‑1 (MQ=255)
   aCCGCACCAAATAACACCGCGCCGAAGAAACCCACCATTGGCACTGGCAGAACGGTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAAccc                                                                                                 >  2:1165295/1‑99 (MQ=255)
      gCACCAAATAACACCGCGCCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGa                                                                                              >  2:1906679/1‑99 (MQ=255)
      gCACCAAATAACACCGCGCCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGan                                                                                             <  2:77430/100‑2 (MQ=255)
       cACCAAATAATACCGCGCCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATa                                                                                            <  2:2708164/100‑1 (MQ=255)
       cACCAAATAACACCGCGCCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATa                                                                                            >  1:1001065/1‑100 (MQ=255)
       cACCAAATAACACCGCGCCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATa                                                                                            <  2:1815304/100‑1 (MQ=255)
        aCCAAATAACACCTCGCCTAAGAAACCCACCAGTGGCACTGGCAGAACGTTATTTACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATa                                                                                            <  2:3876097/99‑1 (MQ=255)
        aCCAAATAACACCGCGCCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATAc                                                                                           >  2:487189/1‑100 (MQ=255)
        aCCAAATAACACCGCGCCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATAc                                                                                           <  1:2614153/100‑1 (MQ=255)
        aCCAAATAACACCGCGCCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATAc                                                                                           <  1:1805647/100‑1 (MQ=255)
          cAAATAACACCGCGCCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATAc                                                                                           <  2:1410258/98‑1 (MQ=255)
          cAAATAACACCGCGCCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAg                                                                                         >  2:503632/1‑100 (MQ=255)
           aaaTAACACCGCGCCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGa                                                                                        <  1:267006/100‑1 (MQ=255)
            aaTAACACCGCGCCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGAt                                                                                       <  2:1102184/100‑1 (MQ=255)
               aaCACCGCGCCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGaa                                                                                    >  1:1275670/1‑100 (MQ=255)
                acacCGCGCCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGaaa                                                                                   >  1:1400051/1‑100 (MQ=255)
                acacCGCGCCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATCTCGGCTTTCGGTAAACCCTGATACAGATGaaa                                                                                   <  1:3920681/100‑1 (MQ=255)
                 cacCGCGCCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAAc                                                                                  <  1:487193/100‑1 (MQ=255)
                  acCGCGCCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGTCCATGTCGGCTTTCGGTAAATCCTGATACAGATGaa                                                                                    >  2:1846011/1‑97 (MQ=255)
                  acCGCGCCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACg                                                                                 >  2:1286718/1‑100 (MQ=255)
                  acCGCGCCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACg                                                                                 >  1:2498341/1‑100 (MQ=255)
                   ccGCGCCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGc                                                                                <  1:1007749/100‑1 (MQ=255)
                     gcgcCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCaa                                                                              >  1:2097106/1‑100 (MQ=255)
                       gcCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATc                                                                            <  2:2230687/100‑1 (MQ=255)
                       gcCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATc                                                                            <  2:2976599/100‑1 (MQ=255)
                       gcCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATc                                                                            <  2:3902258/100‑1 (MQ=255)
                          gaagaaACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAAc                                                                         <  1:3456510/100‑1 (MQ=255)
                           aagaaACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGTCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAAcc                                                                        <  2:3459312/100‑1 (MQ=255)
                            agaaCCCCACCATTGGCACTGGCAGAACGTGATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGCTAAACCCTGATACAGATGAAACGCAATCAAccc                                                                       <  2:2737720/100‑1 (MQ=255)
                             gaaACCCACCATTGGCACTGGCAGAGCGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCt                                                                      <  2:3071631/100‑1 (MQ=255)
                              aaaCCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTg                                                                     <  1:2427545/100‑1 (MQ=255)
                                aCCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGgc                                                                   <  1:3538805/100‑1 (MQ=255)
                                aCCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGgc                                                                   >  2:3617492/1‑100 (MQ=255)
                                 cccACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGgca                                                                  >  1:2971244/1‑100 (MQ=255)
                                 cccACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGgca                                                                  <  2:2840691/100‑1 (MQ=255)
                                 cccACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGgca                                                                  >  1:2120839/1‑100 (MQ=255)
                                  ccaccaTTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGgcag                                                                 <  1:3288324/100‑1 (MQ=255)
                                  ccaccaTTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTgg                                                                    >  2:4132929/1‑97 (MQ=255)
                                   caccaTTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGgcagc                                                                <  2:901736/100‑1 (MQ=255)
                                     ccaTTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCAc                                                              >  1:2706980/1‑100 (MQ=255)
                                      caTTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCAcn                                                             <  1:1907664/100‑2 (MQ=255)
                                      caTTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCAcc                                                             <  1:465155/100‑1 (MQ=255)
                                        ttGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAg                                                           <  1:1484760/100‑1 (MQ=255)
                                          ggCACTGGCAGAACGTTATTGTCCAGCGTCGGGTAGGCCCTGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCCCCAGTa                                                         <  2:3923485/100‑1 (MQ=255)
                                           gCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTAc                                                        >  2:3277033/1‑100 (MQ=255)
                                            cACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTAcc                                                       >  1:3923683/1‑100 (MQ=255)
                                            cACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTAcc                                                       >  1:3781746/1‑100 (MQ=255)
                                            cACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTCTCTGTAAAACCTCATACAGGTGGAACGCAATTTAAACTGACAGCACCAGTAcc                                                       >  2:1195779/1‑100 (MQ=255)
                                               tGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGc                                                    <  2:1086209/100‑1 (MQ=255)
                                                ggCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCg                                                   <  1:2918460/100‑1 (MQ=255)
                                                ggCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCg                                                   >  1:3626119/1‑100 (MQ=255)
                                                  cAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCggg                                                 <  2:2455497/100‑1 (MQ=255)
                                                  cAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCggg                                                 >  1:3359297/1‑100 (MQ=255)
                                                  cAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCggg                                                 <  2:1750756/100‑1 (MQ=255)
                                                  cAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCggg                                                 <  1:1704728/100‑1 (MQ=255)
                                                    gAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTaagatcggaag                                               >  1:2639382/1‑90 (MQ=255)
                                                       cGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAg                                            <  2:4033564/100‑1 (MQ=255)
                                                       cGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAg                                            <  2:1694134/100‑1 (MQ=255)
                                                       cGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAg                                            <  1:636606/100‑1 (MQ=255)
                                                        gTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGGTGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGc                                           >  2:1938146/1‑100 (MQ=255)
                                                        gTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGc                                           <  1:1754429/100‑1 (MQ=255)
                                                              gACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTc                                     >  1:1927039/1‑100 (MQ=255)
                                                                ccAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAg                                   >  2:3114408/1‑100 (MQ=255)
                                                                ccAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCa                                    <  1:3079667/99‑1 (MQ=255)
                                                                ccAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAg                                   <  1:866799/100‑1 (MQ=255)
                                                                 cAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAg                                   <  2:2770737/99‑1 (MQ=255)
                                                                  aGCGTGGGGTAGTCCATGTCGGCTTTCGGGCAATCCGGATACAGATGAAACGCGATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATGTTCa                                    <  2:818438/97‑1 (MQ=255)
                                                                  aGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCGACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCa                                 >  2:2904067/1‑100 (MQ=255)
                                                                  aGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCa                                 >  1:204579/1‑100 (MQ=255)
                                                                          gTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGtt                         >  1:1623739/1‑100 (MQ=255)
                                                                           tAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTa                        >  2:2064807/1‑100 (MQ=255)
                                                                           tAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTa                        <  1:3944871/100‑1 (MQ=255)
                                                                           atctccATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTa                        <  2:955765/96‑1 (MQ=255)
                                                                             ggCCATCTCCGCTTTCTGTAAATCCTGATACAGATCAACCGCACTCAACCCTGGCCCCACCAGTACCAGCCGGTCGACCATTTTCAGCACCACCGTTatc                      >  2:1338601/1‑98 (MQ=255)
                                                                              gCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACa                     <  1:3884929/100‑1 (MQ=255)
                                                                               ccATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAg                    >  2:540549/1‑100 (MQ=255)
                                                                               ccATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTaagat                    >  1:955771/1‑97 (MQ=255)
                                                                               ccATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTAGAGCATTTTCAGCACCAACGTTAACAg                    >  1:2458513/1‑100 (MQ=255)
                                                                                 aTGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAgcg                  <  1:3449664/100‑1 (MQ=255)
                                                                                  tGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAGCCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAgcgc                 <  1:3657103/100‑1 (MQ=255)
                                                                                  tGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAgcgc                 >  2:3368289/1‑100 (MQ=255)
                                                                                       gCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCt            >  1:3056504/1‑100 (MQ=255)
                                                                                        cTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCtt           >  1:2235819/1‑100 (MQ=255)
                                                                                        cTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCtt           >  1:373085/1‑100 (MQ=255)
                                                                                         tttCGGTAAATCCTGATACAGATGAAACGCAAGCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCttt          <  1:1810108/100‑1 (MQ=255)
                                                                                            cGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGc       <  1:1364361/100‑1 (MQ=255)
                                                                                            cGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGc       <  1:1548956/100‑1 (MQ=255)
                                                                                            cGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGc       <  2:1970294/100‑1 (MQ=255)
                                                                                            cGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGc       >  1:769576/1‑100 (MQ=255)
                                                                                             ggTAAGCCCTGATACAGATGATACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGcc      <  2:2326934/100‑1 (MQ=255)
                                                                                             ggTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGcc      >  2:1311958/1‑100 (MQ=255)
                                                                                             ggTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGcc      >  2:1506243/1‑100 (MQ=255)
                                                                                             ggTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGcc      <  1:2152425/100‑1 (MQ=255)
                                                                                              ntAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTCTCAGCACCGCCGTTCACAGCGCCCCCTTTTGccc     >  2:3788126/2‑100 (MQ=255)
                                                                                               tAAATCCTGATACAGAGGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCt    <  1:227067/100‑1 (MQ=255)
                                                                                                aaaTCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCtt   <  1:456169/100‑1 (MQ=255)
                                                                                                 aaTCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTc  >  1:2194422/1‑100 (MQ=255)
                                                                                                   |                                                                                                 
ATCACCGCACCAAATAACACCGCGCCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTC  >  NC_000913/3859092‑3859288

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: ‑