Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,115,627 C→T 100% intergenic (‑137/+35) mdtG ← / ← lpxL putative drug efflux system protein/lauryl‑acyl carrier protein (ACP)‑dependent acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,115,6270CT99.2% 423.4 / ‑4.7 119intergenic (‑137/+35)mdtG/lpxLputative drug efflux system protein/lauryl‑acyl carrier protein (ACP)‑dependent acyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (55/63);  total (55/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 119 total aligned reads displayed.

AAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAACCACATATACT  >  NC_000913/1115531‑1115723
                                                                                                |                                                                                                
aagaaAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCgg                                                                                               >  2:1909033/1‑100 (MQ=255)
 agaaAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGa                                                                                              <  2:710528/100‑1 (MQ=255)
   aaaaGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATg                                                                                            >  2:1993547/1‑100 (MQ=255)
    aaaGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGc                                                                                           >  2:2468263/1‑100 (MQ=255)
    aaaGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGc                                                                                           <  2:516057/100‑1 (MQ=255)
    aaaGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGc                                                                                           <  2:3121497/100‑1 (MQ=255)
     aaGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGcc                                                                                          >  1:2625490/1‑100 (MQ=255)
     aaGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGcc                                                                                          <  2:2059157/100‑1 (MQ=255)
       gATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCAt                                                                                        >  1:362615/1‑100 (MQ=255)
       gATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCAt                                                                                        <  1:1902117/100‑1 (MQ=255)
         tAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCTAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTc                                                                                      <  1:3215388/100‑1 (MQ=255)
         tAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTc                                                                                      <  2:2640879/100‑1 (MQ=255)
           aGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTg                                                                                    >  1:1252544/1‑100 (MQ=255)
              acacTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAg                                                                                 <  2:79314/100‑1 (MQ=255)
              acacTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAg                                                                                 <  2:2618923/100‑1 (MQ=255)
              acacTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAg                                                                                 <  2:1422225/100‑1 (MQ=255)
              acacTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAg                                                                                 <  1:267363/100‑1 (MQ=255)
                   aaTTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATcc                                                                            >  1:2873372/1‑100 (MQ=255)
                   aaTTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATcc                                                                            >  1:493845/1‑100 (MQ=255)
                       gatctcccGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCa                                                                        <  1:896102/95‑1 (MQ=255)
                       aTGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCa                                                                        <  2:2305222/100‑1 (MQ=255)
                        tGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCAt                                                                       <  2:2571947/100‑1 (MQ=255)
                         gcgcCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATg                                                                      <  2:165179/100‑1 (MQ=255)
                         gcgcCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATg                                                                      <  2:2657321/100‑1 (MQ=255)
                          cgcCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATgg                                                                     <  2:452290/100‑1 (MQ=255)
                           gcCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATggg                                                                    >  1:1752606/1‑100 (MQ=255)
                           gcCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATggg                                                                    <  1:794475/100‑1 (MQ=255)
                            cccGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCaagatc                                                                   >  2:896105/1‑95 (MQ=255)
                              cGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGgaga                                                                 <  2:1598447/100‑1 (MQ=255)
                               gACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGAt                                                                <  2:1056457/100‑1 (MQ=255)
                               gACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGGGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGgaga                                                                 >  1:1500434/1‑99 (MQ=255)
                                 cTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGAttt                                                              >  1:1448105/1‑100 (MQ=255)
                                  ttCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTa                                                             <  2:2521595/100‑1 (MQ=255)
                                   tCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTaa                                                            <  1:2505395/100‑1 (MQ=255)
                                    ccAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAAt                                                           <  1:3311085/100‑1 (MQ=255)
                                    ccAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAAt                                                           <  1:290180/100‑1 (MQ=255)
                                    ccAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAAt                                                           >  1:1148958/1‑100 (MQ=255)
                                    ccAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAAt                                                           <  2:1239727/100‑1 (MQ=255)
                                    caggggggCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAAt                                                           >  2:2702666/4‑100 (MQ=255)
                                     cAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATa                                                          >  2:470759/1‑100 (MQ=255)
                                     cAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATa                                                          >  2:923603/1‑100 (MQ=255)
                                     cAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATa                                                          <  1:1607087/100‑1 (MQ=255)
                                       gggggCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGc                                                        <  2:919654/100‑1 (MQ=255)
                                        ggggCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCg                                                       <  2:2764657/100‑1 (MQ=255)
                                        ggggCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCg                                                       >  1:47554/1‑100 (MQ=255)
                                            cgcAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGCTTTAATAGCGTGaa                                                   >  1:135081/1‑100 (MQ=255)
                                            cgcAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGaa                                                   >  2:1277845/1‑100 (MQ=255)
                                             gcAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAg                                                  >  1:513275/1‑100 (MQ=255)
                                              cAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAgg                                                 <  2:221772/100‑1 (MQ=255)
                                               aaTCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGa                                                >  2:1218794/1‑100 (MQ=255)
                                               aaTCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGa                                                <  2:3132341/100‑1 (MQ=255)
                                                 tCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAAc                                              >  2:2998966/1‑100 (MQ=255)
                                                  ccAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACg                                             <  1:2431873/100‑1 (MQ=255)
                                                  ccAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATAGGAGATTTAATAGCGTGAAGGAACg                                             >  1:648768/1‑100 (MQ=255)
                                                   cAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGc                                            >  1:1960045/1‑100 (MQ=255)
                                                         gCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCg                                      >  2:1128867/1‑100 (MQ=255)
                                                         gCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCg                                      >  2:1107304/1‑100 (MQ=255)
                                                           ttttATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCggg                                    <  2:3275469/100‑1 (MQ=255)
                                                           ttttATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCggg                                    >  1:3123222/1‑100 (MQ=255)
                                                           ttttATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCggg                                    <  2:637088/100‑1 (MQ=255)
                                                           cttccgatcTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCggg                                    <  1:563435/92‑1 (MQ=255)
                                                            tttATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGc                                   <  2:3287347/100‑1 (MQ=255)
                                                            tttATCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGc                                   >  2:2480465/1‑100 (MQ=255)
                                                               aTCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATTGCGTGAAGGAACGCCTTCCGGGCgtg                                <  2:2040065/100‑1 (MQ=255)
                                                               aTCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCgtg                                <  2:2229542/100‑1 (MQ=255)
                                                               aTCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCgtg                                <  2:493845/100‑1 (MQ=255)
                                                               aTCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCgtg                                >  2:2311929/1‑100 (MQ=255)
                                                               aTCGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCgtg                                >  2:184386/1‑100 (MQ=255)
                                                                 cGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGtt                              >  2:3285964/1‑100 (MQ=255)
                                                                 cGCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGtt                              <  1:2311912/100‑1 (MQ=255)
                                                                  gCTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGttt                             >  2:1652417/1‑100 (MQ=255)
                                                                   cTAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGtttt                            <  1:1993533/100‑1 (MQ=255)
                                                                    tAAATCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGtttt                            <  2:3191627/99‑1 (MQ=255)
                                                                    tAAAGCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATGCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGttttn                           <  1:583312/100‑2 (MQ=255)
                                                                     aaaTCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTaa                          >  1:2021516/1‑100 (MQ=255)
                                                                       aTCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAg                        <  2:2827679/100‑1 (MQ=255)
                                                                       aTCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAg                        <  2:777925/100‑1 (MQ=255)
                                                                       aTCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAg                        <  2:508433/100‑1 (MQ=255)
                                                                        tCAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGc                       <  1:3018441/100‑1 (MQ=255)
                                                                         cAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCg                      >  1:114190/1‑100 (MQ=255)
                                                                         cAGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTCAAAGCg                      >  2:692369/1‑100 (MQ=255)
                                                                          aGGGGGATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGa                     >  2:837565/1‑100 (MQ=255)
                                                                           gggggATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTGATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGAc                    <  2:231364/100‑1 (MQ=255)
                                                                             gggATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACgg                  <  1:1107301/100‑1 (MQ=255)
                                                                             gggATTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACgg                  >  2:1147053/1‑100 (MQ=255)
                                                                                aTTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGgtgt               >  2:750584/1‑100 (MQ=255)
                                                                                aTTTGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGgtg                >  1:745418/1‑99 (MQ=255)
                                                                                 tttGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTa              <  2:1252548/100‑1 (MQ=255)
                                                                                 tttGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTa              >  2:180804/1‑100 (MQ=255)
                                                                                 tttGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTa              <  2:2834706/100‑1 (MQ=255)
                                                                                 tttGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGgtgt               <  1:1334098/99‑1 (MQ=255)
                                                                                   tGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAAc            >  2:1029756/1‑100 (MQ=255)
                                                                                   tGCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAAc            >  2:3053103/1‑100 (MQ=255)
                                                                                   tGCTGTGGTAATGTCGGATGCAATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAAc            >  1:137218/1‑100 (MQ=255)
                                                                                    gCTGTGGTAATGTCGGATGCGATTCTGAAGCACCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAAcc           <  1:824478/100‑1 (MQ=255)
                                                                                    gCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAAcc           <  1:161140/100‑1 (MQ=255)
                                                                                    gCTGTGGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAAcc           <  2:3316293/100‑1 (MQ=255)
                                                                                       gtgGTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAACcaca        <  2:1498988/100‑1 (MQ=255)
                                                                                          gTAATGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAACCACatat     >  1:2692812/1‑100 (MQ=255)
                                                                                            aaTGTCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAACCACATATAc   >  2:2702471/1‑100 (MQ=255)
                                                                                                |                                                                                                
AAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAACCACATATACT  >  NC_000913/1115531‑1115723

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: