Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,900,005 C→T 5.5% P527S (CCC→TCC)  pdeD → putative c‑di‑GMP phosphodiesterase PdeD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,900,0050CT 5.5% 400.0 / 7.3 128P527S (CCC→TCC) pdeDputative c‑di‑GMP phosphodiesterase PdeD
Reads supporting (aligned to +/- strand):  ref base C (58/62);  new base T (3/4);  total (61/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

Only 100 of 131 total aligned reads displayed.

AAACGCAAGAACAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGA  >  NC_000913/1899856‑1900152
                                                                                                                                                     |                                                                                                                                                   
aaaCGCAAGAACAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGc                                                                                                                                                     >  1:1305744/1‑150 (MQ=255)
 aaCGCAAGAACAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGcc                                                                                                                                                    <  1:2457639/150‑1 (MQ=255)
  aCGCAAGAACAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGccc                                                                                                                                                   <  2:1627258/150‑1 (MQ=255)
   cGCAAGAACAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCa                                                                                                                                                  >  2:227173/1‑150 (MQ=255)
   cGCAAGAACAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCa                                                                                                                                                  >  1:713717/1‑150 (MQ=255)
    gCAAGAACAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCAt                                                                                                                                                 <  2:1405120/150‑1 (MQ=255)
    gCAAGAACAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCAt                                                                                                                                                 >  2:620335/1‑150 (MQ=255)
       aGAACAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTAt                                                                                                                                              >  1:2038316/1‑150 (MQ=255)
        gAACAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATg                                                                                                                                             >  1:1691464/1‑150 (MQ=255)
         aaCAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGc                                                                                                                                            <  2:1853246/150‑1 (MQ=255)
           cAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGGCGTGATATGTCCATTATGCcg                                                                                                                                          <  1:1225529/150‑1 (MQ=255)
           cAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCcg                                                                                                                                          <  1:1750023/150‑1 (MQ=255)
            aGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCcgt                                                                                                                                         >  1:1624874/1‑150 (MQ=255)
            aGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCcgt                                                                                                                                         <  1:2316621/150‑1 (MQ=255)
             ggCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCcgtt                                                                                                                                        >  1:1417671/1‑150 (MQ=255)
             ggCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGCCATAATGGACATATCACGCCCATTATGCcgtt                                                                                                                                        >  2:2451854/1‑150 (MQ=255)
                gAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCcgttacg                                                                                                                                     >  2:2451001/1‑150 (MQ=255)
                   gTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCcgttacgtta                                                                                                                                  <  1:522656/150‑1 (MQ=255)
                     aTTTGCGCAGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAAc                                                                                                                                >  2:1471995/1‑150 (MQ=255)
                          cgcCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTa                                                                                                                           >  2:2409850/1‑150 (MQ=255)
                           gcCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTAc                                                                                                                          >  1:2350140/1‑150 (MQ=255)
                            ccGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGGCGTGATATGTCCATTATGCcg                                                                                                                                          >  2:591621/1‑133 (MQ=38)
                            ccGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGGCGTGATATGTCCATTATGCcg                                                                                                                                          <  1:591621/133‑1 (MQ=38)
                             cGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTc                                                                                                                        >  1:1749114/1‑150 (MQ=255)
                              gTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTTCTCa                                                                                                                       <  2:1481845/150‑1 (MQ=255)
                               tCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCAt                                                                                                                      >  1:42204/1‑150 (MQ=255)
                                  tGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCtt                                                                                                                   >  2:1468501/1‑150 (MQ=255)
                                   ggCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTc                                                                                                                  >  2:717625/1‑150 (MQ=255)
                                       tGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCg                                                                                                              <  1:922546/150‑1 (MQ=255)
                                        gCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGt                                                                                                             <  1:1475605/150‑1 (MQ=255)
                                        gCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGt                                                                                                             >  1:862361/1‑150 (MQ=255)
                                          atatTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTgc                                                                                                           >  2:2305985/1‑150 (MQ=255)
                                             ttttGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGcc                                                                                                        <  1:231508/150‑1 (MQ=255)
                                                 gCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGtt                                                                                                    <  2:786859/150‑1 (MQ=255)
                                                  cAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCTCCGCCCGCCCGGCATAATTGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTg                                                                                                   <  1:2197768/150‑1 (MQ=255)
                                                   aaGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGc                                                                                                  <  1:1714189/150‑1 (MQ=255)
                                                   aaGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGc                                                                                                  <  1:2259176/150‑1 (MQ=255)
                                                    aGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGct                                                                                                 >  1:1131975/1‑150 (MQ=255)
                                                       gTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCtt                                                                                              <  2:2012821/150‑1 (MQ=255)
                                                          tttGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTaa                                                                                           <  2:1800661/150‑1 (MQ=255)
                                                             gtcctcacaTCCGATGCCGGTACTTGATTTTCCCACTGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGGATTATGGACATATGACGCCCATTATGCCGTTACGTTAACGGTTACTCATGTTCATCGTGCGCCGGTTGCTCTTTAACAAt                                                                                       <  1:2415941/145‑1 (MQ=255)
                                                              tACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAAt                                                                                       <  1:2349118/150‑1 (MQ=255)
                                                                cGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGc                                                                                     <  2:2038316/150‑1 (MQ=255)
                                                                cGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGc                                                                                     >  2:1626749/1‑150 (MQ=255)
                                                                cGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGc                                                                                     >  2:2324330/1‑150 (MQ=255)
                                                                 gCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGACATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCg                                                                                    >  1:2110729/1‑150 (MQ=255)
                                                                  cacaGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGa                                                                                   <  1:913550/150‑1 (MQ=255)
                                                                     aGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAAcc                                                                                <  1:1870690/150‑1 (MQ=255)
                                                                      gCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCa                                                                               <  1:2305985/150‑1 (MQ=255)
                                                                        cGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGa                                                                             <  1:1752639/150‑1 (MQ=255)
                                                                        cGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGa                                                                             >  1:2460746/1‑150 (MQ=255)
                                                                         gATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGAt                                                                            >  1:1540337/1‑150 (MQ=255)
                                                                          aTGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATc                                                                           <  1:291001/150‑1 (MQ=255)
                                                                          aTGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATc                                                                           <  1:717625/150‑1 (MQ=255)
                                                                          aTGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATc                                                                           <  2:713717/150‑1 (MQ=255)
                                                                            gCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCaa                                                                         >  2:1245607/1‑150 (MQ=255)
                                                                            gCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCaa                                                                         >  2:991570/1‑150 (MQ=255)
                                                                              cGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAaca                                                                       <  2:770881/150‑1 (MQ=255)
                                                                              cGCTACGTGATTATCCCAAATGGCTGCCCGGCAGCCAACCGCCGCCCGACCGCAATAATGCACATATCACGCCCATTATGCCCTTACCTTAACGCTTACTCATCTTCATCGTGCCCCGGTTGCTCTTTAACAATGCGAACCAGATCAaca                                                                       >  1:1070314/1‑150 (MQ=255)
                                                                               gCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAacac                                                                      <  1:470664/150‑1 (MQ=255)
                                                                                  aCGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGAt                                                                   <  1:2451854/150‑1 (MQ=255)
                                                                                   cGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATa                                                                  >  1:997854/1‑150 (MQ=255)
                                                                                     tGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCAGCCCGGCATAATGGACATATCCCGCCCATTATGCCGTTACGTTAACGCTAACTCATCTTCATAGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATaaa                                                                >  1:1555672/1‑149 (MQ=255)
                                                                                       aTTTTCCCAAATGGCTTGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTTTTTTTTTTTATTTTTAAGTATTCGAACGAGATCTACACGATAATCa                                                              <  2:1608345/150‑1 (MQ=255)
                                                                                       aTTTTCCCAAATGGCTGGCGGGCAGCCAACCTCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGCTACGTTAACGCTTACTCATCTTCATCGTGGGCCGGTTTCTCTTTAACAATTCGAACGAGATCAACACGATAATCa                                                              <  2:1244136/150‑1 (MQ=255)
                                                                                       aTTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAGTCa                                                              >  2:2443593/1‑150 (MQ=255)
                                                                                          ttCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATTGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTg                                                           <  1:1090784/150‑1 (MQ=255)
                                                                                              cAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCGCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCtt                                                       >  1:824785/1‑150 (MQ=255)
                                                                                                     ctggcGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAatgat                                                >  2:268235/1‑150 (MQ=255)
                                                                                                      tggcGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTATAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAatgatg                                               >  2:1402703/1‑150 (MQ=255)
                                                                                                         cGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTg                                            <  2:997854/150‑1 (MQ=255)
                                                                                                         cGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTg                                            >  1:1097869/1‑150 (MQ=255)
                                                                                                         cGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTg                                            <  2:1827939/150‑1 (MQ=255)
                                                                                                            gCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGata                                         <  1:313961/150‑1 (MQ=255)
                                                                                                              aGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATg                                       <  1:1249247/150‑1 (MQ=255)
                                                                                                               gCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGc                                      >  2:971015/1‑150 (MQ=255)
                                                                                                               gCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGc                                      >  1:1433188/1‑150 (MQ=255)
                                                                                                                ccAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCa                                     >  1:659807/1‑150 (MQ=255)
                                                                                                                 cAACCGCCGCCCGCCCGGCAAAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTTTTTTTTTTTATTTTTAAGAATTCGAACGAGATCTACACGATAATCATTGGCTTCAATGATGGTGATATGCAg                                    <  2:832146/150‑1 (MQ=255)
                                                                                                                  aaCCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGt                                   <  1:57705/150‑1 (MQ=255)
                                                                                                                  aaCCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGt                                   <  1:1387093/150‑1 (MQ=255)
                                                                                                                   aCCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGCTCAACACGATAATCATTGGGTTCAAAGAAGGTGATATGCAGTg                                  >  2:2351994/1‑150 (MQ=255)
                                                                                                                    ccgccgCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTAATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTgg                                 >  2:2449564/1‑150 (MQ=255)
                                                                                                                        cgcccgccCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGACCCAGATCAACACGATAATCAGTGGCTTCAATGATGGTGATATGCAGTGGCCCt                             >  1:1467224/1‑150 (MQ=255)
                                                                                                                          cccgcccgGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTAc                           <  2:2382724/150‑1 (MQ=255)
                                                                                                                           ccgcccgGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACa                          <  2:2017325/150‑1 (MQ=255)
                                                                                                                           ccgcccgGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACa                          <  1:2116139/150‑1 (MQ=255)
                                                                                                                             taacgGCATAATGGGCGTGATATGTCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACatc                        >  1:288155/4‑150 (MQ=38)
                                                                                                                             gcccgGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACatc                        >  1:46190/1‑150 (MQ=255)
                                                                                                                              cccgGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACatca                       <  1:1891018/150‑1 (MQ=255)
                                                                                                                               ccgGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCTGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACatcaa                      <  1:1464612/150‑1 (MQ=255)
                                                                                                                               ccgGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACatcaa                      <  1:1324512/150‑1 (MQ=255)
                                                                                                                                    aTAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATcaca                 <  1:1009083/150‑1 (MQ=255)
                                                                                                                                     tAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACAt                >  2:1080461/1‑150 (MQ=255)
                                                                                                                                        tGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTTCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGc             <  2:1417671/150‑1 (MQ=255)
                                                                                                                                          gACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGccc           >  2:463432/1‑150 (MQ=255)
                                                                                                                                            cATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCcac         <  2:1691464/150‑1 (MQ=255)
                                                                                                                                              tatCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCacac       >  1:479087/1‑150 (MQ=255)
                                                                                                                                                tCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATACGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACgg     >  1:837602/1‑150 (MQ=255)
                                                                                                                                                 cACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACggg    >  1:986929/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AAACGCAAGAACAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGA  >  NC_000913/1899856‑1900152

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