Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,597,846 (T)5→4 100% intergenic (‑286/‑138) livK ← / → panM leucine transporter subunit/PanD autocleavage accelerator, panothenate synthesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,597,8420T.97.1% 452.8 / 5.4 103intergenic (‑282/‑142)livK/panMleucine transporter subunit/PanD autocleavage accelerator, panothenate synthesis
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base . (63/37);  total (64/39)
Fisher's exact test for biased strand distribution p-value = 5.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.33e-01

Only 100 of 111 total aligned reads displayed.

TTTATTGATAGATAAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGCCTCCGATTAAAT  >  NC_000913/3597743‑3597937
                                                                                                   |                                                                                               
tttATTGATAGATAAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACAt                                                                                                 >  1:243689/1‑100 (MQ=255)
 ttATTGATAGATAAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACAtt                                                                                                >  1:2165018/1‑100 (MQ=255)
  tATTGATAGTTAAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGAGGAATAATTCACTATCATTCAGGGGATTATGCTGGACAttt                                                                                               >  1:2524691/1‑100 (MQ=255)
   aTTGATAGATAAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTAACTATCATTCAGGGGATTATGCTGGACAtttt                                                                                              >  1:464396/1‑100 (MQ=255)
     tGATAGATAAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACAttttca                                                                                            <  2:2029102/100‑3 (MQ=255)
      gatagataAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACAttttcat                                                                                           <  1:1157807/100‑4 (MQ=255)
      gatagataAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACAttttcat                                                                                           <  2:1835024/100‑4 (MQ=255)
        tagataAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATtc                                                                                        <  1:445628/100‑1 (MQ=255)
          gataAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATtctc                                                                                      >  1:2280507/1‑100 (MQ=255)
          gataAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATtctc                                                                                      >  1:2319892/1‑100 (MQ=255)
          gataAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATtctc                                                                                      >  1:1418246/1‑100 (MQ=255)
          gataAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATtctc                                                                                      >  1:696221/1‑100 (MQ=255)
          gataAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATtct                                                                                       >  2:2147408/1‑99 (MQ=255)
           ataAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTctct                                                                                     >  1:1469335/1‑100 (MQ=255)
           ataAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTctct                                                                                     <  1:192752/100‑1 (MQ=255)
           ataAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTctct                                                                                     >  1:1897683/1‑100 (MQ=255)
           ataAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTctct                                                                                     >  2:1605645/1‑100 (MQ=255)
           ataAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTctct                                                                                     >  1:1016328/1‑100 (MQ=255)
              aaCAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAAt                                                                                  >  1:601429/1‑100 (MQ=255)
              aaCAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAAt                                                                                  <  1:777117/100‑1 (MQ=255)
              aaCAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAAt                                                                                  >  2:2746084/1‑100 (MQ=255)
              aaCAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAAt                                                                                  <  2:601425/100‑1 (MQ=255)
               aCAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATg                                                                                 >  1:1117113/1‑100 (MQ=255)
                cAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGt                                                                                >  1:2254454/1‑100 (MQ=255)
                  gAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGttt                                                                              <  2:2096210/100‑1 (MQ=255)
                  gAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGttt                                                                              >  1:2948126/1‑100 (MQ=255)
                   aaaaaaGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGtttt                                                                             <  2:3016648/100‑1 (MQ=255)
                    aaaaaGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTa                                                                            >  2:1276210/1‑100 (MQ=255)
                     aaaaGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTTCATTCTCTAATGTTTTa                                                                            >  1:2331210/1‑100 (MQ=255)
                         gTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAAtttt                                                                       >  2:90692/1‑100 (MQ=255)
                           gCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGAGATAATGCTGAACATTTTCATTCTCTAATGTTTTAATTTTGt                                                                     >  1:493732/1‑100 (MQ=255)
                            ccTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTa                                                                    >  2:265939/1‑100 (MQ=255)
                             cTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTTCATTCTCTAATGTTTTAATTTTGTa                                                                    <  1:1856516/100‑1 (MQ=255)
                                     gCATAAAATAACGGCACAAAGGGCGGCATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCt                                                           >  1:1851346/1‑100 (MQ=255)
                                     gCATAAAATAACGGCACAAAGGGCGGAATAATTCCCTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCt                                                           >  2:404179/1‑100 (MQ=255)
                                     gCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTACTGTTTTAATTTTGTAATTATTGCt                                                           >  2:1891047/1‑100 (MQ=255)
                                     gCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCt                                                           <  2:1064146/100‑1 (MQ=255)
                                     gCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCt                                                           <  2:2652566/100‑1 (MQ=255)
                                     gCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCt                                                           <  2:456770/100‑1 (MQ=255)
                                     gCATAAAATAACGGCACAAAGAGCGGAATACTTCACTATAATTCAGGGGATTATGCTGGACATTTTCATCATCTAATGTTTTAATTTTGTAATTATTGc                                                            >  2:569833/1‑99 (MQ=255)
                                      cATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTg                                                          >  1:2033160/1‑100 (MQ=255)
                                        tAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGtt                                                        >  2:3024720/1‑100 (MQ=255)
                                        tAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGtt                                                        >  1:1821966/1‑100 (MQ=255)
                                        tAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGt                                                         >  2:1507443/1‑99 (MQ=255)
                                        tAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATATTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGtt                                                        >  1:314453/1‑100 (MQ=255)
                                          aaaTAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTATAATTTTGTAATTATTGCTGTTaa                                                      >  1:1532257/1‑100 (MQ=255)
                                           aaTAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGTACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTa                                                       >  2:68416/1‑98 (MQ=255)
                                           aaTAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCAGTCTCTCATGTTTTACTTTTGTAATTATTGCTGTTaaa                                                     >  2:2042662/1‑100 (MQ=255)
                                             tAACGGCTCAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTa                                                       <  1:656830/96‑1 (MQ=255)
                                              aaCGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTaaaaaa                                                  <  1:2444873/100‑1 (MQ=255)
                                              aaCGGAACAAAGGGCGGAATAATTCACTATCATTTAGGGGATTATTCTTAACATTTTCATTCTCTAATGTTTTAATTTTTTAATTATTGCTGTTaaaaaa                                                  <  2:848466/100‑1 (MQ=255)
                                               aCGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAAt                                                 <  2:2017284/100‑1 (MQ=255)
                                               aCGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAAt                                                 >  1:2219187/1‑100 (MQ=255)
                                                cGGCACAAAGGGCGGAATAATTCACTATCATTCAGGTTATTATGCTGGATATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAAtt                                                <  1:449260/100‑1 (MQ=255)
                                                cGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAAtt                                                <  1:1218126/100‑1 (MQ=255)
                                                cGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAAtt                                                <  2:243686/100‑1 (MQ=255)
                                                  gCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTaa                                              >  2:2817832/1‑100 (MQ=255)
                                                   cacaAATGGCGCAATAATTCACTATCATTCAGGGTATTATGCTGGACATTTTCATTTTCTAATGTTTTAATTTTTTAATTATTGCTGTTAAAAAAt                                                 <  2:2606654/96‑1 (MQ=255)
                                                      aaaGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCAc                                          <  1:1374650/100‑1 (MQ=255)
                                                       aaGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCAcc                                         >  1:2263767/1‑100 (MQ=255)
                                                         gggCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTg                                       >  2:1258984/1‑100 (MQ=255)
                                                          ggCGGAATAATTCACTATCATTCAGGGGCTTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCt                                        >  1:817368/1‑98 (MQ=255)
                                                            cGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCa                                    >  1:527832/1‑100 (MQ=255)
                                                            cGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCa                                    >  1:1169955/1‑100 (MQ=255)
                                                                 taatTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGa                               <  1:566108/100‑1 (MQ=255)
                                                                 taatTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGa                               >  1:1283508/1‑100 (MQ=255)
                                                                   atTCACTATCATTCAGGGGATTATGATGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATAGTTGTGGTTAAAAAATTAATTACCAACTAAAAGa                               >  2:2712729/1‑98 (MQ=255)
                                                                     tCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATa                           <  2:639311/100‑1 (MQ=255)
                                                                     tCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATa                           <  2:852419/100‑1 (MQ=255)
                                                                     tCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATa                           <  2:852753/100‑1 (MQ=255)
                                                                     tCAATATCACTCAGTGCATTATGTTTGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGTTGTTAATTTTTTAATCACCTGCCAAAAGAAATa                           <  2:353379/100‑1 (MQ=255)
                                                                       aCTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATaaa                         >  2:184262/1‑100 (MQ=255)
                                                                        cTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCTAAAAGAAATaaaa                        >  2:2061627/1‑100 (MQ=255)
                                                                          aTCATTCAGGTGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATaaaaaa                      <  1:679796/100‑1 (MQ=255)
                                                                           tCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAATTTAATCACCTGCCAAAAGAAATAAAAAag                     <  1:337625/100‑1 (MQ=255)
                                                                           tCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAag                     <  1:394032/100‑1 (MQ=255)
                                                                             aTTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAagag                   <  1:2201256/100‑1 (MQ=255)
                                                                             aTTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAagag                   >  2:228405/1‑100 (MQ=255)
                                                                             aTTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAagag                   <  2:1723770/100‑1 (MQ=255)
                                                                             aTTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAagag                   >  2:1434604/1‑100 (MQ=255)
                                                                               tCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGaa                 >  2:2532427/1‑100 (MQ=255)
                                                                                cAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTCATCACCTGCCAAACGAAATAAAAAAGAGaaa                >  2:2390362/1‑100 (MQ=255)
                                                                                cAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGaaa                >  2:2895922/1‑100 (MQ=255)
                                                                                 agaggATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCAAAAAAAAAAAAAAAAGAGAAAg               <  2:1093577/97‑1 (MQ=255)
                                                                                  ggggATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGc              >  1:608816/1‑100 (MQ=255)
                                                                                  ggggATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGc              >  1:1568419/1‑100 (MQ=255)
                                                                                  ggggATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAAATAATCACATACAAAAAGAAATAAAAAAGAGAAAg               >  1:988401/1‑99 (MQ=255)
                                                                                   gggATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGcc             <  1:2792735/100‑1 (MQ=255)
                                                                                        tATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGCCTCCGa        <  1:1457322/100‑1 (MQ=255)
                                                                                          tGCTGGATATTTTCATTCTCTAATTTTTTAATTTTGTAATTATGGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGCCTCCg         <  2:508424/97‑1 (MQ=255)
                                                                                          tGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGCCTCCGAtt      >  1:386511/1‑100 (MQ=255)
                                                                                          tGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGCCTCCGAtt      >  1:697082/1‑100 (MQ=255)
                                                                                          tGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGCCTCCGAtt      >  1:567356/1‑100 (MQ=255)
                                                                                           gCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGCCTCCGATTa     >  1:20180/1‑100 (MQ=255)
                                                                                            cTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGCCTCCGATTaa    >  2:608667/1‑100 (MQ=255)
                                                                                            cTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGCCTCCGATTaa    <  1:81597/100‑1 (MQ=255)
                                                                                             tGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGCCTCCGATTaaa   <  1:1293831/100‑1 (MQ=255)
                                                                                             tGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGCCTCCGATTa     <  1:3016860/98‑1 (MQ=255)
                                                                                              ggACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGCCTCCGATTAAat  >  1:1287664/1‑100 (MQ=255)
                                                                                              ggACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGCCTCCGATTAAat  >  2:2952316/1‑100 (MQ=255)
                                                                                                   |                                                                                               
TTTATTGATAGATAAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGCCTCCGATTAAAT  >  NC_000913/3597743‑3597937

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