Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,162,147 T→A 9.6% intergenic (‑47/+30) PP_3663 ← / ← pssA GGDEF domain‑containing protein/CDP‑diacylglycerol‑‑serine O‑phosphatidyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,162,1470TA 9.6% 367.7 / 27.4 127intergenic (‑47/+30)PP_3663/pssAGGDEF domain‑containing protein/CDP‑diacylglycerol‑‑serine O‑phosphatidyltransferase
Reads supporting (aligned to +/- strand):  ref base T (57/56);  new base A (6/6);  total (63/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.27e-01

Only 100 of 128 total aligned reads displayed.

CCCGACTGTATGCGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCCTCTGCCGAACCGATTCTCGG  >  NC_002947/4162000‑4162296
                                                                                                                                                   |                                                                                                                                                     
 ccGACTGTATGCGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGc                                                                                                                                                    >  2:2485733/1‑150 (MQ=255)
   gACTGTATGCGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCg                                                                                                                                                  >  2:1825562/1‑150 (MQ=255)
           gCGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGACCCTGTCGCCGGCAAGCCAGCGACAgg                                                                                                                                          <  2:210885/150‑1 (MQ=255)
           gCGATCGTCAACACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGCCATCCACCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAgg                                                                                                                                          <  1:2970878/150‑1 (MQ=255)
             gATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGt                                                                                                                                        >  2:1971936/1‑150 (MQ=255)
             gATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGt                                                                                                                                        <  1:2104974/150‑1 (MQ=255)
               tCGTGATCACGTTCAGAGGGGCCTTGAACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCa                                                                                                                                      >  1:2147206/1‑150 (MQ=255)
                cGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCCGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAg                                                                                                                                     <  2:2943775/150‑1 (MQ=255)
                 gTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCTAGGCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTTTGGCCCTGTCGCTGGATTGACGGCGACAGGGTCAGt                                                                                                                                    >  2:1482140/1‑150 (MQ=255)
                   gATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTAc                                                                                                                                  >  1:451698/1‑150 (MQ=255)
                      cACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAgg                                                                                                                               <  1:2393311/150‑1 (MQ=255)
                       aCGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGt                                                                                                                                    >  2:566604/1‑144 (MQ=255)
                       aCGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGt                                                                                                                                    <  1:566604/144‑1 (MQ=255)
                       aCGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGc                                                                                                                              <  1:1319078/150‑1 (MQ=255)
                         gTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCTATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGc                                                                                                                                                    <  2:480902/126‑1 (MQ=255)
                         gTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCTATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGc                                                                                                                                                    >  1:480902/1‑126 (MQ=255)
                         gTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGa                                                                                                                            >  2:1136246/1‑150 (MQ=255)
                         gTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGa                                                                                                                            <  2:2314581/150‑1 (MQ=255)
                         gTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCAGGCCACAGGGTCAGTACAGGCGa                                                                                                                            >  1:158661/1‑150 (MQ=255)
                          ttCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGAt                                                                                                                           <  1:1992613/150‑1 (MQ=255)
                          ttCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGAt                                                                                                                           <  1:283171/150‑1 (MQ=255)
                           tCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGACCCTGTCGCCGGCAAGCCAGCGACAGGGCCa                                                                                                                                      <  1:2484714/138‑1 (MQ=38)
                                ggggCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCaga                                                                                                                     <  1:1042555/150‑1 (MQ=255)
                                ggggCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCaga                                                                                                                     >  1:2088717/1‑150 (MQ=255)
                                    ccTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCgg                                                                                                                                                 >  2:1760188/1‑118 (MQ=255)
                                    ccTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCgg                                                                                                                                                 <  1:1760188/118‑1 (MQ=255)
                                      ttGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCAccc                                                                                                               >  2:395529/1‑150 (MQ=255)
                                      ggtcACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCAccc                                                                                                               <  2:2451091/147‑1 (MQ=255)
                                        gCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGAt                                                                                                                           >  1:2616297/1‑136 (MQ=255)
                                        gCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGAt                                                                                                                           <  2:2616297/136‑1 (MQ=255)
                                        gCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGc                                                                                                             <  2:3056320/150‑1 (MQ=255)
                                        gCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGc                                                                                                             >  1:152673/1‑150 (MQ=255)
                                        gCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCCGGCAAGCCAGCGACa                                                                                                                                            >  1:1125370/1‑119 (MQ=255)
                                        gCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCCGGCAAGCCAGCGACa                                                                                                                                            <  2:1125370/119‑1 (MQ=255)
                                          aCTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCtt                                                                                                           <  2:251728/150‑1 (MQ=255)
                                             gTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCAACACTCTGGACTCGTCGCAGGCTTCTCGCTGACTGTCTCAGTACAGGCTATGTCAGGGCATCCGCTTGAt                                                                                                        <  2:31311/150‑1 (MQ=255)
                                              ttGCCTGCACATTGTAATCCAGCCCCTGACAAAAGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATc                                                                                                       <  1:1482140/150‑1 (MQ=255)
                                                 ccTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAAc                                                                                                    <  2:152673/150‑1 (MQ=255)
                                                 ccTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAAc                                                                                                    >  2:2608625/1‑150 (MQ=255)
                                                   tGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACgg                                                                                                  <  1:306561/150‑1 (MQ=255)
                                                    gCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGc                                                                                                 >  2:88925/1‑150 (MQ=255)
                                                        aTTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCa                                                                                             >  1:2707310/1‑150 (MQ=255)
                                                            tAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCg                                                                                         >  1:2629632/1‑150 (MQ=255)
                                                             aaTCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCgg                                                                                        <  2:2629632/150‑1 (MQ=255)
                                                             aaTCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCgg                                                                                        >  1:548878/1‑150 (MQ=255)
                                                               tCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCtt                                                                                                           >  1:922389/1‑129 (MQ=255)
                                                               tCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCtt                                                                                                           <  2:922389/129‑1 (MQ=255)
                                                                    ccccTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCagag                                                                                                                    <  1:121612/115‑1 (MQ=255)
                                                                    ccccTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCagag                                                                                                                    >  2:121612/1‑115 (MQ=255)
                                                                      ccTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAgcgc                                                                               >  1:2389724/1‑150 (MQ=255)
                                                                          aCAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGAt                                                                                                        <  1:1524843/121‑1 (MQ=255)
                                                                            aaaGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGa                                                                         >  1:975603/1‑150 (MQ=255)
                                                                              aGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACGGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACg                                                                       >  1:2000791/1‑150 (MQ=255)
                                                                              aGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACg                                                                       >  2:3005869/1‑150 (MQ=255)
                                                                               gggCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGc                                                                      <  2:2389724/150‑1 (MQ=255)
                                                                                  ccGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCCCTAAGATGCTGCAGGCGCGGTGACCTTGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGCCGCAGa                                                                   <  1:1136246/150‑1 (MQ=255)
                                                                                      gggCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGAt                                                                                                        <  1:502659/109‑1 (MQ=255)
                                                                                      gggCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGAt                                                                                                        >  2:502659/1‑109 (MQ=255)
                                                                                               tatCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAAc                                                                <  1:1435864/140‑1 (MQ=255)
                                                                                               tatCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAAc                                                                >  2:1435864/1‑140 (MQ=255)
                                                                                                 tCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCActtgc                                                    <  2:2147206/150‑1 (MQ=255)
                                                                                                  cATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCATAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCActtgct                                                   >  1:500025/1‑150 (MQ=255)
                                                                                                   aTAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGGGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCActtgctt                                                  >  1:2907272/1‑150 (MQ=255)
                                                                                                   aTAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCaa                                                                             >  2:800179/1‑123 (MQ=255)
                                                                                                   aTAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCaa                                                                             <  1:800179/123‑1 (MQ=255)
                                                                                                   aTAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCActtgctt                                                  <  1:3039380/150‑1 (MQ=255)
                                                                                                   aTAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCCGGCAAGCCAGCGACAGGGCCAgag                                                                                                                                   <  2:1376246/69‑3 (MQ=17)
                                                                                                   aTAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCCGGCAAGCCAGCGACAGGGCCAgag                                                                                                                                   >  1:1376246/1‑67 (MQ=17)
                                                                                                    tAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTcc                                                                                              >  1:2875444/1‑105 (MQ=255)
                                                                                                    tAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTcc                                                                                              <  2:2875444/105‑1 (MQ=255)
                                                                                                       tGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGAGAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGa                                              >  1:873139/1‑150 (MQ=255)
                                                                                                       tGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCaat                                                                                                    <  2:770836/96‑2 (MQ=255)
                                                                                                       tGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCaat                                                                                                    <  2:771926/96‑2 (MQ=255)
                                                                                                       tGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCaat                                                                                                    >  1:771926/1‑95 (MQ=255)
                                                                                                       tGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGa                                              <  1:1825562/150‑1 (MQ=255)
                                                                                                           aTCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCActtgctt                                                  <  2:2265733/142‑1 (MQ=255)
                                                                                                            tCAATGCGCGGTGCGCACCACTCTGGCCCTGTGGATGAGGTCAAGTGGACCGACGAACTAGAGGCTGTGGCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGAGAAACCGACGCAGAAACTTGCGCACTTGCTTTTGTTGTTc                                         <  2:2968494/150‑1 (MQ=255)
                                                                                                            tCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGGTTGCCGCCGAGATGGTCAGTACATGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGAAACCGACGCAGAAACTTGCGCACTTGCTTGTGATGTTc                                         <  2:2968100/150‑1 (MQ=255)
                                                                                                            tCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTgcgc                                                          >  1:181710/1‑133 (MQ=255)
                                                                                                             cAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAgc                                                                                 >  2:155476/1‑109 (MQ=255)
                                                                                                                        gcgcACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTcccc                             <  2:1902525/150‑1 (MQ=255)
                                                                                                                            accacTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGc                                                                                                 >  1:2532471/1‑78 (MQ=255)
                                                                                                                            accacTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGc                                                                                                 <  2:2532471/78‑1 (MQ=255)
                                                                                                                            accacTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGa                                              <  2:558806/129‑1 (MQ=255)
                                                                                                                            accacTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGa                                              >  1:558806/1‑129 (MQ=255)
                                                                                                                              cacTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCa                                                                                                                       >  1:816242/1‑54 (MQ=255)
                                                                                                                              cacTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAgcgc                       >  2:700937/1‑150 (MQ=255)
                                                                                                                              cacTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAgcgc                       >  1:1102320/1‑150 (MQ=255)
                                                                                                                              cacTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAgagc                       >  2:489947/1‑150 (MQ=255)
                                                                                                                                     acccTGTCGCCGGCAAGCCAGCGACAGGGTCAGTACAGGCAATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCGGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCCTCTGcc                >  1:2933500/2‑150 (MQ=255)
                                                                                                                                      cccTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCGGCTTGATCAACGGCACCCACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTGGCTTGGGATGTTCGACCGCCTCCCCCAGCGCctatttcg               >  2:545630/1‑144 (MQ=255)
                                                                                                                                       ccTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGc                                                                      >  1:897093/1‑94 (MQ=255)
                                                                                                                                       ccTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGc                                                                      <  2:897093/94‑1 (MQ=255)
                                                                                                                                       ccTGTCGCCGGCAAGCCAGCGAGAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAGCCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCCTCTGCCGa              >  2:1089706/1‑150 (MQ=255)
                                                                                                                                         tGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCCTCTGCCGAAc            <  2:1102320/150‑1 (MQ=255)
                                                                                                                                         tGTCGCCGGCAAGCCAGCGACAGGGCCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAAc                                                                                                    >  2:214158/1‑62 (MQ=18)
                                                                                                                                         tGTCGCCGGCAAGCCAGCGACAGGGCCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAAc                                                                                                    <  1:214158/62‑1 (MQ=18)
                                                                                                                                          gTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCCTCTGCCGAAcc           >  1:1958006/1‑150 (MQ=255)
                                                                                                                                                 gCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCCTCTGCCGAACCGATtctc    >  2:470843/1‑150 (MQ=255)
                                                                                                                                                   ttGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCCTCTGCCGAACCGATTCTCgg  <  2:1839136/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CCCGACTGTATGCGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCCTCTGCCGAACCGATTCTCGG  >  NC_002947/4162000‑4162296

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