Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,898,316 T→C 6.2% intergenic (+20/‑111) 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,3160TC 6.2% 443.9 / 13.5 145intergenic (+20/‑111)sdaA/pdeDL‑serine deaminase I/putative c‑di‑GMP phosphodiesterase PdeD
Reads supporting (aligned to +/- strand):  ref base T (74/62);  new base C (2/7);  total (76/69)
Fisher's exact test for biased strand distribution p-value = 8.60e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.14e-01

Only 100 of 147 total aligned reads displayed.

CCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTAAAACCTATCGTC  >  NC_000913/1898167‑1898463
                                                                                                                                                     |                                                                                                                                                   
ccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCt                                                                                                                                                     >  1:373866/1‑150 (MQ=255)
ccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCt                                                                                                                                                     <  2:2142176/150‑1 (MQ=255)
 cgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTg                                                                                                                                                    <  2:2419967/150‑1 (MQ=255)
     ccAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAAc                                                                                                                                                <  1:628703/150‑1 (MQ=255)
        gTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGa                                                                                                                                             >  1:2002175/1‑150 (MQ=255)
           cACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATgg                                                                                                                                          <  1:705000/150‑1 (MQ=255)
            aCCGCGCGTCTCGCTGTATAATGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTTGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATggg                                                                                                                                         <  2:855180/150‑1 (MQ=255)
             ccGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGc                                                                                                                                        <  2:2231304/150‑1 (MQ=255)
             ccGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGc                                                                                                                                        >  2:198066/1‑150 (MQ=255)
              cgcgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCg                                                                                                                                       >  1:1137695/1‑150 (MQ=255)
              cgcgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCCGTTGCAGATGGGCg                                                                                                                                       <  2:1347241/150‑1 (MQ=255)
              cgcgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCCGTTGCAGATGGGCg                                                                                                                                       <  2:573505/150‑1 (MQ=255)
              cgcgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCCGTTGCAGATGGGCg                                                                                                                                       <  1:330015/150‑1 (MQ=255)
               gcgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCCGTTGCAGATGGGCGa                                                                                                                                      <  2:597089/150‑1 (MQ=255)
                cgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCCGTTGCAGATGGGCGag                                                                                                                                     <  1:397953/150‑2 (MQ=255)
                cgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCCGTTGCAGATGGGCGag                                                                                                                                     <  2:2335336/150‑2 (MQ=255)
                  cgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAAtt                                                                                                                                   >  1:189000/1‑150 (MQ=255)
                   gTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAAttt                                                                                                                                  <  2:2262735/150‑1 (MQ=255)
                     ctcGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTtat                                                                                                                                <  1:1802860/150‑1 (MQ=255)
                        gCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATAcc                                                                                                                             >  1:1478246/1‑150 (MQ=255)
                        gCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATAcc                                                                                                                             >  1:1734479/1‑150 (MQ=255)
                        gCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATAcc                                                                                                                             <  1:1801319/150‑1 (MQ=255)
                             aTAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCtt                                                                                                                        <  2:913751/150‑1 (MQ=255)
                              tAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCttt                                                                                                                       >  1:2168087/1‑150 (MQ=255)
                                 ggTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCATGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTtctc                                                                                                                    >  1:801891/1‑150 (MQ=255)
                                 ggTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACCTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTtctc                                                                                                                    >  2:1360757/1‑150 (MQ=255)
                                  gTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCg                                                                                                                   <  1:2290546/150‑1 (MQ=255)
                                    cATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCCCGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTc                                                                                                                 >  2:204385/1‑150 (MQ=255)
                                    cATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTc                                                                                                                 >  2:2029005/1‑150 (MQ=255)
                                        gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTctgct                                                                                                             >  1:322309/1‑150 (MQ=255)
                                         aaaCGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTctgctg                                                                                                            <  1:1563201/150‑1 (MQ=255)
                                                gTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATAtt                                                                                                     >  1:2120593/1‑150 (MQ=255)
                                                    gAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTcccc                                                                                                 >  1:366854/1‑150 (MQ=255)
                                                     aaaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCa                                                                                                <  1:2088429/150‑1 (MQ=255)
                                                     aaaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCa                                                                                                >  1:1342891/1‑150 (MQ=255)
                                                          ggTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTCCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTaca                                                                                           >  1:1496353/1‑150 (MQ=255)
                                                                   aTGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTg                                                                                  <  2:363181/150‑1 (MQ=255)
                                                                    tGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGt                                                                                 >  1:814488/1‑150 (MQ=255)
                                                                    tGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGt                                                                                 >  1:1661249/1‑150 (MQ=255)
                                                                     gAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGtt                                                                                <  2:296135/150‑1 (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)
                                                                             aaTACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAgg                                                                        <  1:888256/150‑1 (MQ=255)
                                                                              aTACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGc                                                                       <  1:1078197/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                                                                   >  2:1020263/1‑150 (MQ=255)
                                                                                   gcgAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTg                                                                  >  1:566848/1‑150 (MQ=255)
                                                                                   gcgAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTg                                                                  >  2:1825840/1‑150 (MQ=255)
                                                                                       aaCCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGc                                                              <  1:600289/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)
                                                                                                    ggtCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAg                                                 <  1:1847463/150‑1 (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:1020263/150‑1 (MQ=255)
                                                                                                           cAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGcc                                          >  1:1134583/1‑150 (MQ=255)
                                                                                                           cAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGcc                                          >  1:862433/1‑150 (MQ=255)
                                                                                                           cAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGcc                                          <  2:1226047/150‑1 (MQ=255)
                                                                                                            aaTCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCATCGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGccc                                         <  2:696283/150‑1 (MQ=255)
                                                                                                            aaTCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGccc                                         >  2:2410071/1‑150 (MQ=255)
                                                                                                             aTCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCg                                        >  2:78145/1‑150 (MQ=255)
                                                                                                                aaaGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCAt                                     >  2:2363538/1‑150 (MQ=255)
                                                                                                                 aaGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATg                                    >  2:2339904/1‑150 (MQ=255)
                                                                                                                   gTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGCCGCCATACGCTTACAGGGTGGCCCGCATGCa                                  >  2:2076083/1‑150 (MQ=255)
                                                                                                                       aGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCaaaaa                              >  1:185775/1‑150 (MQ=255)
                                                                                                                       aGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCaaaaa                              <  2:899798/150‑1 (MQ=255)
                                                                                                                        gtgtGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAg                             <  2:1758585/150‑1 (MQ=255)
                                                                                                                         tgtgACTAATACTTCTTACTCGCCCATCTGCAACGGCTGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGc                            <  2:314156/150‑1 (MQ=255)
                                                                                                                         tgtgACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGc                            <  1:2379100/150‑1 (MQ=255)
                                                                                                                            gACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGcaca                         >  2:2064675/1‑150 (MQ=255)
                                                                                                                             aCTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCAGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACaa                        >  2:1713372/1‑150 (MQ=255)
                                                                                                                                  tACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGAt                   <  2:1894439/150‑1 (MQ=255)
                                                                                                                                  tACTTCTTACTCGACCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCAACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATAAGCTTAAAGGGTGGCCCGCATGCAAAAAGCACAACGGAt                   >  2:1496103/1‑150 (MQ=255)
                                                                                                                                   aCTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATc                  <  2:2438017/150‑1 (MQ=255)
                                                                                                                                    cttcttACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCa                 <  1:1251764/150‑1 (MQ=255)
                                                                                                                                     ttcttACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCAt                <  1:1277688/150‑1 (MQ=255)
                                                                                                                                     ttcttACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCAt                >  1:975888/1‑150 (MQ=255)
                                                                                                                                      tcttACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCAtt               >  2:1378871/1‑150 (MQ=255)
                                                                                                                                      tcttACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATAATCCCAAATACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCAtt               >  1:1637694/1‑150 (MQ=255)
                                                                                                                                       cttACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTa              <  1:1621434/150‑1 (MQ=255)
                                                                                                                                          attCGCCCATCCGTTGCAGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTaaaa           >  1:1752151/3‑150 (MQ=255)
                                                                                                                                           cTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTAAAAc          >  2:131678/1‑150 (MQ=255)
                                                                                                                                           cTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTAAAAc          <  1:2410071/150‑1 (MQ=255)
                                                                                                                                            tCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTAAAAcc         <  1:1526542/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      >  1:855040/1‑150 (MQ=255)
                                                                                                                                               cccATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTAAAACCTAt      >  2:1981840/1‑150 (MQ=255)
                                                                                                                                                ccATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTAAAACCTAtc     >  1:1522291/1‑150 (MQ=255)
                                                                                                                                                 cATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTAAAACCTAtcg    <  2:1459649/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTAAAACCTATCGTC  >  NC_000913/1898167‑1898463

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: