Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,899,995 A→G 5.4% G523G (GGA→GGG 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,899,9950AG 5.4% 395.0 / 9.1 129G523G (GGA→GGGpdeDputative c‑di‑GMP phosphodiesterase PdeD
Reads supporting (aligned to +/- strand):  ref base A (56/66);  new base G (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 = 4.49e-01

Only 100 of 130 total aligned reads displayed.

GAGGGCGTGGAAACGCAAGAACAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCC  >  NC_000913/1899846‑1900143
                                                                                                                                                     |                                                                                                                                                    
gAGGGCGTGGAAACGCAAGAACAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGGGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGa                                                                                                                                                      <  2:1425222/150‑1 (MQ=255)
 aGGGCGTGGAAACGCAAGAACAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGAc                                                                                                                                                     <  1:497508/150‑1 (MQ=255)
        ggAAACGCAAGAACAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCAc                                                                                                                                              >  2:1870690/1‑150 (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                                                                                                                                         >  1:713717/1‑150 (MQ=255)
             cGCAAGAACAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCa                                                                                                                                         <  2:1669721/150‑1 (MQ=255)
             cGCAAGAACAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCa                                                                                                                                         >  2:227173/1‑150 (MQ=255)
              gCAAGAACAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCAt                                                                                                                                        >  2:620335/1‑150 (MQ=255)
              gCAAGAACAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCAt                                                                                                                                        <  2:1405120/150‑1 (MQ=255)
               cAAGAACAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCAtt                                                                                                                                       >  1:2275370/1‑150 (MQ=255)
                 aGAACAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTAt                                                                                                                                     >  1:2038316/1‑150 (MQ=255)
                  gAACAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATg                                                                                                                                    >  1:1691464/1‑150 (MQ=255)
                     cAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGGCGTGATATGTCCATTATGCcg                                                                                                                                 <  1:1225529/150‑1 (MQ=255)
                     cAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCcg                                                                                                                                 <  1:1750023/150‑1 (MQ=255)
                      aGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCcgt                                                                                                                                <  1:2316621/150‑1 (MQ=255)
                      aGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCcgt                                                                                                                                >  1:1624874/1‑150 (MQ=255)
                       ggCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCcgtt                                                                                                                               >  1:1417671/1‑150 (MQ=255)
                       ggCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTAAGCcgtt                                                                                                                               >  2:1466001/1‑150 (MQ=255)
                       ggCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGCCATAATGGACATATCACGCCCATTATGCcgtt                                                                                                                               >  2:2451854/1‑150 (MQ=255)
                             gTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCcgttacgtta                                                                                                                         <  1:522656/150‑1 (MQ=255)
                                    cgcCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTa                                                                                                                  >  2:2418245/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)
                                            tGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCtt                                                                                                          >  2:1468501/1‑150 (MQ=255)
                                                  gCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGt                                                                                                    >  1:862361/1‑150 (MQ=255)
                                                  gCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGt                                                                                                    <  1:1475605/150‑1 (MQ=255)
                                                     tatTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTgcg                                                                                                 <  1:77099/150‑1 (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)
                                                              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)
                                                                      tGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACa                                                                                >  2:1013902/1‑150 (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:2324330/1‑150 (MQ=255)
                                                                            cacaGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGa                                                                          <  1:2061564/150‑1 (MQ=255)
                                                                            cacaGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGa                                                                          <  1:913550/150‑1 (MQ=255)
                                                                               aGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGCCGGGCAGCCAACCGCCGGCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAAcc                                                                       <  2:1324597/150‑1 (MQ=255)
                                                                                gCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCa                                                                      <  1:2305985/150‑1 (MQ=255)
                                                                                  cGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGa                                                                    >  1:2460746/1‑150 (MQ=255)
                                                                                  cGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGa                                                                    <  1:1752639/150‑1 (MQ=255)
                                                                                   gATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGAt                                                                   >  1:1540337/1‑150 (MQ=255)
                                                                                    aTGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATc                                                                  <  1:717625/150‑1 (MQ=255)
                                                                                    aTGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATc                                                                  <  1:291001/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)
                                                                                              gTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATaa                                                        >  1:1795226/1‑150 (MQ=255)
                                                                                               tGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCAGCCCGGCATAATGGACATATCCCGCCCATTATGCCGTTACGTTAACGCTAACTCATCTTCATAGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATaaa                                                       >  1:1555672/1‑149 (MQ=255)
                                                                                                 aTTTTCCCAAATGGCTTGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTTTTTTTTTTTATTTTTAAGTATTCGAACGAGATCTACACGATAATCa                                                     <  2:1608345/150‑1 (MQ=255)
                                                                                                 aTTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAGTCa                                                     >  2:2443593/1‑150 (MQ=255)
                                                                                                    ttCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATTGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTg                                                  <  1:1090784/150‑1 (MQ=255)
                                                                                                               ctggcGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAatgat                                       >  2:268235/1‑150 (MQ=255)
                                                                                                               ctggcGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTAAGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAatgat                                       <  1:1466001/150‑1 (MQ=255)
                                                                                                                tggcGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTATAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAatgatg                                      >  2:1402703/1‑150 (MQ=255)
                                                                                                                   cGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTg                                   <  2:1827939/150‑1 (MQ=255)
                                                                                                                   cGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTg                                   >  1:1097869/1‑150 (MQ=255)
                                                                                                                      gCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGata                                <  1:313961/150‑1 (MQ=255)
                                                                                                                      gCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGata                                >  1:865466/1‑150 (MQ=255)
                                                                                                                        aGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATg                              <  1:1249247/150‑1 (MQ=255)
                                                                                                                         gCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGc                             >  2:602547/1‑150 (MQ=255)
                                                                                                                         gCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGc                             >  2:971015/1‑150 (MQ=255)
                                                                                                                         gCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGc                             >  1:1433188/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)
                                                                                                                                    cccgcccgGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTAc                  <  2:2382724/150‑1 (MQ=255)
                                                                                                                                     ccgcccgGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACa                 <  1:2116139/150‑1 (MQ=255)
                                                                                                                                      cgcccgGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGATCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACat                <  2:1515430/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)
                                                                                                                                                    gACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGccc  >  2:463432/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GAGGGCGTGGAAACGCAAGAACAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCC  >  NC_000913/1899846‑1900143

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: