Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,654,763 A→G I27M (ATA→ATG dctR → Putative LuxR family repressor for dicarboxylate transport

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,654,7630AG98.6% 257.4 / ‑3.8 70I27M (ATA→ATGdctRPutative LuxR family repressor for dicarboxylate transport
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (19/49);  total (20/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

CCCTGTCCTGTTAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGTCGGTTATGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATACAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCATTATTGATACATTAATGAATAATGTATTTCATTCTAATTTTCTCACTCAAATTGAACGATTAAAACCTGTCCATGTCATTATTTTCTCCCCCTTTA  >  NC_000913/3654618‑3654908
                                                                                                                                                 |                                                                                                                                                 
cccTGTCCTGTTAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGTCGGTTATGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGaa                                                                                                                                              <  2:736779/151‑1 (MQ=255)
  ctgtcctgtTAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGTCGGTTATGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACg                                                                                                                                            >  2:761501/1‑151 (MQ=255)
    gtcctgtTAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGTCGGTTATGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACga                                                                                                                                           <  2:753965/150‑1 (MQ=255)
     tcctgtTAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGTCGGTTATGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACga                                                                                                                                           >  2:370656/1‑149 (MQ=255)
      cctgtTAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGTCGGTTATGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACgaag                                                                                                                                         <  2:664820/150‑1 (MQ=255)
      cctgtTAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGTCGGTTATGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACgaag                                                                                                                                         <  1:207248/150‑1 (MQ=255)
      cctgtTAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGTCGGTTATGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACga                                                                                                                                           >  1:690854/1‑148 (MQ=255)
         gtTAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGTCGGTTATGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAagaaga                                                                                                                                     <  1:699434/151‑1 (MQ=255)
           tAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGTCGGTTATGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGaagaag                                                                                                                                      >  1:201457/1‑148 (MQ=255)
           tAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGTCGGTTATGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAagag                                                                                                                                    <  2:30102/150‑1 (MQ=255)
           tAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGTCGGTTATGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAagag                                                                                                                                    <  2:30025/150‑1 (MQ=255)
           tAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGTCGGTTATGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAagag                                                                                                                                    <  1:857551/150‑1 (MQ=255)
           tAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGTCGGTTATGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAagag                                                                                                                                    <  1:283129/150‑1 (MQ=255)
            aaTTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGTCGGTTATGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAgaga                                                                                                                                   <  2:348989/150‑1 (MQ=255)
            aaTTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGTCGGTTATGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATGCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAgaga                                                                                                                                   <  2:543330/150‑1 (MQ=255)
            aaTTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGTCGGGTATGATTCTTATAATTATCAGGGATACGAAGTTCTTCACCGCGATGAAAAACATTCTGAGCAAAGGGAATGTCGTTCATATTCCGAAc                                                                                                                                             >  2:531390/1‑140 (MQ=255)
                aCCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGTCGGTTATGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCg                                                                                                                                <  2:779749/149‑1 (MQ=255)
                aCCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGTCGGTTATGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGAc                                                                                                                              >  2:766949/1‑151 (MQ=255)
                   tttGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGTCGGTTATGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGt                                                                                                                            >  2:295045/1‑150 (MQ=255)
                    ttGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGTCGGTTATGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTa                                                                                                                           >  1:494695/1‑150 (MQ=255)
                             tGATTATAAAGTTAATGTCCGCACCAGGAGTCGGTTATGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATc                                                                                                                 >  1:757248/1‑151 (MQ=255)
                              gATTATAAAGTTAATGTCCGCACCAGGAGTCGGTTATGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCAt                                                                                                                  <  2:82697/149‑1 (MQ=255)
                              gATTATAAAGTTAATGTCCGCACCAGGAGTCGGTTATGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATCa                                                                                                                <  1:205434/151‑1 (MQ=255)
                              gATTATAAAGTTAATGTCCGCACCAGGAGGCGGTTATGTTTCTTATAATTACCCGGTATACAATGTTCTTCACCGCGATGAAAAACACTCTGAGTAAATGTAATTTCGTTCATATGCAGAACGAAGAAGAGATCGGCGTAATGTTGc                                                                                                                    <  2:772076/147‑1 (MQ=255)
                               aTTATATAGTTAATGTCCGCACCAGGAGTCGGTTATGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATc                                                                                                                 >  1:207188/1‑149 (MQ=255)
                               aTTATAAAGTTAATGTCCGCACCAGGAGTCGGTTATGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATCAg                                                                                                               <  2:962157/151‑1 (MQ=255)
                                ttATAAAGTTAATGTCCGCACCAGGAGTCGGTTATGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATc                                                                                                                 >  1:422019/1‑148 (MQ=255)
                                  ataAAGTTAATGTCCGCACCAGGAGTCGGTTATGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAAt                                                                                                            <  2:569263/151‑1 (MQ=255)
                                           aTGTCCGCACCAGGAGTCGGTTATGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTc                                                                                                   <  1:1004490/151‑1 (MQ=255)
                                             gTCCGCACCAGGAGTCGGTTATGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCa                                                                                                  <  1:834194/150‑1 (MQ=255)
                                             gTCCGCACCAGGAGTCGGTTATGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATACAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCat                                                                                                 <  1:982660/151‑1 (MQ=255)
                                               ccGCACCAGGAGTCGGTTATGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCatt                                                                                                <  1:75922/150‑1 (MQ=255)
                                                    ccAGGAGTCGGTTATGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCATTATTGa                                                                                           <  2:440415/150‑1 (MQ=255)
                                                    ccAGGAGTCGGTTATGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCATTATTGAt                                                                                          <  1:721399/151‑1 (MQ=255)
                                                          gTCGGTTATGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCATTATTGATACAtt                                                                                     >  2:14981/1‑150 (MQ=255)
                                                            cGGTTATGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCATTATTGATACATTAAt                                                                                  <  2:1028811/151‑1 (MQ=255)
                                                            cGGTTATGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCATTATTGATACATTAAt                                                                                  <  1:498203/151‑1 (MQ=255)
                                                             ggTTATGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCATTATTGATACATTAATg                                                                                 <  2:4245/151‑1 (MQ=255)
                                                              gTTATGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCATTATTGATACATTAATGa                                                                                <  1:239915/151‑1 (MQ=255)
                                                                tATGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCATTATTGATACATTAATGaa                                                                               <  1:375302/150‑1 (MQ=255)
                                                                 aTGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCTTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCATTATTGATACATTAATGaat                                                                              <  2:783056/150‑1 (MQ=255)
                                                                   gTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCATTATTGATACATTAATGaataa                                                                            <  2:717954/150‑1 (MQ=255)
                                                                      tCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCATTATTGATACATTAATGAATAATGt                                                                         <  1:125493/150‑1 (MQ=255)
                                                                          ataATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCATTATTGATACATTAATGAATAATGTAttt                                                                     <  2:605583/150‑1 (MQ=255)
                                                                          ataATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCATTATTGATACATTAATGAATAATGTATTTc                                                                    <  1:837611/151‑1 (MQ=255)
                                                                                         aCGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCATTATTGATACATTAATGAATAATGTATTTCATTCTAATTTtctc                                                      <  1:444568/150‑1 (MQ=255)
                                                                                         aCGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCATTATTGATACATTAATGAATAATGTATTTCATTCTAATTTTctca                                                     <  1:677807/151‑1 (MQ=255)
                                                                                           gATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCATTATTGATACATTAATGAATAATGTATTTCATTCTAATTTTctcact                                                   <  1:136286/151‑1 (MQ=255)
                                                                                            aTGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCATTATTGATACATTAATGAATAATGTATTTCATTCTAATTTTctcacac                                                  <  2:337672/151‑3 (MQ=255)
                                                                                                tcttcACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGGTCATATGCAGAACGAAGAAGAGATCGACGGAATGCTGCATCAGAATGCCTTCGTCATTATTGATACATTAATGAATAATGTATTTCATTCTAATTTTCTCACTCAAAt                                              >  1:501168/1‑151 (MQ=255)
                                                                                                  ttcACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCATTATTGATACATTAATGAATAATGTATTTCGTTCTAATTTTCTCACTCAAATTg                                            >  2:646545/1‑151 (MQ=255)
                                                                                                   tcACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCATTATTGATACATTAATGAATAATGTATTTCATTCTAATTTTCTCACTCAAAtt                                             <  2:690820/149‑1 (MQ=255)
                                                                                                   tcACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCATTATTGATACATTAATGAATAATGTATTTCATTCTAATTTTCTCACTCAAAt                                              <  2:291451/148‑1 (MQ=255)
                                                                                                    cACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCATTATTGATACATTAATGAATAATGTATTTCATTCTAATTTTCTCACTCAAATTGaa                                          <  1:650822/151‑1 (MQ=255)
                                                                                                     aCCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCATTATTGATACATTAATGAATAATGTATTTCATTCTAATTTTCTCACTCAAATTGAAc                                         <  2:1051989/151‑1 (MQ=255)
                                                                                                      ccGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCATTATTGATACATTAATGAATAATGTATTTCATTCTAATTTTCTCACTCAAATTGAACg                                        <  1:370682/151‑1 (MQ=255)
                                                                                                         cgATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCATTATTGATACATTAATGAATAATGTATTTCATTCTAATTTTCTCACTCAAATTGAACGAtt                                     >  1:750786/1‑151 (MQ=255)
                                                                                                             gAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCATTATTGATACATTAATGAATAATGTATTTCATTCTAATTTTCTCACTCAAATTGAACGATTaaaa                                 <  1:744840/151‑1 (MQ=255)
                                                                                                             gAAAAACAGTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAATAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCATTATTGATACATTAATGAATAATGTATTTCATTCTAATTTTCTCACTCAAATTGAACGATTaa                                   <  2:551242/149‑1 (MQ=255)
                                                                                                              aaaaaCATTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCATTATTGATACATTAATGAATAATGTATTTCATTCTAATTTTCTCACTCAAATTGAACGATTAAAAc                                >  2:982089/1‑151 (MQ=255)
                                                                                                                    aTTCTGAGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCATTATTGATACATTAATGAATAATGTATTTCATTCTAATTTTCTCACTCAAATTGAACGATTAAAACCTGTcc                          >  1:149246/1‑151 (MQ=255)
                                                                                                                          aGTAAAGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCATTATTGATACATTAATGAATAATGTATTTCATTCTAATTTTCTCACTCAAATTGAACGATTAAAACCTGTCCATGTc                     >  1:1060572/1‑150 (MQ=255)
                                                                                                                             aaaGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCATTATTGATACATTAATGAATAATGTATTTCATTCTAATTTTCTCACTCAAATTGAACGATTAAAACCTGTCCATGTCatt                  <  2:570284/150‑1 (MQ=255)
                                                                                                                               aGGTAATGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCATTATTGATACATTAATGAATAATGTATTTCATTCTAATTTTCTCACTCAAATTGAACGATTAAAACCTGTCCATGTCattatt               <  1:135397/151‑1 (MQ=255)
                                                                                                                                   aaTTTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCATTATTGATACATTAATGAATAATGTATTTCATTCTAATTTTCTCACTCAAATTGAACGATTAAAACCTGTCCATGTCATTATTTtct           >  2:302597/1‑151 (MQ=255)
                                                                                                                                    aTGTCGTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCATTATTGATACATTAATGAATAATGTATTTCATTCTAATTTTCTCACTCAAATTGAACGATTAAAACCTGTCCATGTCATTATTTtctc          >  2:295616/1‑151 (MQ=255)
                                                                                                                                         gTTCATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATCTGAATGCCTTCGTCATTATTGATACATTAATGAATAATGTATTTCATTCTAATTTTCTCACTCAAATTGAACGATTAAATCCTGTCCATGTCATTATTTTCTCCCCCt     <  2:927826/151‑1 (MQ=255)
                                                                                                                                            cATATGCAGATCGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCATTATTGATGCGTTAATGAATAATGTATTTCATTCTAATTTTCTCTCTCAAATTGAACGATTAAAACCTGTCCATGTCATTATTTTCTCCCCCGTTa  <  2:704156/151‑1 (MQ=255)
                                                                                                                                            cATATGCAGAACGAAGAAGAGATCGTCGTAATGTTGCATCAGAATGCCTTCGTCATTATTGATACATTAATGAATAATGTATTTCATTCTAATTTTCTCACTCAAATTGAACGATTAAAACCTGTCCGTGTGATTATTTTCTCCCCCtt  <  2:653878/150‑2 (MQ=255)
                                                                                                                                            cATATGCAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCATTATTGATACATTAATGAATAATGTATTTCATTCTAATTTTCTCACTCAAATTGAACGATTAAAACCTGTCCATGTCATTATTTTCTCCCCCTTTa  <  1:989507/151‑1 (MQ=255)
                                                                                                                                              tatgCAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCATTATTGATACATTAATGAATAATGTATTTCATTCTAATTTTCTCACTCAAATTGAACGATTAAAACCTGTCCa                         >  2:546898/1‑126 (MQ=255)
                                                                                                                                              tatgCAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCATTATTGATACATTAATGAATAATGTATTTCATTCTAATTTTCTCACTCAAATTGAACGATTAAAACCTGTCCa                         <  1:546929/126‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CCCTGTCCTGTTAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGTCGGTTATGTTTCTTATAATTACCAGGGATACGATGTTCTTCACCGCGATGAAAAACATTCTGAGTAAAGGTAATGTCGTTCATATACAGAACGAAGAAGAGATCGACGTAATGTTGCATCAGAATGCCTTCGTCATTATTGATACATTAATGAATAATGTATTTCATTCTAATTTTCTCACTCAAATTGAACGATTAAAACCTGTCCATGTCATTATTTTCTCCCCCTTTA  >  NC_000913/3654618‑3654908

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