Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,516,177 C→A 23.3% D177Y (GAT→TAT)  sgcR ← putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,516,1770CA23.3% 725.9 / 170.4 283D177Y (GAT→TAT) sgcRputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (108/109);  new base A (26/40);  total (134/149)
Fisher's exact test for biased strand distribution p-value = 1.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 284 total aligned reads displayed.

GTGTTGCTACAGCATGCGGCTCAACGGCATCAAATTTTGAGTGATCGGCCATGAGAATAACTTCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGTCCTCCGATAATGTAACTTTCCACATGGGGAAAACAGCCCAGTTCATTGGCGATTTTAA  >  CP009273/4516030‑4516326
                                                                                                                                                   |                                                                                                                                                     
  gtTGCTACAGCATGCGGCTCAACGGCATCAAATTTTGAGTGATCGGCCATGAGAATAACTTCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAg                                                                                                                                                   <  2:506032/150‑1 (MQ=255)
    tGCTACAGCATGCGGCTCAACGGCATCAAATTTTGAGTGATCGGCCATGAGAATAACTTCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGga                                                                                                                                                 >  1:2996278/1‑150 (MQ=255)
       tACAGCATGCGGCTCAACGGCATCAAATTTTGAGTGATCGGCCATGAGAATAACTTCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGAt                                                                                                                                              >  2:2114490/1‑150 (MQ=255)
        aCAGCATGCGGCTCAACGGCATCAAATTTTGAGTGATCGGCCATGAGAATAACTTCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATa                                                                                                                                             <  1:2896022/150‑1 (MQ=255)
        aCAGCATGCGGCTCAACGGCATCAAATTTTGAGTGATCGGCCATGAGAATAACTTCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATa                                                                                                                                             >  1:1001748/1‑150 (MQ=255)
             aTGCGGCTCAACGGCATCAAATTTTGAGTGATCGGCCATGAGAATAACTTCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGc                                                                                                                                        >  2:2487958/1‑150 (MQ=255)
             aTGCGGCTCAACGGCATCAAATTTTGAGTGATCGGCCATGAGAATAACTTCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGc                                                                                                                                        >  2:196147/1‑150 (MQ=255)
               gCGGCTCAACGGCATCAAATTTTGAGTGATCGGCCATGAGAATAACTTCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCAc                                                                                                                                      <  1:3250482/150‑1 (MQ=255)
               gCGGCTCAACGGCATCAAATTTTGAGTGATCGGCCATGAGAATAACTTCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCAc                                                                                                                                      <  2:13766/150‑1 (MQ=255)
               gCGGCTCAACGGCATCAAATTTTGAGTGATCGGCCATGAGAATAACTTCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATAGCAGGAGATAAACGCAc                                                                                                                                      <  2:1747337/150‑1 (MQ=255)
                 ggCTCAACGGCATCAAATTTTGAGTGATCGGCCATGAGAATAACTTCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGt                                                                                                                                    <  2:1484463/150‑1 (MQ=255)
                 ggCTCAACGGCATCAAATTTTGAGTGATCGGCCATGAGAATAACTTCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGt                                                                                                                                    <  1:3459999/150‑1 (MQ=255)
                 ggCTCAACGGCATCAAATTTTGAGTGATCGGCCATGAGAATAACTTCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGt                                                                                                                                    >  2:3286891/1‑150 (MQ=255)
                 ggCTCAACGGCATCAAATTTTGAGTGATCGGCCATGAGAATAACTTCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGt                                                                                                                                    >  2:3307830/1‑150 (MQ=255)
                  gCTCAACGGCATCAAATTTTGAGTGATCGGCCATGAGAATAACTTCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGtt                                                                                                                                   <  2:581878/150‑1 (MQ=255)
                    tCAACGGCATCAAATTTTGAGTGATCGGCCATGAGAATAACTTCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCt                                                                                                                                 <  2:3532624/150‑1 (MQ=255)
                     cAACGGCATCAAATTTTGAGTGATCGGCCATGAGAATAACTTCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTa                                                                                                                                <  1:2692932/150‑1 (MQ=255)
                      aaCGGCATCAAATTTTGAGTGATCGGCCATGAGAATAACTTCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATAGCAGGAGATAAACGCACGTTCTAc                                                                                                                               >  1:2874219/1‑150 (MQ=255)
                       aCGGCATCAAATTTTGAGTGATCGGCCATGAGAATAACTTCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATAGCAGGAGATAAACGCACGTTCTACg                                                                                                                              <  2:3095416/150‑1 (MQ=255)
                         ggCATCAAATTTTGAGTGATCGGCCATGAGAATAACTTCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGa                                                                                                                            >  2:3705613/1‑150 (MQ=255)
                         ggCATCAAATTTTGAGTGATCGGCCATGAGAATAACTTCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGa                                                                                                                            >  2:1275262/1‑150 (MQ=255)
                          gCATCAAATTTTGAGTGATCGGCCATGAGAATAACTTCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATAGCAGGAGATAAACGCACGTTCTACGGaa                                                                                                                           >  2:1803380/1‑150 (MQ=255)
                           cATCAAATTTTGAGTGATCGGCCATGAGAATAACTTCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGaaa                                                                                                                          >  2:924976/1‑150 (MQ=255)
                              cAAATTTTGAGTGATCGGCCATGAGAATAACTTCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAgg                                                                                                                       >  2:1172078/1‑150 (MQ=255)
                              cAAATTTTGAGTGATCGGCCATGAGAATAACTTCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAgg                                                                                                                       <  1:3020752/150‑1 (MQ=255)
                               aaaTTTTGAGTGATCGGCCATGAGAATAACTTCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGc                                                                                                                      <  1:689186/150‑1 (MQ=255)
                               aaaTTTTGAGTGATCGGCCATGAGAATAACTTCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGc                                                                                                                      <  1:502033/150‑1 (MQ=255)
                               aaaTTTTGAGTGATCGGCCATGAGAATAACTTCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGc                                                                                                                      <  1:1791235/150‑1 (MQ=255)
                                    ttGAGTGATCGGCCATGAGAATAACTTCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGa                                                                                                                 >  1:777859/1‑150 (MQ=255)
                                     tGAGTGATCGGCCATGAGAATAACTTCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGat                                                                                                                <  1:916836/150‑1 (MQ=255)
                                      gAGTGATCGGCCATGAGAATAACTTCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATAGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGatc                                                                                                               <  2:2963271/150‑1 (MQ=255)
                                       aGTGATCGGCCATGAGAATAACTTCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGatca                                                                                                              <  2:917150/150‑1 (MQ=255)
                                       aGTGATCGGCCATGAGAATAACTTCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATAGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGatca                                                                                                              <  2:1895610/150‑1 (MQ=255)
                                             cGGCCATGAGAATAACTTCACGCGAGCGTTGATTAATCCGTGTTTTCACCCCGACCTCATACTTGGTGGCATTGGTGATCGCGGTCTCCAGCGACTGCGCATCGCAGGAGATAATCGGACGTTCTACGGAAAAGGCATTGATCATCTCCa                                                                                                        <  1:4127711/150‑1 (MQ=255)
                                             cGGCCATGAGAATAACTTCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCa                                                                                                        <  1:1875778/150‑1 (MQ=255)
                                              ggCCATGAGAATAACTTCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAg                                                                                                       >  1:3197402/1‑150 (MQ=255)
                                                          aaCTTCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGActctc                                                                                           >  1:3544192/1‑150 (MQ=255)
                                                           aCTTCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATAGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGActctct                                                                                          <  2:3791483/150‑1 (MQ=255)
                                                           aCTTCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATAGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGActctct                                                                                          >  1:253517/1‑150 (MQ=255)
                                                              tCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCt                                                                                       <  1:1066407/150‑1 (MQ=255)
                                                              tCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCt                                                                                       <  2:4112328/150‑1 (MQ=255)
                                                                aCGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTAc                                                                                     >  1:3840460/1‑150 (MQ=255)
                                                                 cgcgAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTAcc                                                                                    >  2:3498809/1‑150 (MQ=255)
                                                                  gcgAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCg                                                                                   <  2:1005785/150‑1 (MQ=255)
                                                                          tGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATaa                                                                           <  1:115772/150‑1 (MQ=255)
                                                                             ataatCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCt                                                                        >  2:2542851/1‑150 (MQ=255)
                                                                             ataatCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATAGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCt                                                                        <  1:2736839/150‑1 (MQ=255)
                                                                             ataatCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATAGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCt                                                                        <  2:3981109/150‑1 (MQ=255)
                                                                             ataatCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATAGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCt                                                                        >  1:3443203/1‑150 (MQ=255)
                                                                              taatCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTg                                                                       <  1:1852674/150‑1 (MQ=255)
                                                                                atCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTggg                                                                     <  1:3427061/150‑1 (MQ=255)
                                                                                atCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTggg                                                                     <  1:1493600/150‑1 (MQ=255)
                                                                                 tCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGc                                                                    >  2:3836245/1‑150 (MQ=255)
                                                                                  ccGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCg                                                                   >  2:2761296/1‑150 (MQ=255)
                                                                                      gtTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATc                                                               <  1:3327196/150‑1 (MQ=255)
                                                                                      gtTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATc                                                               >  2:471065/1‑150 (MQ=255)
                                                                                      gtTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATAGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATc                                                               <  1:2891264/150‑1 (MQ=255)
                                                                                       ttttCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATAGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCa                                                              <  1:3079588/150‑1 (MQ=255)
                                                                                         ttCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGt                                                            >  2:1182900/1‑150 (MQ=255)
                                                                                           cACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATAGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGtcc                                                          >  1:1367646/1‑150 (MQ=255)
                                                                                            aCCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGtcct                                                         <  1:1334975/150‑1 (MQ=255)
                                                                                            aCCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATAGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGtcct                                                         <  2:3025997/150‑1 (MQ=255)
                                                                                            aCCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATAGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGtcct                                                         <  1:2938341/150‑1 (MQ=255)
                                                                                            aCCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATAGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGtcct                                                         <  1:3608311/150‑1 (MQ=255)
                                                                                               ccGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGTCCTCCg                                                      >  2:3922143/1‑150 (MQ=255)
                                                                                                 gACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGTCCTCCGAt                                                    <  2:824252/150‑1 (MQ=255)
                                                                                                 gACCTCAAACATGGTGGCATTGGTGATCCCCGTCACCAGCTACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGAACATCTCCAGCTCCAGACTCTCACCTACCGAAAAATAACCTGGGCGAATCAGACCTCCGAt                                                    >  2:2015945/1‑150 (MQ=255)
                                                                                                  aCCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGTCCTCCGATa                                                   <  2:2242915/150‑1 (MQ=255)
                                                                                                    cTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGTCCTCCGATAAt                                                 >  1:571090/1‑150 (MQ=255)
                                                                                                    cTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATAGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGTCCTCCGATAAt                                                 <  1:1039234/150‑1 (MQ=255)
                                                                                                      cAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGTCCTCCGATAATGt                                               >  1:1634945/1‑150 (MQ=255)
                                                                                                      cAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATAGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGTCCTCCGATAATGt                                               >  2:1107352/1‑150 (MQ=255)
                                                                                                      aaaaCATGGTGGCATTGGTGAACCCCGTCTCCAGCGACAGCGCATCGAAGGAGATAAAAGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGTCCTCCGATAATGt                                               >  2:2268679/2‑150 (MQ=255)
                                                                                                       aaaCATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATAGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGTCCTCCGATAATGTa                                              >  2:4051835/1‑150 (MQ=255)
                                                                                                       aaaCATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATAGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGTCCTCCGATAATGTa                                              >  2:1506216/1‑150 (MQ=255)
                                                                                                           aTGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGTCCTCCGATAATGTAACtt                                          >  1:1766097/1‑150 (MQ=255)
                                                                                                            tggtggCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGTCCTCCGATAATGTAACttt                                         >  1:2562725/1‑150 (MQ=255)
                                                                                                                   aTTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGTCCTCCGATAATGTAACTTTCCACATg                                  >  2:3194588/1‑150 (MQ=255)
                                                                                                                      ggTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGTCCTCCGATAATGTAACTTTCCACATgggg                               <  2:1414360/150‑1 (MQ=255)
                                                                                                                      ggTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGTCCTCCGATAATGTAACTTTCCACATgggg                               <  2:3444328/150‑1 (MQ=255)
                                                                                                                      ggTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGTCCTCCGATAATGTAACTTTCCACATgggg                               >  1:2036848/1‑150 (MQ=255)
                                                                                                                        tGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGTCCTCCGATAATGTAACTTTCCACATGGGGaa                             <  2:2913043/150‑1 (MQ=255)
                                                                                                                         gATCCCCGTCTCCAGCGACAGCGCATAGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGTCCTCCGATAATGTAACTTTCCACATGGGGaaa                            >  1:704257/1‑150 (MQ=255)
                                                                                                                          aTCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGATAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGTCCTCCGATAATGTAACTTTCCACATGGGGaaaa                           <  1:3277879/150‑1 (MQ=255)
                                                                                                                          aTCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGTCCTCCGATAATGTAACTTTCCACATGGGGaaaa                           >  2:130157/1‑150 (MQ=255)
                                                                                                                          aTCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGTCCTCCGATAATGTAACTTTCCACATGGGGaaaa                           <  1:2226061/150‑1 (MQ=255)
                                                                                                                          aTCCCCGTCTCCAGCGACAGCGCATAGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGTCCTCCGATAATGTAACTTTCCACATGGGGaaaa                           <  1:2468506/150‑1 (MQ=255)
                                                                                                                          aTCCCCGTCTCCAGCGACAGCGCATAGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGTCCTCCGATAATGTAACTTTCCACATGGGGaaaa                           <  1:154551/150‑1 (MQ=255)
                                                                                                                           tCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGTCCTCCGATAATGTAACTTTCCACATGGGGAAAAc                          <  1:1976609/150‑1 (MQ=255)
                                                                                                                              ccGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGTCCTCCGATAATGTAACTTTCCACATGGGGAAAACAGc                       >  1:447097/1‑150 (MQ=255)
                                                                                                                                      aGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGTCCTCCGATAATGTAACTTTCCACATGGGGAAAACAGCCCAGTTCa               >  1:1107722/1‑150 (MQ=255)
                                                                                                                                      aGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGTCCTCCGATAATGTAACTTTCCACATGGGGAAAACAGCCCAGTTCa               >  1:3883831/1‑150 (MQ=255)
                                                                                                                                        cGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGTCCTCCGATAATGTAACTTTCCACATGGGGAAAACAGCCCAGTTCAtt             <  2:1676546/150‑1 (MQ=255)
                                                                                                                                        cGACAGCGCATAGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGTCCTCCGATAATGTAACTTTCCACATGGGGAAAACAGCCCAGTTCAtt             >  2:719501/1‑150 (MQ=255)
                                                                                                                                            aGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGTCCTCCGATAATGTAACTTTCCACATGGGGAAAACAGCCCAGTTCATTGGCg         >  1:1189563/1‑150 (MQ=255)
                                                                                                                                              cgcATAGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGTCCTCCGATAATGTAACTTTCCACATGGGGAAAACAGCCCAGTTCATTGGCGAt       <  1:4065949/150‑1 (MQ=255)
                                                                                                                                               gcATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGTCCTCCGATAATGTAACTGTCCACATGGGGAAAACAGCCCAGTTCATTGGCGAtt      <  2:3085548/150‑1 (MQ=255)
                                                                                                                                               gcATAGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGAGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGTCCTCCGATAATGTAACTTTCCACATGGGGAAAACAGCCCAGTTCATTGGCGAtt      >  2:540847/1‑150 (MQ=255)
                                                                                                                                                cATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGTCCTCCGATAATGTAACTTTCCACATGGGGAAAACAGCCCAGTTCATTGGCGAttt     <  1:644566/150‑1 (MQ=255)
                                                                                                                                                cATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGTCCTCCGATAATGTAACTTTCCACATGGGGAAAACAGCCCAGTTCATTGGCGAttt     >  2:634421/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GTGTTGCTACAGCATGCGGCTCAACGGCATCAAATTTTGAGTGATCGGCCATGAGAATAACTTCACGCGAGCGTTGAATAATCCGTGTTTTCACCCCGACCTCAAACATGGTGGCATTGGTGATCCCCGTCTCCAGCGACAGCGCATCGCAGGAGATAAACGCACGTTCTACGGAAAAGGCATTGATCATCTCCAGCGCCAGACTCTCTCCTACCGAAAAATAACCTGGGCGAATCAGTCCTCCGATAATGTAACTTTCCACATGGGGAAAACAGCCCAGTTCATTGGCGATTTTAA  >  CP009273/4516030‑4516326

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