Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,900,003 C→T 6.1% T526M (ACG→ATG)  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,0030CT 6.1% 413.0 / 10.7 131T526M (ACG→ATG) pdeDputative c‑di‑GMP phosphodiesterase PdeD
Reads supporting (aligned to +/- strand):  ref base C (58/65);  new base T (3/5);  total (61/70)
Fisher's exact test for biased strand distribution p-value = 7.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.35e-01

Only 100 of 134 total aligned reads displayed.

GGAAACGCAAGAACAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGA  >  NC_000913/1899854‑1900152
                                                                                                                                                     |                                                                                                                                                     
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                                                                                                                                                  >  1:1265390/1‑150 (MQ=255)
     cGCAAGAACAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCa                                                                                                                                                  >  2:227173/1‑150 (MQ=255)
     cGCAAGAACAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCa                                                                                                                                                  <  2:1669721/150‑1 (MQ=255)
      gCAAGAACAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCAt                                                                                                                                                 <  2:1405120/150‑1 (MQ=255)
         aGAACAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTAt                                                                                                                                              >  1:2038316/1‑150 (MQ=255)
           aaCAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGc                                                                                                                                            <  1:1117398/150‑1 (MQ=255)
           aaCAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGc                                                                                                                                            <  2:1853246/150‑1 (MQ=255)
             cAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGGCGTGATATGTCCATTATGCcg                                                                                                                                          <  1:1225529/150‑1 (MQ=255)
              aGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCcgt                                                                                                                                         >  1:1624874/1‑150 (MQ=255)
              aGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCcgt                                                                                                                                         <  1:2316621/150‑1 (MQ=255)
               ggCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGGCGTGATATGTCCATTATGCcggg                                                                                                                                        <  1:1414277/150‑3 (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)
                       aTTTGCGCAGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAAc                                                                                                                                >  2:1471995/1‑150 (MQ=255)
                            cgcCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTa                                                                                                                           >  2:2418245/1‑150 (MQ=255)
                             gcCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTAc                                                                                                                          >  1:2350140/1‑150 (MQ=255)
                              ccGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGGCGTGATATGTCCATTATGCcg                                                                                                                                          >  2:591621/1‑133 (MQ=38)
                               cGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTc                                                                                                                        >  1:1749114/1‑150 (MQ=255)
                                gTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTTCTCa                                                                                                                       <  2:1481845/150‑1 (MQ=255)
                                 tCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCAt                                                                                                                      >  1:42204/1‑150 (MQ=255)
                                     ggCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTc                                                                                                                  >  2:717625/1‑150 (MQ=255)
                                          gCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGt                                                                                                             >  1:862361/1‑150 (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:1619311/150‑1 (MQ=255)
                                                   gCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGtt                                                                                                    <  2:786859/150‑1 (MQ=255)
                                                    cAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCTCCGCCCGCCCGGCATAATTGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTg                                                                                                   <  1:2197768/150‑1 (MQ=255)
                                                     aaGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGc                                                                                                  <  1:2259176/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                                                                                         >  2:1013902/1‑150 (MQ=255)
                                                               gtcctcacaTCCGATGCCGGTACTTGATTTTCCCACTGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGGATTATGGACATATGACGCCCATTATGCCGTTACGTTAACGGTTACTCATGTTCATCGTGCGCCGGTTGCTCTTTAACAAt                                                                                       <  1:2415941/145‑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)
                                                                    cacaGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGa                                                                                   >  1:1742085/1‑150 (MQ=255)
                                                                    cacaGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGa                                                                                   <  1:2061564/150‑1 (MQ=255)
                                                                       aGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAAcc                                                                                <  1:1870690/150‑1 (MQ=255)
                                                                       aGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGCCGGGCAGCCAACCGCCGGCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAAcc                                                                                <  2:1324597/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)
                                                                            aTGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATc                                                                           <  1:717625/150‑1 (MQ=255)
                                                                            aTGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATc                                                                           <  2:713717/150‑1 (MQ=255)
                                                                            aTGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATc                                                                           <  1:291001/150‑1 (MQ=255)
                                                                              gCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCaa                                                                         >  2:1245607/1‑150 (MQ=255)
                                                                              gCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCaa                                                                         >  2:991570/1‑150 (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)
                                                                                         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)
                                                                                                       ctggcGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTAAGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAatgat                                                <  1:1466001/150‑1 (MQ=255)
                                                                                                           cGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTg                                            <  2:997854/150‑1 (MQ=255)
                                                                                                           cGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTg                                            <  2:1827939/150‑1 (MQ=255)
                                                                                                              gCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGata                                         <  1:313961/150‑1 (MQ=255)
                                                                                                              gCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGata                                         >  1:865466/1‑150 (MQ=255)
                                                                                                                 gCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGc                                      >  1:1433188/1‑150 (MQ=255)
                                                                                                                 gCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGc                                      >  2:602547/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)
                                                                                                                     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)
                                                                                                                              cgcccgGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGATCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACat                         <  2:1515430/150‑1 (MQ=255)
                                                                                                                               taacgGCATAATGGGCGTGATATGTCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACatc                        >  1:288155/4‑150 (MQ=38)
                                                                                                                               gcccgGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACatc                        <  2:2350140/150‑1 (MQ=255)
                                                                                                                               gcccgGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACatc                        >  1:46190/1‑150 (MQ=255)
                                                                                                                               gcccgGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAAAGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACatc                        >  2:1957297/1‑150 (MQ=255)
                                                                                                                                 ccgGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACatcaa                      <  1:1324512/150‑1 (MQ=255)
                                                                                                                                  cgGCATAATGGGCGTGATATGTCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACAAGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACatcaat                     >  2:2238536/1‑150 (MQ=255)
                                                                                                                                      aTAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATcaca                 <  1:1009083/150‑1 (MQ=255)
                                                                                                                                          tGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTTCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGc             <  2:1417671/150‑1 (MQ=255)
                                                                                                                                                   cACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACggg    >  1:986929/1‑150 (MQ=255)
                                                                                                                                                    aCGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACgggg   <  2:2460746/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGAAACGCAAGAACAGGCGAAGTATTTGCGCCGTCATGGCGTGCATATTTTGCAAGGGTATTTGTACGCACAGCCGATGCCGCTACGTGATTTTCCCAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTTACGTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGA  >  NC_000913/1899854‑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: