Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,880,731 G→T 100% D182Y (GAT→TAT)  mnmE → tRNA U34 5‑methylaminomethyl‑2‑thiouridine modification GTPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,880,7310GT97.3% 812.6 / 9.8 256D182Y (GAT→TAT) mnmEtRNA U34 5‑methylaminomethyl‑2‑thiouridine modification GTPase
Reads supporting (aligned to +/- strand):  ref base G (0/7);  new base T (94/155);  total (94/162)
Fisher's exact test for biased strand distribution p-value = 4.95e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 256 total aligned reads displayed.

TCGACGCCAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCGATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGTAAATC  >  CP009273/3880585‑3880876
                                                                                                                                                  |                                                                                                                                                 
tCGACGCCAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTAt                                                                                                                                                 <  2:3966915/149‑1 (MQ=255)
  gACGCCAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTAttt                                                                                                                                               <  1:1003241/149‑1 (MQ=255)
  gACGCCAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTAttt                                                                                                                                               >  2:87410/1‑149 (MQ=255)
  gACGCCAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGAGGAAGAGATCTAttt                                                                                                                                               >  2:119487/1‑149 (MQ=255)
   aCGCCAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTc                                                                                                                                              >  1:345648/1‑149 (MQ=255)
   aCGCCAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAAATGATTTCCCCGATGAAGAGATCTATTTc                                                                                                                                              >  2:2385504/1‑149 (MQ=255)
    cGCCAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTcc                                                                                                                                             <  2:3496360/149‑1 (MQ=255)
        aGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCtctc                                                                                                                                         <  1:328835/149‑1 (MQ=255)
        aGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCtctc                                                                                                                                         <  1:143944/149‑1 (MQ=255)
         gTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTcc                                                                                                                                        <  1:2709592/149‑1 (MQ=255)
         gTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTcc                                                                                                                                        <  1:650454/149‑1 (MQ=255)
              gAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACgg                                                                                                                                   >  2:3761261/1‑149 (MQ=255)
               aaCAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGa                                                                                                                                  <  1:1514632/149‑1 (MQ=255)
                aCAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGaa                                                                                                                                 >  1:2176107/1‑149 (MQ=255)
                   ggcggcCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGaaaaa                                                                                                                              >  1:3686591/1‑149 (MQ=255)
                       gcCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGa                                                                                                                          <  1:3703136/149‑1 (MQ=255)
                           gTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGcc                                                                                                                      >  2:931646/1‑149 (MQ=255)
                            ttCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGccc                                                                                                                     <  1:1656340/149‑1 (MQ=255)
                              cGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGGTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAg                                                                                                                   >  1:3005601/1‑149 (MQ=255)
                                   cTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCaa                                                                                                              <  2:1549943/149‑1 (MQ=255)
                                      aaCTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGa                                                                                                           >  1:1747034/1‑149 (MQ=255)
                                      aaCTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGa                                                                                                           <  1:3933260/149‑1 (MQ=255)
                                          cGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGtt                                                                                                       <  1:2871972/149‑1 (MQ=255)
                                            cTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTAt                                                                                                     <  2:302000/149‑1 (MQ=255)
                                               cAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCGATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGc                                                                                                  <  2:2228899/149‑1 (MQ=255)
                                                aaGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGcc                                                                                                 <  1:1545872/149‑1 (MQ=255)
                                                    cgcATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATc                                                                                             <  1:4025397/149‑1 (MQ=255)
                                                    cgcATTCTCCGCACGGGTTAATCATCAGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATc                                                                                             <  2:3443171/149‑1 (MQ=255)
                                                     gcATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCCATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCt                                                                                            >  2:2932577/1‑149 (MQ=255)
                                                      cATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCtt                                                                                           <  1:30482/149‑1 (MQ=255)
                                                        ttCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGa                                                                                         >  2:4028516/1‑149 (MQ=255)
                                                          ctcCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATg                                                                                       <  1:564800/149‑1 (MQ=255)
                                                           tcCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGc                                                                                      <  1:859022/149‑1 (MQ=255)
                                                            ccGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCa                                                                                     <  2:1638985/149‑1 (MQ=255)
                                                             cGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAg                                                                                    <  1:1009794/149‑1 (MQ=255)
                                                             cGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAg                                                                                    <  1:2462178/149‑1 (MQ=255)
                                                                aCGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAgtgc                                                                                 <  1:2140622/149‑1 (MQ=255)
                                                                aCGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAgtgc                                                                                 >  1:1939211/1‑149 (MQ=255)
                                                                 cGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAgtgcg                                                                                <  1:1118651/149‑1 (MQ=255)
                                                                  gggTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAgtgcgt                                                                               <  2:1380993/149‑1 (MQ=255)
                                                                   ggTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAgtgcgtg                                                                              >  1:74019/1‑149 (MQ=255)
                                                                       aaTCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGa                                                                          >  1:134495/1‑149 (MQ=255)
                                                                       aaTCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGa                                                                          >  2:4002055/1‑149 (MQ=255)
                                                                       aaTCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGa                                                                          >  2:1902753/1‑149 (MQ=255)
                                                                         tcatcTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAg                                                                        <  2:2137954/149‑1 (MQ=255)
                                                                         tcatcTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAg                                                                        <  2:4113237/149‑1 (MQ=255)
                                                                          catcTTGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGc                                                                       <  1:2539899/149‑1 (MQ=255)
                                                                          catcTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGc                                                                       <  1:2965791/149‑1 (MQ=255)
                                                                          catcTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGc                                                                       <  2:3431619/149‑1 (MQ=255)
                                                                           atcTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCa                                                                      <  2:1465359/149‑1 (MQ=255)
                                                                            tcTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCAc                                                                     >  2:3703349/1‑149 (MQ=255)
                                                                               ggTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTc                                                                  >  2:1503997/1‑149 (MQ=255)
                                                                               ggTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTc                                                                  <  2:99085/149‑1 (MQ=255)
                                                                                  aGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAgg                                                               >  2:1981163/1‑149 (MQ=255)
                                                                                   gAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAggg                                                              >  1:2431593/1‑149 (MQ=255)
                                                                                     aGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTa                                                            >  2:2633572/1‑149 (MQ=255)
                                                                                      gCCCTCACCCACTTGCGCATGTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAg                                                           <  1:2593315/149‑1 (MQ=255)
                                                                                        ccTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGtt                                                         >  2:2280248/1‑149 (MQ=255)
                                                                                        ccTCACCCACTTGCGCATGTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGtt                                                         <  2:2666684/149‑1 (MQ=255)
                                                                                            aCCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTttgtt                                                     >  2:2996849/1‑149 (MQ=255)
                                                                                            aCCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTttgtt                                                     <  2:472104/149‑1 (MQ=255)
                                                                                              ccACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTgc                                                   >  1:3534066/1‑149 (MQ=255)
                                                                                                 cTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGcgcg                                                <  1:2979947/149‑1 (MQ=255)
                                                                                                 cTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGcgcg                                                >  2:1788646/1‑149 (MQ=255)
                                                                                                 cTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGcgcg                                                >  1:609204/1‑149 (MQ=255)
                                                                                                 cTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGcgcg                                                >  1:3221884/1‑149 (MQ=255)
                                                                                                    gcgcATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCGATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAg                                             <  1:7789/149‑1 (MQ=255)
                                                                                                     cgcAGTTACGTCGACGCGTCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCTAAgt                                            <  2:3580238/149‑2 (MQ=255)
                                                                                                      gcATCTACGTCGAAGCGGCAAGTGATGTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAggg                                           <  2:2173044/149‑1 (MQ=255)
                                                                                                       cATTTACGTCTAAGCGGCATTTGATTTCCCCGATGAAGAGATCTTTTTCCTCTCCTACGGAAAAATTGAAGACAAGATCACTGACGTTTTTGCCGATCTTGATGATTTGATTGGTGAAGCACGTCAGGTTAGTTTGTTGCGCGAAGGGa                                          <  1:2995747/149‑1 (MQ=255)
                                                                                                       cATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCTATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGa                                          <  1:3986683/149‑1 (MQ=255)
                                                                                                        aTTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGAt                                         >  2:737112/1‑149 (MQ=255)
                                                                                                             cGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAg                                    <  1:1037974/149‑1 (MQ=255)
                                                                                                             cGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAg                                    <  1:1331279/149‑1 (MQ=255)
                                                                                                              gTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAgt                                   <  1:1159394/149‑1 (MQ=255)
                                                                                                              gTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAgt                                   <  1:2436131/149‑1 (MQ=255)
                                                                                                                cGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAgtgg                                 <  1:1231054/149‑1 (MQ=255)
                                                                                                                cGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAgtgg                                 <  2:3174028/149‑1 (MQ=255)
                                                                                                                 gAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAgtggt                                <  1:1148701/149‑1 (MQ=255)
                                                                                                                   aGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGa                              >  2:1979186/1‑149 (MQ=255)
                                                                                                                    gCGGCAATTGATTTCCCAGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCCGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGAt                             <  1:157829/149‑1 (MQ=255)
                                                                                                                     cGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGAtt                            >  1:1398987/1‑149 (MQ=255)
                                                                                                                     cGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGAtt                            <  1:2847257/149‑1 (MQ=255)
                                                                                                                       gCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGc                          <  2:544818/149‑1 (MQ=255)
                                                                                                                       gCAATTGATGTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGc                          <  2:1062563/149‑1 (MQ=255)
                                                                                                                        cAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGcc                         <  2:3306034/149‑1 (MQ=255)
                                                                                                                         aaTTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCg                        >  1:3286589/1‑149 (MQ=255)
                                                                                                                         aaTTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCg                        >  1:3872445/1‑149 (MQ=255)
                                                                                                                          aTTGTTTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCgg                       >  2:2497667/1‑149 (MQ=255)
                                                                                                                          aTTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCgg                       <  1:2280248/149‑1 (MQ=255)
                                                                                                                           ttGATGTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGa                      <  1:1210019/149‑1 (MQ=255)
                                                                                                                                   cccGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTaa              <  1:3595813/149‑1 (MQ=255)
                                                                                                                                   cccGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTaa              <  2:797036/149‑1 (MQ=255)
                                                                                                                                   cccGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTaa              <  1:87410/149‑1 (MQ=255)
                                                                                                                                    ccGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAAc             <  2:600031/149‑1 (MQ=255)
                                                                                                                                        tGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCg         >  1:3792276/1‑149 (MQ=255)
                                                                                                                                         gAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCgg        >  1:44455/1‑149 (MQ=255)
                                                                                                                                          aaGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTCGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGt       >  2:2762880/1‑149 (MQ=255)
                                                                                                                                           agagATCTATTTCCTCTCCTTCGGAAATATTGAAGCCCTGATCAAGGTGTTTATTGCCGATCTTGATGCTGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGATGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGTa      <  1:87280/149‑1 (MQ=255)
                                                                                                                                              gaTCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGTAAAt   <  2:1119272/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TCGACGCCAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCGATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGTAAATC  >  CP009273/3880585‑3880876

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: