Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,900,000 T→A 5.0% I525N (ATC→AAC) ‡ 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,0000TA 5.0% 416.4 / 7.1 129I525N (ATC→AAC) ‡pdeDputative c‑di‑GMP phosphodiesterase PdeD
Reads supporting (aligned to +/- strand):  ref base T (56/66);  new base A (3/4);  total (59/70)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 131 total aligned reads displayed.

GGAAACGCAAGAACAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACGG  >  NC_000913/1899854‑1900149
                                                                                                                                                  |                                                                                                                                                     
ggAAACGCAAGAACAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCAc                                                                                                                                                    <  2:580748/150‑1 (MQ=255)
 gAAACGCAAGAACAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAACTGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACg                                                                                                                                                   <  1:2204442/150‑1 (MQ=255)
  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)
       cAAGAACAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCAtt                                                                                                                                             >  1:2275370/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)
             cAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGGCGTGATATGTCCATTATGCcg                                                                                                                                       <  1:2479410/150‑1 (MQ=255)
             cAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCcg                                                                                                                                       <  1:1750023/150‑1 (MQ=255)
              aGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCcgt                                                                                                                                      >  1:1624874/1‑150 (MQ=255)
               ggCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCcgtt                                                                                                                                     >  1:1417671/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:2418245/1‑150 (MQ=255)
                             gcCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTAc                                                                                                                       >  1:2350140/1‑150 (MQ=255)
                               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)
                                         tGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCg                                                                                                           <  1:922546/150‑1 (MQ=255)
                                          gCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGt                                                                                                          <  1:1475605/150‑1 (MQ=255)
                                            atatTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTgc                                                                                                        >  2:2305985/1‑150 (MQ=255)
                                            atatTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTgc                                                                                                        <  1:620335/150‑1 (MQ=255)
                                             tatTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTgcg                                                                                                       <  1:77099/150‑1 (MQ=255)
                                               ttttGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGcc                                                                                                     <  1:231508/150‑1 (MQ=255)
                                                   gCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGtt                                                                                                 >  1:2330517/1‑150 (MQ=255)
                                                   gCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGtt                                                                                                 <  2:786859/150‑1 (MQ=255)
                                                    cAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCTCCGCCCGCCCGGCATAATTGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTg                                                                                                <  1:2197768/150‑1 (MQ=255)
                                                     aaGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGc                                                                                               <  1:1714189/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)
                                                             ttGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGTATTTATTTTTATGATTATGCGTTTACGTTTATGTTTACTGATCTTCATTGTTTGCCTTTTGGTCTTTatc                                                                                       <  2:721556/150‑3 (MQ=255)
                                                             ttGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAAc                                                                                       <  2:950582/150‑1 (MQ=255)
                                                              tGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACa                                                                                      <  1:420891/150‑1 (MQ=255)
                                                                tACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAAt                                                                                    <  1:2349118/150‑1 (MQ=255)
                                                                  cGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGc                                                                                  >  2:1626749/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)
                                                                         ccGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAg                                                                           <  1:2426244/150‑1 (MQ=255)
                                                                          cGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGa                                                                          <  1:1752639/150‑1 (MQ=255)
                                                                           gATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGAt                                                                         >  1:1540337/1‑150 (MQ=255)
                                                                            aTGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATc                                                                        <  2:713717/150‑1 (MQ=255)
                                                                            aTGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATc                                                                        <  1:291001/150‑1 (MQ=255)
                                                                              gCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCaa                                                                      >  2:991570/1‑150 (MQ=255)
                                                                              gCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCaa                                                                      >  2:1245607/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                                         >  1:1097869/1‑150 (MQ=255)
                                                                                                           cGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTg                                         <  2:1827939/150‑1 (MQ=255)
                                                                                                              gCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGata                                      >  1:865466/1‑150 (MQ=255)
                                                                                                                aGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATg                                    <  1:1249247/150‑1 (MQ=255)
                                                                                                                 gCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGc                                   >  1:1433188/1‑150 (MQ=255)
                                                                                                                 gCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGc                                   >  2:602547/1‑150 (MQ=255)
                                                                                                                 gCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGc                                   >  2:971015/1‑150 (MQ=255)
                                                                                                                  ccAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCa                                  >  1:659807/1‑150 (MQ=255)
                                                                                                                    aaCCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGt                                <  1:1387093/150‑1 (MQ=255)
                                                                                                                    aaCCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGt                                <  1:57705/150‑1 (MQ=255)
                                                                                                                     aCCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGCTCAACACGATAATCATTGGGTTCAAAGAAGGTGATATGCAGTg                               >  2:2351994/1‑150 (MQ=255)
                                                                                                                      ccgccgCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTAATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTgg                              >  2:2449564/1‑150 (MQ=255)
                                                                                                                         ccgcccgcCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGccc                           <  1:519300/150‑1 (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)
                                                                                                                               gcccgGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACatc                     <  2:2350140/150‑1 (MQ=255)
                                                                                                                               gcccgGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAAAGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACatc                     >  2:1957297/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)
                                                                                                                                  cgGCATAATGGGCGTGATATGTCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACAAGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACatcaat                  >  2:2238536/1‑150 (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)
                                                                                                                                                  tCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATACGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACgg  >  1:837602/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GGAAACGCAAGAACAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACGG  >  NC_000913/1899854‑1900149

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: ‑