Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,898,319 A→T 6.1% intergenic (+23/‑108) sdaA → / → pdeD L‑serine deaminase I/putative c‑di‑GMP phosphodiesterase PdeD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,898,3190AT 6.1% 472.3 / 14.3 148intergenic (+23/‑108)sdaA/pdeDL‑serine deaminase I/putative c‑di‑GMP phosphodiesterase PdeD
Reads supporting (aligned to +/- strand):  ref base A (76/63);  new base T (2/7);  total (78/70)
Fisher's exact test for biased strand distribution p-value = 8.50e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.55e-01

Only 100 of 150 total aligned reads displayed.

CCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTAAAACCTATCGTCGTAAT  >  NC_000913/1898172‑1898468
                                                                                                                                                   |                                                                                                                                                     
 cAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACg                                                                                                                                                    >  1:296135/1‑150 (MQ=255)
     gCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATg                                                                                                                                                >  2:2261771/1‑150 (MQ=255)
      cACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATgg                                                                                                                                               <  1:705000/150‑1 (MQ=255)
       aCCGCGCGTCTCGCTGTATAATGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTTGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATggg                                                                                                                                              <  2:855180/150‑1 (MQ=255)
        ccGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGc                                                                                                                                             <  2:2351649/150‑1 (MQ=255)
        ccGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGc                                                                                                                                             <  2:2231304/150‑1 (MQ=255)
         cgcgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCCGTTGCAGATGGGCg                                                                                                                                            <  1:330015/150‑1 (MQ=255)
         cgcgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCCGTTGCAGATGGGCg                                                                                                                                            <  2:573505/150‑1 (MQ=255)
         cgcgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCCGTTGCAGATGGGCg                                                                                                                                            <  2:1347241/150‑1 (MQ=255)
          gcgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCCGTTGCAGATGGGCGa                                                                                                                                           <  2:597089/150‑1 (MQ=255)
           cgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCCGTTGCAGATGGGCGag                                                                                                                                          <  2:2335336/150‑2 (MQ=255)
             cgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAAtt                                                                                                                                        <  1:1871120/150‑1 (MQ=255)
                ctcGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTtat                                                                                                                                     <  1:1802860/150‑1 (MQ=255)
                   gCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATAcc                                                                                                                                  <  1:1178861/150‑1 (MQ=255)
                        aTAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCtt                                                                                                                             <  2:913751/150‑1 (MQ=255)
                            ggTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTtctc                                                                                                                         >  1:1807143/1‑150 (MQ=255)
                            ggTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACCTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTtctc                                                                                                                         >  2:1360757/1‑150 (MQ=255)
                              tCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGt                                                                                                                       >  1:1133142/1‑150 (MQ=255)
                               cATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCCCGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTc                                                                                                                      >  2:204385/1‑150 (MQ=255)
                                aTCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTct                                                                                                                     >  1:175323/1‑150 (MQ=255)
                                  cGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTctgc                                                                                                                   <  2:547933/150‑1 (MQ=255)
                                   gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTctgct                                                                                                                  >  1:2271066/1‑150 (MQ=255)
                                       cGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTaa                                                                                                              <  2:1986949/150‑1 (MQ=255)
                                           gTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATAtt                                                                                                          >  1:2120593/1‑150 (MQ=255)
                                               gAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTcccc                                                                                                      >  1:366854/1‑150 (MQ=255)
                                                aaaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCa                                                                                                     >  1:1342891/1‑150 (MQ=255)
                                                  aCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACt                                                                                                   >  1:2153514/1‑150 (MQ=255)
                                                     ggTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTCCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTaca                                                                                                >  1:1496353/1‑150 (MQ=255)
                                                               tGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGt                                                                                      >  1:814488/1‑150 (MQ=255)
                                                               tGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGt                                                                                      >  2:1951066/1‑150 (MQ=255)
                                                               tGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGt                                                                                      >  1:1661249/1‑150 (MQ=255)
                                                                  aCGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGc                                                                                   >  1:218892/1‑150 (MQ=255)
                                                                   cGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCg                                                                                  <  2:1782623/150‑1 (MQ=255)
                                                                      cAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCa                                                                               <  2:175323/150‑1 (MQ=255)
                                                                      cAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACAGTTGCGTCa                                                                               >  1:2033767/1‑150 (MQ=255)
                                                                        aaTACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAgg                                                                             <  1:261288/150‑1 (MQ=255)
                                                                         aTACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGc                                                                            <  1:1977962/150‑1 (MQ=255)
                                                                          tACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCg                                                                           >  1:497194/1‑150 (MQ=255)
                                                                            ccGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGtt                                                                         >  2:1522722/1‑150 (MQ=255)
                                                                            ccGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGtt                                                                         <  2:1233073/150‑1 (MQ=255)
                                                                             cgcgAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGttt                                                                        >  1:547743/1‑150 (MQ=255)
                                                                             cgcgAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGttt                                                                        >  2:1020263/1‑150 (MQ=255)
                                                                                  aaCCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGc                                                                   <  1:600289/150‑1 (MQ=255)
                                                                                   aCCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGcc                                                                  <  1:601972/150‑1 (MQ=255)
                                                                                    ccTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCa                                                                 >  2:2438653/1‑150 (MQ=255)
                                                                                     cTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCAt                                                                <  2:1647749/150‑1 (MQ=255)
                                                                                            ggtggtCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTa                                                         >  1:2292214/1‑150 (MQ=255)
                                                                                              tggtCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACa                                                       >  2:1290144/1‑150 (MQ=255)
                                                                                                 tCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAggg                                                    >  1:1057497/1‑150 (MQ=255)
                                                                                                   tGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTg                                                  >  2:1846354/1‑150 (MQ=255)
                                                                                                    ggCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTgg                                                 >  2:534254/1‑150 (MQ=255)
                                                                                                      cAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGcc                                               >  1:862433/1‑150 (MQ=255)
                                                                                                      cAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGcc                                               <  2:1226047/150‑1 (MQ=255)
                                                                                                      cAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGcc                                               >  1:853812/1‑150 (MQ=255)
                                                                                                      cAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGcc                                               <  1:1020263/150‑1 (MQ=255)
                                                                                                      cAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGcc                                               >  1:1134583/1‑150 (MQ=255)
                                                                                                       aaTCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCATCGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGccc                                              <  2:696283/150‑1 (MQ=255)
                                                                                                       aaTCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGccc                                              >  2:2410071/1‑150 (MQ=255)
                                                                                                        aTCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCg                                             >  2:1830168/1‑150 (MQ=255)
                                                                                                        aTCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCg                                             >  2:78145/1‑150 (MQ=255)
                                                                                                        aTCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCg                                             >  1:870809/1‑150 (MQ=255)
                                                                                                         tCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGc                                            >  2:2045027/1‑150 (MQ=255)
                                                                                                          cAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCa                                           >  2:2436381/1‑150 (MQ=255)
                                                                                                             aGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGc                                        <  1:78145/150‑1 (MQ=255)
                                                                                                              gTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGCCGCCATACGCTTACAGGGTGGCCCGCATGCa                                       >  2:2076083/1‑150 (MQ=255)
                                                                                                                ccAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTCTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCaaa                                     >  1:1643626/1‑150 (MQ=255)
                                                                                                                  aGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCaaaaa                                   <  2:899798/150‑1 (MQ=255)
                                                                                                                  aGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCaaaaa                                   >  1:185775/1‑150 (MQ=255)
                                                                                                                   gtgtGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAg                                  <  2:1758585/150‑1 (MQ=255)
                                                                                                                   gtgtGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAg                                  >  1:1296933/1‑150 (MQ=255)
                                                                                                                    tgtgACTAATACTTCTTACTCGCCCATCTGCAACGGCTGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGc                                 <  2:314156/150‑1 (MQ=255)
                                                                                                                    tgtgACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGc                                 <  1:1566166/150‑1 (MQ=255)
                                                                                                                        aCTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCAGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACaa                             >  2:1713372/1‑150 (MQ=255)
                                                                                                                             tACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGAt                        <  2:1894439/150‑1 (MQ=255)
                                                                                                                             tACTTCTTACTCGACCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCAACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATAAGCTTAAAGGGTGGCCCGCATGCAAAAAGCACAACGGAt                        >  2:1496103/1‑150 (MQ=255)
                                                                                                                              aCTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGCCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATc                       <  1:2076083/150‑1 (MQ=255)
                                                                                                                               cttcttACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCa                      >  2:292961/1‑150 (MQ=255)
                                                                                                                               cttcttACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCa                      <  1:2363538/150‑1 (MQ=255)
                                                                                                                               cttcttACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCa                      <  1:1251764/150‑1 (MQ=255)
                                                                                                                               cttcttACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCa                      >  2:2314872/1‑150 (MQ=255)
                                                                                                                               cttcttACTCGCCCATCTGAAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCATTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCa                      >  1:1492172/1‑150 (MQ=255)
                                                                                                                                ttcttACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCAt                     >  1:975888/1‑150 (MQ=255)
                                                                                                                                 tcttACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCAtt                    >  2:1378871/1‑150 (MQ=255)
                                                                                                                                  cttACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTa                   <  1:1621434/150‑1 (MQ=255)
                                                                                                                                     attCGCCCATCCGTTGCAGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTaaaa                >  1:1752151/3‑150 (MQ=255)
                                                                                                                                      ttCGCCCATCCGTTGCAGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTAAAAc               >  1:379220/2‑150 (MQ=255)
                                                                                                                                      cTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTAAAAc               >  2:131678/1‑150 (MQ=255)
                                                                                                                                      cTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTAAAAc               >  1:1525557/1‑150 (MQ=255)
                                                                                                                                      cTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTAAAAc               <  1:2410071/150‑1 (MQ=255)
                                                                                                                                       tCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTAAAAcc              <  1:1526542/150‑1 (MQ=255)
                                                                                                                                       tCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTAAAAcc              <  2:322309/150‑1 (MQ=255)
                                                                                                                                        cGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTAAAAccg             >  2:675319/1‑149 (MQ=255)
                                                                                                                                        cGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTAAAACCt             >  2:1776457/1‑150 (MQ=255)
                                                                                                                                        cGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATACTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTAAAACCt             >  2:22903/1‑150 (MQ=255)
                                                                                                                                          cccATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTAAAACCTAt           >  2:1981840/1‑150 (MQ=255)
                                                                                                                                            cATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTAAAACCTAtcg         <  2:1459649/150‑1 (MQ=255)
                                                                                                                                              tCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTAAAACCTAtcgtc       <  1:198066/150‑1 (MQ=255)
                                                                                                                                                 gCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTAAAACCTATCGTCGTa    >  1:2428813/1‑150 (MQ=255)
                                                                                                                                                 gCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTAAAACCTATCGTCGTa    >  2:2100720/1‑150 (MQ=255)
                                                                                                                                                   aaCGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTAAAACCTATCGTCGTAAt  >  1:330320/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTAAAACCTATCGTCGTAAT  >  NC_000913/1898172‑1898468

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