Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 980,801 C→T 5.5% intergenic (+15/‑246) 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,8010CT 5.5% 231.7 / 2.2 73intergenic (+15/‑246)mukB/ldtDchromosome partitioning protein MukB/L,D‑transpeptidase LdtD
Reads supporting (aligned to +/- strand):  ref base C (35/34);  new base T (2/2);  total (37/36)
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.53e-01

GGCACCACCTATAAACTGGTGCGTAAAGTCTTCCAGAATACCGAACACGTTCATGTCGTCGGCCTGCGAGGATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTAAAAACGTATCGGCGTTTATATACTGAAGATAAGCCTGATGAGTAACAGGC  >  NC_000913/980652‑980950
                                                                                                                                                     |                                                                                                                                                     
ggCACCACCTATAAACTGGTGCGTAAAGTCTTCCAGAATACCGAACACGTTCATGTCGTCGGCCTGCGAGGATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGgcagc                                                                                                                                                       >  2:1098613/1‑150 (MQ=255)
   accaccTATAAACTGGTGCGTAAAGTCTTCCAGAATACCGAACACGTTCATGTCGTCGGCCTGCGAGGATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAAt                                                                                                                                                    <  1:79548/150‑1 (MQ=255)
       ccTATAAACTGGTGCGTAAAGTCTTCCAGAATACCGAACACGTTCATGTCGTCGGCCTGCGAGGATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATgccg                                                                                                                                                <  1:1067446/150‑1 (MQ=255)
       ccTATAAACTGGTGCGTAAAGTCTTCCAGAATACCGAACACGTTCATGTCGTCGGCCTGCGAGGATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATgccg                                                                                                                                                >  1:638091/1‑150 (MQ=255)
           taAACTGGTGCGTAAAGTCTTCCAGAATACCGAACACGTTCATGTCGTCGGCCTGCGAGGATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCATTGCTGCCGCctt                                                                                                                                            <  1:966388/150‑1 (MQ=255)
           taAACTGGTGCGTAAAGTCTTCCAGAATACCGAACACGTTCATGTCGTCGGCCTGCGAGGATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCATTGCTGCCGCctt                                                                                                                                            <  1:505951/150‑1 (MQ=255)
               cTGGTGCGTAAAGTCTTCCAGAATACCGAACACGTTCATGTCGTCGGCCTGCGAGGATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCcttttct                                                                                                                                        <  2:1225327/150‑1 (MQ=255)
               cTGGTGCGTAAAGTCTTCCAGAATACCGAACACGTTCATGTCGTCGGCCTGCGAGGATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCcttttct                                                                                                                                        <  2:410125/150‑1 (MQ=255)
               cTGGTGCGTAAAGTCTTCCAGAATACCGAACACGTTCATGTCGTCGGCCTGCGAGGATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCcttttct                                                                                                                                        <  2:1093908/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:1229428/150‑1 (MQ=255)
                                                                     ggATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCttt                                                                                  >  1:830758/1‑150 (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:1010089/1‑150 (MQ=255)
                                                                                         tGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACAtt                                                              >  1:535975/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                                                    <  1:558905/150‑1 (MQ=255)
                                                                                                   ccAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTa                                                    <  2:1933644/150‑1 (MQ=255)
                                                                                                    cAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTaa                                                   >  2:1483327/1‑150 (MQ=255)
                                                                                                    cAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTaa                                                   >  2:222020/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:104665/1‑150 (MQ=255)
                                                                                                              cGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTAAAAACGTATCg                                         >  1:1444637/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:1211695/150‑1 (MQ=255)
                                                                                                                                             ggcggcAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTAAAAACGTATCGGCGTTTATATACTGAAGATAAGCCTGATGAg          <  1:396354/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)
                                                                                                                                                     caATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTAAAAACGTATCGGCGTTTATATACTGAAGATAAGCCTGATGAGTAACAGGc  <  1:1001600/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCACCACCTATAAACTGGTGCGTAAAGTCTTCCAGAATACCGAACACGTTCATGTCGTCGGCCTGCGAGGATTTGCGCCGCAACTCCCTGAAACGCTTCCAGGAACTGACGAAGCGCCTTCTCAGGCGAGTTAAAATTAAGGCGGCAGCAATGCCGCCTTTTCTTTTTCCGAAAACTCCGTTTCTGCACTAAAAAAGTGGCACATTACGGCGCGCTTTTCTTTAAACTTCTTTACATTAGGTTATGTAAAAACGTATCGGCGTTTATATACTGAAGATAAGCCTGATGAGTAACAGGC  >  NC_000913/980652‑980950

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: