Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,165,628 A→T 5.2% intergenic (+107/‑40) yegP → / → yegQ DUF1508 domain‑containing protein YegP/putative peptidase YegQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,165,6280AT 5.2% 249.2 / 2.9 77intergenic (+107/‑40)yegP/yegQDUF1508 domain‑containing protein YegP/putative peptidase YegQ
Reads supporting (aligned to +/- strand):  ref base A (41/32);  new base T (2/2);  total (43/34)
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.15e-01

TAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCGGGCCAGCCGCGTTATTCCCT  >  NC_000913/2165479‑2165774
                                                                                                                                                     |                                                                                                                                                  
aaaaGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGGTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTtcatca                                                                                                                                                    <  1:1848128/149‑1 (MQ=255)
         tGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCg                                                                                                                                           <  2:382899/150‑1 (MQ=255)
            cacaAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGAt                                                                                                                                        <  1:202533/150‑1 (MQ=255)
             acaAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATc                                                                                                                                       <  1:894216/150‑1 (MQ=255)
             acaAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCACCTGATGAAGCggtta                                                                                                                                       <  2:1245017/150‑4 (MQ=255)
               aaGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTg                                                                                                                                     <  1:1664747/150‑1 (MQ=255)
                aGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGa                                                                                                                                    <  2:765661/150‑1 (MQ=255)
                aGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGa                                                                                                                                    >  1:1006821/1‑150 (MQ=255)
                       cGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGtca                                                                                                                             >  1:512113/1‑150 (MQ=255)
                        ggTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGtcat                                                                                                                            <  1:1691440/150‑1 (MQ=255)
                         gTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGtcatc                                                                                                                           >  1:1433810/1‑150 (MQ=255)
                          tGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGtcatca                                                                                                                          >  2:1795278/1‑150 (MQ=255)
                                aCAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACga                                                                                                                    >  2:205996/1‑150 (MQ=255)
                                   aTACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACgagag                                                                                                                 >  1:1748305/1‑150 (MQ=255)
                                     aCGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGaa                                                                                                               <  1:1098744/150‑1 (MQ=255)
                                         aaTTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTAt                                                                                                           >  1:1272108/1‑150 (MQ=255)
                                             aCGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGttt                                                                                                       >  2:565860/1‑150 (MQ=255)
                                               gATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTaa                                                                                                     >  2:1335087/1‑150 (MQ=255)
                                                  gCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAAc                                                                                                   <  2:1244235/149‑1 (MQ=255)
                                                  gCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAAc                                                                                                   >  1:1244235/1‑149 (MQ=255)
                                                   ccGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCg                                                                                                 >  2:590324/1‑150 (MQ=255)
                                                       gCGTGGGTTGTTGCCGTCCGGCGCTTTAACAGAGCCCACCACAGCGCTACAATGCCCGACCTTAAAGTGGTGGCACTCCCCTAACCGCTTCATCAGTTGAAGAGGATCTGACCTGTCATAAGAACGAGAGAATTATCTTTAAACCGGAAc                                                                                             >  2:381140/1‑150 (MQ=255)
                                                         gTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCCCAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAACCCGGACCTc                                                                                           >  1:1433881/1‑150 (MQ=255)
                                                           gggTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCt                                                                                         >  1:1421459/1‑150 (MQ=255)
                                                                gttgGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTccc                                                                                    >  1:638318/1‑150 (MQ=255)
                                                                 ttgGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTcccc                                                                                   >  1:1693537/1‑150 (MQ=255)
                                                                      gTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCcggcgg                                                                              <  1:1795278/150‑1 (MQ=255)
                                                                      gTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCcggcgg                                                                              <  1:1513993/150‑1 (MQ=255)
                                                                      gTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCcggcgg                                                                              <  2:1369776/150‑1 (MQ=255)
                                                                       tCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCggg                                                                             <  2:196423/150‑1 (MQ=255)
                                                                       tCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCggg                                                                             <  2:1256368/150‑1 (MQ=255)
                                                                           gcgcTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACg                                                                         >  2:869013/1‑150 (MQ=255)
                                                                            cgcTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGc                                                                        <  1:590324/150‑1 (MQ=255)
                                                                                 tACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGaaa                                                                   >  2:226989/1‑150 (MQ=255)
                                                                                 tACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGaaa                                                                   <  1:1634534/150‑1 (MQ=255)
                                                                                  aCCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGaaaa                                                                  >  1:1049407/1‑150 (MQ=255)
                                                                                       agCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATg                                                             <  1:1590821/150‑1 (MQ=255)
                                                                                        gCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGc                                                            >  1:622374/1‑150 (MQ=255)
                                                                                        gCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGc                                                            <  2:1006821/150‑1 (MQ=255)
                                                                                          cccccACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGt                                                          <  1:268263/150‑1 (MQ=255)
                                                                                            cccACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTa                                                        >  1:360476/1‑150 (MQ=255)
                                                                                             ccACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTAc                                                       >  1:1964144/1‑150 (MQ=255)
                                                                                                caGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCt                                                    >  1:78221/1‑150 (MQ=255)
                                                                                                   cgcTACAATGCCCGCCCTTAAAGTGGCGGGACTCCCCTAACCGCGTCATCAGGTGAAGCGGATGTGACGTGTCATCAGAACGAGAGATTTATGTTTAAACCTGAACTCCTTTCCCCGGCGTGAACGCTGAAAAATATGCGTTACGCTTTc                                                 <  2:479616/150‑1 (MQ=255)
                                                                                                    aggggagtGCCCCCACTTTAAGGGCGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCg                                                >  2:140931/8‑150 (MQ=255)
                                                                                                     cTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCgg                                                                                                >  1:555553/1‑101 (MQ=255)
                                                                                                     cTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCgg                                                                                                <  2:555553/101‑1 (MQ=255)
                                                                                                     cTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGc                                               >  1:221542/1‑150 (MQ=255)
                                                                                                      tACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCt                                              <  1:461956/150‑1 (MQ=255)
                                                                                                       aCAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCta                                             >  2:140719/1‑149 (MQ=255)
                                                                                                          aTGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATg                                          <  1:1335087/150‑1 (MQ=255)
                                                                                                          aTGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATg                                          <  1:205996/150‑1 (MQ=255)
                                                                                                            gcccccacTTTAAGGGCGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGgc                                        >  1:1927546/1‑150 (MQ=255)
                                                                                                            acccgcccTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGgc                                        >  1:1757559/2‑150 (MQ=255)
                                                                                                               cgcccTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCa                                     >  2:141741/1‑150 (MQ=255)
                                                                                                                 cccTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGa                                   >  2:336965/1‑150 (MQ=255)
                                                                                                                   cTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATg                                 <  2:1748305/150‑1 (MQ=255)
                                                                                                                   cTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATg                                 >  2:1879658/1‑150 (MQ=255)
                                                                                                                   cTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATg                                 <  2:1693537/150‑1 (MQ=255)
                                                                                                                   cTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATg                                 >  2:1640652/1‑150 (MQ=255)
                                                                                                                    ttAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGc                                <  1:377772/150‑1 (MQ=255)
                                                                                                                         gTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGttt                           >  2:302221/1‑150 (MQ=255)
                                                                                                                          tGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTa                          <  1:1239963/150‑1 (MQ=255)
                                                                                                                               gCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCgg                     <  2:512113/150‑1 (MQ=255)
                                                                                                                               gCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCgg                     <  1:1640652/150‑1 (MQ=255)
                                                                                                                                   tCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCGGGCCa                 >  2:1827144/1‑150 (MQ=255)
                                                                                                                                    ccccTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCGGGCCAg                >  2:245666/1‑150 (MQ=255)
                                                                                                                                    ccccTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTACCCGGCGGGAACGCTGATAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCGGGCCAg                >  2:505442/1‑150 (MQ=255)
                                                                                                                                     cccTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCGGGCCAGc               <  1:1638747/150‑1 (MQ=255)
                                                                                                                                        gatccGCTTCACCTGATGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATgg                                         >  2:1107458/4‑121 (MQ=255)
                                                                                                                                        gatccGCTTCACCTGATGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATgg                                         <  1:1107458/118‑1 (MQ=255)
                                                                                                                                         aaCCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCGGGCCAGCcgcg           <  2:1757559/150‑1 (MQ=255)
                                                                                                                                         aaCCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCGGGCCAGCcgcg           >  2:1414306/1‑150 (MQ=255)
                                                                                                                                         aaCCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCGGGCCAGCcgcg           >  1:73302/1‑150 (MQ=255)
                                                                                                                                          tccGCTTCACCTGATGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCGGGCCAGCCGCGt          >  2:93436/2‑150 (MQ=255)
                                                                                                                                            cGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCGGGCCAGCCGCGTTa        >  2:1799732/1‑150 (MQ=255)
                                                                                                                                            cGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCGGGCCAGCCGCGTTa        <  2:1421459/150‑1 (MQ=255)
                                                                                                                                             gCTTCATCAGGTGAAGCGGATCTGAACAGTCAACAGAACGAAAAAATTATGATTTAAACGTAAAAACTTTCCACGCCGTGAAATATTAAAAATATGCGTTACGCTTTGGCTTATGGCGCTGATGCTGTTTATGCGGGCCAGCCGCGTTAt       >  2:1839097/1‑150 (MQ=255)
                                                                                                                                                tcatcaGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCGGGCCAGCCGCGTTATTcc    <  1:1799732/150‑1 (MQ=255)
                                                                                                                                                 catcaGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCGGGCCAGCCGCGTTATTccc   >  2:1956439/1‑150 (MQ=255)
                                                                                                                                                 catcaGGTGAAGCGGATCTGACCTCTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCGGGCCAGCCGCGTTATTcc  >  2:830878/1‑149 (MQ=255)
                                                                                                                                                  atcaGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCGGGCCAGCCGCGTTATTCCCt  <  1:465075/150‑1 (MQ=255)
                                                                                                                                                  atcaGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCGGGCCAGCCGCGTTATTCCCt  >  1:1105514/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCGGGCCAGCCGCGTTATTCCCT  >  NC_000913/2165479‑2165774

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: