Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 425,256 C→T 6.6% P99L (CCG→CTG)  secF → SecYEG protein translocase auxillary subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273425,2560CT 6.6% 1116.5 / 44.8 318P99L (CCG→CTG) secFSecYEG protein translocase auxillary subunit
Reads supporting (aligned to +/- strand):  ref base C (160/137);  new base T (6/15);  total (166/152)
Fisher's exact test for biased strand distribution p-value = 3.97e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 330 total aligned reads displayed.

GCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATG  >  CP009273/425107‑425404
                                                                                                                                                     |                                                                                                                                                    
gCTTGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATgcc                                                                                                                                                      <  2:4605803/150‑1 (MQ=255)
gCTGGATGTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCAGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATgcc                                                                                                                                                      <  2:5239736/150‑1 (MQ=255)
  tGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATgccgc                                                                                                                                                    >  1:3160408/1‑150 (MQ=255)
   ggATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATgccgcc                                                                                                                                                   <  2:2557398/150‑1 (MQ=255)
          aCCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGaa                                                                                                                                            <  1:4094685/150‑1 (MQ=255)
            cGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAgg                                                                                                                                          >  1:4762062/1‑150 (MQ=255)
             ggtggtACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGc                                                                                                                                         >  1:1355581/1‑150 (MQ=255)
              gtggtACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCg                                                                                                                                        <  1:5615359/150‑1 (MQ=255)
              gtggtACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCg                                                                                                                                        >  2:5472526/1‑150 (MQ=255)
              gtggtACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCg                                                                                                                                        <  1:1905257/150‑1 (MQ=255)
               tggtACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGa                                                                                                                                       <  2:2319831/150‑1 (MQ=255)
                  tACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAAc                                                                                                                                    >  2:126127/1‑150 (MQ=255)
                   aCGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAAcc                                                                                                                                   <  2:1232868/150‑1 (MQ=255)
                   aCGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAAcc                                                                                                                                   <  1:2582165/150‑1 (MQ=255)
                     ggTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACcgg                                                                                                                                 >  2:5510069/1‑150 (MQ=255)
                      gTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACcggc                                                                                                                                >  1:2103132/1‑150 (MQ=255)
                         aTTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGt                                                                                                                             >  1:1964437/1‑150 (MQ=255)
                             aaaTTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAgg                                                                                                                         >  1:3611436/1‑150 (MQ=255)
                              aaTTACGCTCGAAACACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGgt                                                                                                                        <  2:5054879/150‑1 (MQ=255)
                                 tACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTGTGGTAGCAGCCATGACATCATGGTCCGAATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGgtggt                                                                                                                     >  2:3753594/1‑148 (MQ=255)
                                 tACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGAAAAAAGCCGGTTATGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGtt                                                                                                                     >  2:908868/1‑150 (MQ=255)
                                  aCGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTg                                                                                                                    <  2:1685091/150‑1 (MQ=255)
                                        gAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCTGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGGTGGGCAGc                                                                                                              >  1:931506/1‑150 (MQ=255)
                                          aaaaCCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCa                                                                                                            <  1:4465779/150‑1 (MQ=255)
                                          aaaaCCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCa                                                                                                            <  1:3900050/150‑1 (MQ=255)
                                          aaaaCCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCa                                                                                                            >  1:5515148/1‑150 (MQ=255)
                                          aaaaCCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCa                                                                                                            <  1:3014149/150‑1 (MQ=255)
                                             aCCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGt                                                                                                         <  1:684130/150‑1 (MQ=255)
                                                ggCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCt                                                                                                      <  1:838792/150‑1 (MQ=255)
                                                 gCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCTGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTg                                                                                                     <  2:1609773/150‑1 (MQ=255)
                                                      aaTTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCTGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGt                                                                                                >  2:1784281/1‑150 (MQ=255)
                                                       aTTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTg                                                                                               >  2:2162768/1‑150 (MQ=255)
                                                        ttGACGTAATGCGTGATGCATTTCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGa                                                                                              <  2:3803530/150‑1 (MQ=255)
                                                        ttGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGa                                                                                              >  2:2937527/1‑150 (MQ=255)
                                                        ttGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGa                                                                                              >  1:1552980/1‑150 (MQ=255)
                                                         tGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGAt                                                                                             >  2:4499169/1‑150 (MQ=255)
                                                          gACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCTGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGAtt                                                                                            <  1:775349/150‑1 (MQ=255)
                                                          gACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGAtt                                                                                            <  1:1193314/150‑1 (MQ=255)
                                                          gACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGAtt                                                                                            >  1:5351207/1‑150 (MQ=255)
                                                            cGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTaa                                                                                          >  2:1107105/1‑150 (MQ=255)
                                                            cGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTaa                                                                                          >  1:226674/1‑150 (MQ=255)
                                                            cGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTaa                                                                                          <  2:3921564/150‑1 (MQ=255)
                                                              tAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACg                                                                                        >  1:114136/1‑150 (MQ=255)
                                                                aTGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGaa                                                                                      >  2:2334115/1‑150 (MQ=255)
                                                                  gCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATc                                                                                    >  2:1323323/1‑150 (MQ=255)
                                                                  gCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGAAGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACTAATc                                                                                    >  1:1623236/1‑150 (MQ=255)
                                                                    gTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATcca                                                                                  >  2:205816/1‑150 (MQ=255)
                                                                    gTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATcca                                                                                  >  2:1192936/1‑150 (MQ=255)
                                                                        tGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCaa                                                                              <  2:2987970/150‑1 (MQ=255)
                                                                         gCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAAt                                                                             <  2:623305/150‑1 (MQ=255)
                                                                         gCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAAt                                                                             <  2:1495299/150‑1 (MQ=255)
                                                                            ttGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAg                                                                          >  2:1068580/1‑150 (MQ=255)
                                                                            ttGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGGTCTGAAGGTGATTAACGAATCCACCAATCAg                                                                          >  2:759249/1‑150 (MQ=255)
                                                                              gCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGaa                                                                        >  2:4012730/1‑150 (MQ=255)
                                                                                aaaaaGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATg                                                                      >  1:125819/1‑150 (MQ=255)
                                                                                aaaaaGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATg                                                                      >  1:2369431/1‑150 (MQ=255)
                                                                                     gCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATgcagca                                                                 <  2:3977730/150‑1 (MQ=255)
                                                                                      ccGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGcagcag                                                                >  1:5247040/1‑150 (MQ=255)
                                                                                         gTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGa                                                             >  1:972179/1‑150 (MQ=255)
                                                                                          ttttGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCTGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGaa                                                            <  1:547894/150‑1 (MQ=255)
                                                                                             tGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCg                                                         <  1:2520769/150‑1 (MQ=255)
                                                                                               aaGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTa                                                       >  2:3115068/1‑150 (MQ=255)
                                                                                                agagCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTAt                                                      >  2:565793/1‑150 (MQ=255)
                                                                                                 gagCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTAtt                                                     <  1:4147259/150‑1 (MQ=255)
                                                                                                 gagCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTAtt                                                     <  2:683728/150‑1 (MQ=255)
                                                                                                   gCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGa                                                   >  1:1028647/1‑150 (MQ=255)
                                                                                                     cGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGt                                                 >  2:2689437/1‑150 (MQ=255)
                                                                                                     cGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGt                                                 >  1:4691953/1‑150 (MQ=255)
                                                                                                      gATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGtt                                                >  1:2135035/1‑150 (MQ=255)
                                                                                                          ctgcAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTtcgtc                                            >  2:242439/1‑150 (MQ=255)
                                                                                                            gcAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCgg                                          <  2:4970587/150‑1 (MQ=255)
                                                                                                             cAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGt                                         >  2:5134744/1‑150 (MQ=255)
                                                                                                                aaCTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCg                                      <  2:3096158/150‑1 (MQ=255)
                                                                                                                aaCTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCg                                      >  1:4258602/1‑150 (MQ=255)
                                                                                                                 aCTTTGGTAGCAGCCATGACATCATGGTCCGTATGCTGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGa                                     <  1:4870679/150‑1 (MQ=255)
                                                                                                                  cTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTCTTGAGTTCGTCGGTCCGAg                                    >  2:430302/1‑150 (MQ=255)
                                                                                                                  cTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAg                                    >  1:1708472/1‑150 (MQ=255)
                                                                                                                   tttGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGATACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGc                                   <  1:5004918/150‑1 (MQ=255)
                                                                                                                   tttGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGc                                   <  2:1003831/150‑1 (MQ=255)
                                                                                                                    ttGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCg                                  >  2:847570/1‑150 (MQ=255)
                                                                                                                       gTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTgg                               >  2:5316452/1‑150 (MQ=255)
                                                                                                                       gTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTgg                               <  1:2343741/150‑1 (MQ=255)
                                                                                                                        tAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTggg                              >  1:4307785/1‑150 (MQ=255)
                                                                                                                        tAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTggg                              <  1:2872411/150‑1 (MQ=255)
                                                                                                                        tAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTggg                              <  2:4130355/150‑1 (MQ=255)
                                                                                                                         agcagcCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTgggg                             <  2:349221/150‑1 (MQ=255)
                                                                                                                         agcagcCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTgggg                             >  2:319125/1‑150 (MQ=255)
                                                                                                                          gcagcCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTggggg                            >  2:1570109/1‑150 (MQ=255)
                                                                                                                           cagcCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGc                           <  1:1849410/150‑1 (MQ=255)
                                                                                                                             gcCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAg                         <  1:3506484/150‑1 (MQ=255)
                                                                                                                                     atcatGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTgcg                 >  2:4913186/1‑150 (MQ=255)
                                                                                                                                       catGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCa               <  1:2685499/150‑1 (MQ=255)
                                                                                                                                         tGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCaaa             >  1:146627/1‑150 (MQ=255)
                                                                                                                                            tCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCg          >  1:5156768/1‑150 (MQ=255)
                                                                                                                                              cGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGGGAAGCGTATGGAGTTCGTCGGTCCGAGCGGGGGGGCAGTCCTTGCGCAAACCGGt        >  2:3860694/1‑150 (MQ=255)
                                                                                                                                               gTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTg       <  1:384932/150‑1 (MQ=255)
                                                                                                                                               gTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTg       >  2:1908351/1‑150 (MQ=255)
                                                                                                                                                 aTGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCg     >  2:5472829/1‑150 (MQ=255)
                                                                                                                                                  tGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGa    <  1:3520861/150‑1 (MQ=255)
                                                                                                                                                    ccgccTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATg  <  1:3548987/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATG  >  CP009273/425107‑425404

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: