Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 980,800 G→T 6.0% intergenic (+14/‑247) mukB → / → ldtD chromosome partitioning protein MukB/L,D‑transpeptidase LdtD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913980,8000GT 6.0% 243.9 / 3.3 73intergenic (+14/‑247)mukB/ldtDchromosome partitioning protein MukB/L,D‑transpeptidase LdtD
Reads supporting (aligned to +/- strand):  ref base G (35/33);  new base T (2/3);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGCACCACCTATAAACTGGTGCGTAAAGTCTTCCAGAATACCGAACACGTTCATGTCGTCGGCCTGCGAGGATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTAAAAACGTATCGGCGTTTATATACTGAAGATAAGCCTGATGAGTAACAGG  >  NC_000913/980651‑980949
                                                                                                                                                     |                                                                                                                                                     
aGGCACCACCTATAAACTGGTGCGTAAAGTCTTCCAGAATACCGAACACGTTCATGTCGTCGGCCTGCGAGGATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGgcag                                                                                                                                                       <  2:733887/150‑1 (MQ=255)
 ggCACCACCTATAAACTGGTGCGTAAAGTCTTCCAGAATACCGAACACGTTCATGTCGTCGGCCTGCGAGGATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGgcagc                                                                                                                                                      >  2:1098613/1‑150 (MQ=255)
    accaccTATAAACTGGTGCGTAAAGTCTTCCAGAATACCGAACACGTTCATGTCGTCGGCCTGCGAGGATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAAt                                                                                                                                                   <  1:79548/150‑1 (MQ=255)
        ccTATAAACTGGTGCGTAAAGTCTTCCAGAATACCGAACACGTTCATGTCGTCGGCCTGCGAGGATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATgccg                                                                                                                                               >  1:638091/1‑150 (MQ=255)
        ccTATAAACTGGTGCGTAAAGTCTTCCAGAATACCGAACACGTTCATGTCGTCGGCCTGCGAGGATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATgccg                                                                                                                                               <  1:1067446/150‑1 (MQ=255)
            taAACTGGTGCGTAAAGTCTTCCAGAATACCGAACACGTTCATGTCGTCGGCCTGCGAGGATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCATTGCTGCCGCctt                                                                                                                                           <  1:505951/150‑1 (MQ=255)
            taAACTGGTGCGTAAAGTCTTCCAGAATACCGAACACGTTCATGTCGTCGGCCTGCGAGGATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCATTGCTGCCGCctt                                                                                                                                           <  1:966388/150‑1 (MQ=255)
                cTGGTGCGTAAAGTCTTCCAGAATACCGAACACGTTCATGTCGTCGGCCTGCGAGGATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCcttttct                                                                                                                                       <  2:410125/150‑1 (MQ=255)
                cTGGTGCGTAAAGTCTTCCAGAATACCGAACACGTTCATGTCGTCGGCCTGCGAGGATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCcttttct                                                                                                                                       <  2:1093908/150‑1 (MQ=255)
                cTGGTGCGTAAAGTCTTCCAGAATACCGAACACGTTCATGTCGTCGGCCTGCGAGGATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCcttttct                                                                                                                                       <  2:1225327/150‑1 (MQ=255)
                   gTGCGTAAAGTCTTCCAGAATACCGAACACGTTCATGTCGTCGGCCTGCGAGGATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGAGGCAGCAATGCCGCattttctttt                                                                                                                                    >  2:1113938/1‑150 (MQ=255)
                       gTAAAGTCTTCCAGAATACCGAACACGTTCATGTCGTCGGCCTGCGAGGATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCg                                                                                                                                <  2:286341/150‑1 (MQ=255)
                          aaGTCTTCCAGAATACCGAACACGTTCATGTCGTCGGCCTGCGAGGATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGaaa                                                                                                                             >  1:832255/1‑150 (MQ=255)
                               ttCCAGAATACCGAACACGTTCATGTCGTCGGCCTGCGAGGATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGAGAGTTAAAATTAAGGCTGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTcc                                                                                                                        <  2:1799898/150‑1 (MQ=255)
                                tCCAGAATACCGAACACGTTCATGTCGTCGGCCTGCGAGGATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCg                                                                                                                       >  2:396354/1‑150 (MQ=255)
                                 ccAGAATACCGAACACGTTCATGTCGTCGGCCTGCGAGGATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGt                                                                                                                      >  1:1616751/1‑150 (MQ=255)
                                   aGAATACCGAACACGTTCATGTCGTCGGCCTGCGAGGATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGttt                                                                                                                    >  1:1171240/1‑150 (MQ=255)
                                    gAATACCGAACACGTTCATGTCGTCGGCCTGCGAGGATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTc                                                                                                                   >  2:1332927/1‑150 (MQ=255)
                                     aaTACCGAACACGTTCATGTCGTCGGCCTGCGAGGATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCt                                                                                                                  <  2:1335349/150‑1 (MQ=255)
                                         ccGAACACGTTCATGTCGTCGGCCTGCGAGGATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCAc                                                                                                              >  1:405298/1‑150 (MQ=255)
                                                    cATGTCGTCGGCCTGCGAGGATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTgg                                                                                                   <  1:763619/150‑1 (MQ=255)
                                                      tGTCGTCGGCCTGCGAGGATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGca                                                                                                 >  2:558905/1‑150 (MQ=255)
                                                            cgGCCTGCGAGGATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTAc                                                                                           >  2:1640082/1‑150 (MQ=255)
                                                               ccTGCGAGGATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGgc                                                                                        >  2:1222356/1‑150 (MQ=255)
                                                                  gCGAGGATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGgcgcg                                                                                     >  1:1497583/1‑150 (MQ=255)
                                                                      ggATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCttt                                                                                 >  1:830758/1‑150 (MQ=255)
                                                                      ggATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCttt                                                                                 <  1:1229428/150‑1 (MQ=255)
                                                                        aTTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTc                                                                               <  1:1332927/150‑1 (MQ=255)
                                                                         tttGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCt                                                                              >  2:1001600/1‑150 (MQ=255)
                                                                               ccgcAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAAc                                                                        >  1:790594/1‑150 (MQ=255)
                                                                                    aCTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAActtctt                                                                   >  1:1185069/1‑150 (MQ=255)
                                                                                     cTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCttt                                                                  >  2:1211695/1‑150 (MQ=255)
                                                                                        ccTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACa                                                               >  1:1400447/1‑150 (MQ=255)
                                                                                        ccTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACa                                                               <  1:1098613/150‑1 (MQ=255)
                                                                                          tGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACAtt                                                             >  1:535975/1‑150 (MQ=255)
                                                                                          tGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACAtt                                                             >  1:1010089/1‑150 (MQ=255)
                                                                                           gAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTa                                                            <  1:339246/150‑1 (MQ=255)
                                                                                               cGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGtt                                                        <  2:304754/150‑1 (MQ=255)
                                                                                                 cTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTAt                                                      >  1:954636/1‑150 (MQ=255)
                                                                                                 cTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTAt                                                      <  2:830758/150‑1 (MQ=255)
                                                                                                   tCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGt                                                    >  2:771557/1‑150 (MQ=255)
                                                                                                    ccccgAACTGACGCCGCGCCTTCTCAGGCGAGTTAAACTTAAGGCGGCAGCAATGCCGCCGTTTCTTTTTCCGAAACCTCCGCTTCTGCACTAAAAAAGCGGCACATTACGGCGCGCCTCTCTTTAAACTTCTTTACATTAGGTTATGTa                                                   <  2:535975/146‑1 (MQ=255)
                                                                                                    ccAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTa                                                   <  2:1933644/150‑1 (MQ=255)
                                                                                                    ccAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTa                                                   <  1:558905/150‑1 (MQ=255)
                                                                                                     cAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTaa                                                  >  2:222020/1‑150 (MQ=255)
                                                                                                     cAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTaa                                                  >  2:1483327/1‑150 (MQ=255)
                                                                                                         aaCTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTAAAAAc                                              <  2:1185069/150‑1 (MQ=255)
                                                                                                          aCTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCATCAATGCCGCCTTTTCTTTTTCCGAAAAATCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTAAAAACg                                             <  2:1884044/150‑1 (MQ=255)
                                                                                                          aCTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTAAAAACg                                             <  1:1745922/150‑1 (MQ=255)
                                                                                                             gACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTAAAAACGTAt                                          >  2:1028242/1‑150 (MQ=255)
                                                                                                              aCGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTAAAAACGTATc                                         <  2:1171240/150‑1 (MQ=255)
                                                                                                               cGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTAAAAACGTATCg                                        >  1:1444637/1‑150 (MQ=255)
                                                                                                               cGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTAAAAACGTATCg                                        >  1:104665/1‑150 (MQ=255)
                                                                                                                gAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTAAAAACGTATCgg                                       <  2:1616751/150‑1 (MQ=255)
                                                                                                                gAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTAAAAACGTATCgg                                       <  2:638091/150‑1 (MQ=255)
                                                                                                                 aaGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTAAAAACGTATCGGc                                      >  2:380908/1‑150 (MQ=255)
                                                                                                                      ccTTCTCATGCTAGTTAAAATTAAGGCGGCAGCAATGCCGACTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTAAAAACGTATCGGCGTTta                                 <  2:224697/150‑1 (MQ=255)
                                                                                                                      ccTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTAAAAACGTATCGGCGTTta                                 <  2:779583/150‑1 (MQ=255)
                                                                                                                       cTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCGTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTAAAAACGTATCGGCGTTtat                                >  2:452561/1‑150 (MQ=255)
                                                                                                                        ttCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTAAAAACGTATCGGCGTTtata                               <  2:832255/150‑1 (MQ=255)
                                                                                                                         tctcAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTAAAAACGTATCGGCGTTtatat                              >  2:1646732/1‑150 (MQ=255)
                                                                                                                                  aGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTTTGTAAAAACGTATCGGCGTTTATATACTGAAGAt                     >  2:1768506/1‑150 (MQ=255)
                                                                                                                                   gTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTAAAAACGTATCGGCGTTTATATACTGAAGATa                    <  2:1444637/150‑1 (MQ=255)
                                                                                                                                     tAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTAAAAACGTATCGGCGTTTATATACTGAAGATAAg                  >  1:23241/1‑150 (MQ=255)
                                                                                                                                       aaaTTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCCCATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTAAAAACGTATCGGCGTTTATATACTGAAGATAAGcc                >  1:1615448/1‑150 (MQ=255)
                                                                                                                                       aaaTTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTAAAAACGTATCGGCGTTTATATACTGAAGATAATcc                <  2:832889/150‑1 (MQ=255)
                                                                                                                                          aaaaGGCGGCATTGCTGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTAAAAACGTATCGGCGTTTATATACTGAAGATAAGCCtga             >  2:1808300/3‑150 (MQ=255)
                                                                                                                                           tAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTAAAAACGTATCGGCGTTTATATACTGAAGATAAGCCtgat            >  1:1022653/1‑150 (MQ=255)
                                                                                                                                           aaaGGCGGCATTGCTGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTAAAAACGTATCGGCGTTTATATACTGAAGATAAGCCtgat            >  1:1848153/2‑150 (MQ=255)
                                                                                                                                            aaGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTAAAAACGTATCGGCGTTTATATACTGAAGATAAGCCtgatg           <  2:104665/150‑1 (MQ=255)
                                                                                                                                              ggcggcAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTAAAAACGTATCGGCGTTTATATACTGAAGATAAGCCTGATGAg         <  2:1930138/150‑1 (MQ=255)
                                                                                                                                              ggcggcAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTAAAAACGTATCGGCGTTTATATACTGAAGATAAGCCTGATGAg         <  1:396354/150‑1 (MQ=255)
                                                                                                                                              ggcggcAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTAAAAACGTATCGGCGTTTATATACTGAAGATAAGCCTGATGAg         <  1:1211695/150‑1 (MQ=255)
                                                                                                                                               gcggcAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTAAAAACGTATCGGCGTTTATATACTGAAGATAAGCCTGATGAGt        <  1:1028242/150‑1 (MQ=255)
                                                                                                                                                cggcAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTAAAAACGTATCGGCGTTTATATACTGAAGATAAGCCTGATGAGTa       >  2:1837150/1‑150 (MQ=255)
                                                                                                                                                   cagcaATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTAAAAACGTATCGGCGTTTATATACTGAAGATAAGCCTGATGAGTAACa    <  2:545255/150‑1 (MQ=255)
                                                                                                                                                   cagcaATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTAAAAACGTATCGGCGTTTATATACTGAAGATAAGCCTGATGAGTAACa    >  1:334183/1‑150 (MQ=255)
                                                                                                                                                     gcaATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTAAAAACGTATCGGCGTTTATATACTg                          >  2:274185/1‑126 (MQ=255)
                                                                                                                                                     gcaATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTAAAAACGTATCGGCGTTTATATACTg                          <  1:274185/126‑1 (MQ=255)
                                                                                                                                                     gcaATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTAAAAACGTATCGGCGTTTATATACTGAAGATAAGCCTGATGAGTAACAgg  >  2:1895593/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGGCACCACCTATAAACTGGTGCGTAAAGTCTTCCAGAATACCGAACACGTTCATGTCGTCGGCCTGCGAGGATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTAAAAACGTATCGGCGTTTATATACTGAAGATAAGCCTGATGAGTAACAGG  >  NC_000913/980651‑980949

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