Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 568,685 T→A 16.4% intergenic (+38/‑217) emrE → / → ybcK DLP12 prophage; multidrug resistance protein/DLP12 prophage; putative recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913568,6850TA16.4% 388.8 / 61.4 152intergenic (+38/‑217)emrE/ybcKDLP12 prophage; multidrug resistance protein/DLP12 prophage; putative recombinase
Reads supporting (aligned to +/- strand):  ref base T (66/61);  new base A (15/10);  total (81/71)
Fisher's exact test for biased strand distribution p-value = 5.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

Only 100 of 183 total aligned reads displayed.

TGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATA  >  NC_000913/568586‑568784
                                                                                                   |                                                                                                   
tgatgTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGctct                                                                                                     <  1:747014/100‑1 (MQ=255)
   tgTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCttat                                                                                                  >  1:3581025/1‑100 (MQ=255)
   tgTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCttac                                                                                                  <  2:1092265/100‑2 (MQ=255)
    gTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCttata                                                                                                 >  2:3158501/1‑100 (MQ=255)
     ttGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCttatat                                                                                                >  2:2572364/1‑100 (MQ=255)
      tGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCttatatt                                                                                               >  2:1837381/1‑100 (MQ=255)
      tGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCttatatt                                                                                               <  2:2700074/100‑1 (MQ=255)
       gATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCttatattt                                                                                              >  2:624174/1‑100 (MQ=255)
        aTTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGCAGGCTCttatattta                                                                                             >  2:4205735/1‑100 (MQ=255)
           tgtgCTGATGTGTTGATCATTAACTTATCGTCACGAAGCACACCACATTAACATTATTTGTTTCTCATCGACGAAACTATGGAGGCTCttatatttatat                                                                                          <  1:2830739/100‑1 (MQ=255)
            gtgCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATg                                                                                         <  1:2572357/100‑1 (MQ=255)
            gtgCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATg                                                                                         >  2:436015/1‑100 (MQ=255)
            gtgCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATg                                                                                         >  2:369159/1‑100 (MQ=255)
              gCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAg                                                                                       <  1:1478450/100‑1 (MQ=255)
               ccGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGc                                                                                      >  1:744667/1‑100 (MQ=255)
                cGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCATATAAATATAAGAGcc                                                                                     <  1:2525637/100‑1 (MQ=255)
                cGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCATATAAATATAAGAGcc                                                                                     <  1:2525663/100‑1 (MQ=255)
                cGGTGTGTTGATTATTAATATATTATCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGcc                                                                                     >  1:3729789/1‑100 (MQ=255)
                  gtgtgtTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCATATAAATATAAGAGCCTc                                                                                   <  2:1921699/100‑1 (MQ=38)
                   tgtgtTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCATATAAATATAAGAGCCTcc                                                                                  <  1:3884728/100‑2 (MQ=38)
                     tgtTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGtt                                                                                >  1:4206728/1‑100 (MQ=255)
                     tgtTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGt                                                                                 >  1:1849753/1‑99 (MQ=255)
                      gtTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGttt                                                                               <  1:2696143/100‑1 (MQ=255)
                      gtTGATAAATAATTTATTATCACGAAAAACACAAAATTAAAATAATTTGATACTAAACGACTAAAATATGGAATCTCATATAAATATAAGAGccaccat                                                                                <  1:1944051/99‑6 (MQ=255)
                        tGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTa                                                                             <  2:1036163/100‑1 (MQ=255)
                        tGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTa                                                                             <  2:1126058/100‑1 (MQ=255)
                        tGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTa                                                                             >  2:788850/1‑100 (MQ=255)
                         gATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTAt                                                                            <  2:1857260/100‑1 (MQ=255)
                           ttattaATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGc                                                                          <  2:2426130/100‑1 (MQ=255)
                            tattaATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGTCTCATATAAATATAAGAGCCTccatattttag                                                                         <  2:4051036/100‑11 (MQ=25)
                            tattaATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCt                                                                         >  1:3401260/1‑100 (MQ=255)
                            tattaATTTATTGTCAAGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAGAATATGGAGGCTCATATAAATATAAGAGCCTccatattttag                                                                         <  1:2554828/100‑11 (MQ=21)
                              ttaATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCttt                                                                       >  1:1966351/1‑100 (MQ=255)
                               taATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCtttt                                                                      <  1:116791/100‑1 (MQ=255)
                                aaTTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCttttt                                                                     <  1:1217662/100‑1 (MQ=255)
                                  tttATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGt                                                                   <  2:646792/100‑1 (MQ=255)
                                   ttATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGtt                                                                  <  1:2746453/100‑1 (MQ=255)
                                    tATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTa                                                                 >  1:876027/1‑100 (MQ=255)
                                     aTTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTaa                                                                >  1:2735885/1‑100 (MQ=255)
                                      ttGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAAt                                                               <  2:2519109/100‑1 (MQ=255)
                                       tGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATg                                                              <  1:3111587/100‑1 (MQ=255)
                                         tCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTCTTTGTTAATGTc                                                            <  2:2801897/100‑1 (MQ=255)
                                          cACGAAGCACACCACATTAAAATAATTTGTATCGAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTAGTTAATGTCt                                                           >  1:2810688/1‑100 (MQ=255)
                                            cGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCttt                                                         >  2:215037/1‑100 (MQ=255)
                                            cGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCttt                                                         >  2:1086739/1‑100 (MQ=255)
                                                 cacaCCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTAtttt                                                    <  2:2833842/100‑1 (MQ=255)
                                                 cacaCCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTAtttt                                                    <  2:2833715/100‑1 (MQ=255)
                                                   caccacATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTAtttttt                                                  <  1:4233060/100‑1 (MQ=255)
                                                   caccacATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTAtttttt                                                  <  2:4140279/100‑1 (MQ=255)
                                                     ccacATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTAt                                                >  2:1570318/1‑100 (MQ=255)
                                                     ccacATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTAt                                                <  1:1604143/100‑1 (MQ=255)
                                                      cacaTTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTAt                                                <  1:709152/99‑1 (MQ=255)
                                                        caTTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTa                                             >  1:3531296/1‑100 (MQ=255)
                                                           tAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTc                                          >  1:1036510/1‑100 (MQ=255)
                                                            aaaaTAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCt                                         >  1:2304126/1‑100 (MQ=255)
                                                             aaaTAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCtt                                        >  2:332597/1‑100 (MQ=255)
                                                                taatTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTtg                                     <  2:3185439/100‑1 (MQ=255)
                                                                taatTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTtg                                     <  1:3915159/100‑1 (MQ=255)
                                                                  atTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTtgtg                                   >  2:3932528/1‑100 (MQ=255)
                                                                  atTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTtgtg                                   <  2:3451378/100‑1 (MQ=255)
                                                                      gTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCttt                               >  2:81391/1‑100 (MQ=255)
                                                                        ttCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTc                              >  2:388217/1‑99 (MQ=255)
                                                                        ttCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCa                             <  1:640858/100‑1 (MQ=255)
                                                                          cTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAg                           <  2:1293274/100‑1 (MQ=255)
                                                                            aaaCGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGAt                         >  2:875153/1‑100 (MQ=255)
                                                                            aaaCGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGAt                         >  1:4163896/1‑100 (MQ=255)
                                                                             aaCGACTAAATTATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGAtt                        <  1:1719221/100‑1 (MQ=255)
                                                                             aaCGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGAtt                        >  2:379739/1‑100 (MQ=255)
                                                                              aCGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTa                       <  2:2459746/100‑1 (MQ=255)
                                                                               cGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTAt                      >  1:470891/1‑100 (MQ=255)
                                                                                 aCTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATgg                    <  1:4126472/100‑1 (MQ=255)
                                                                                 aCTAAAAGATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATgg                    <  2:362485/100‑1 (MQ=255)
                                                                                  cTAAAATATGGAGGCTCTTATATTTATCTGCGCCTCGTTTTATGCTTTTTGTTAATGTATTTATTTTTTATGTATTCTTTTGAGCTTTGAAGATTATg                     >  1:1677261/1‑98 (MQ=255)
                                                                                  cTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGc                   >  2:1984654/1‑100 (MQ=255)
                                                                                  cTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGc                   >  2:1003279/1‑100 (MQ=255)
                                                                                  cTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGc                   >  2:3581148/1‑100 (MQ=255)
                                                                                   gcataaaacgAGGCTCATATAAATATAAGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAATATTATGTCg                  >  1:2051994/10‑100 (MQ=21)
                                                                                    cataaaacgAGGCTCATATAAATATAAGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGt                 >  2:1655990/9‑100 (MQ=25)
                                                                                    cataaaacgAGGCTCATATAAATATAAGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGt                 >  1:4094025/9‑100 (MQ=25)
                                                                                    aaaaTATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGt                 <  2:470893/100‑1 (MQ=255)
                                                                                    aaaaTATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTATTTTGTGCTTTCAAGATTATGGCGt                 >  2:535006/1‑100 (MQ=255)
                                                                                     ataaaacgAGGCTCATATAAATATAAGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTa                >  2:2935011/8‑100 (MQ=25)
                                                                                      taaaacgAGGCTCATATAAATATAAGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTaa               >  2:4013217/7‑100 (MQ=37)
                                                                                       aaaacgAGGCTCATATAAATATAAGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTTCTTTCAAGATTATGTCGt                 >  1:2769822/6‑97 (MQ=25)
                                                                                       aaaacgAGGCTCATATAAATATAAGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTTAAGATTATGGCGTAAg              >  1:31769/6‑100 (MQ=37)
                                                                                        tatGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGa             >  1:2981964/1‑100 (MQ=255)
                                                                                        aaacgAGGCTCATATAAATATAAGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGa             >  2:2100869/5‑100 (MQ=37)
                                                                                        aaacgAGGCTCATATAAATATAAGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGa             >  2:3273829/5‑100 (MQ=37)
                                                                                         atGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGaa            >  1:2933029/1‑100 (MQ=255)
                                                                                         aacgAGGCTCATATAAATATAAGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATTTCGTAAGa             >  1:2157193/4‑99 (MQ=37)
                                                                                         aacgAGGCTCATATAAATATAAGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGaa            >  1:2405626/4‑100 (MQ=38)
                                                                                          acgAGGCTCATATAAATATAAGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTTCAGATTATGGCGTTAGa             >  1:1011939/3‑98 (MQ=37)
                                                                                           ggAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGaaaa          >  1:3028615/1‑100 (MQ=255)
                                                                                           ggAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGaaaa          >  1:1435231/1‑100 (MQ=255)
                                                                                              ggCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTg       >  2:3802936/1‑100 (MQ=255)
                                                                                               gCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGc      <  1:1542986/100‑1 (MQ=255)
                                                                                                ctctTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAATATTATGGCGTAAGAAAATTGCa     >  2:2985500/1‑100 (MQ=255)
                                                                                                ctctTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCa     <  2:4091258/100‑1 (MQ=255)
                                                                                                ctctTATATTTATATGAGCATCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGc      <  1:4144640/99‑1 (MQ=255)
                                                                                                  ctTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAAt   >  2:761190/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATA  >  NC_000913/568586‑568784

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