Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 158,287 C→A 94.9% D78Y (GAC→TAC)  ligT ← 2'‑5' RNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273158,2870CA94.9% 1112.5 / 41.0 336D78Y (GAC→TAC) ligT2'‑5' RNA ligase
Reads supporting (aligned to +/- strand):  ref base C (7/10);  new base A (171/147);  total (178/158)
Fisher's exact test for biased strand distribution p-value = 3.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 343 total aligned reads displayed.

TAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGA  >  CP009273/158140‑158433
                                                                                                                                                   |                                                                                                                                                  
tAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGtaat                                                                                                                                                  <  1:1935767/150‑1 (MQ=255)
 aaTATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGtaatc                                                                                                                                                 <  2:701537/150‑1 (MQ=255)
  atatGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGtaatca                                                                                                                                                >  1:44821/1‑150 (MQ=255)
   tatGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCaa                                                                                                                                               >  1:23916/1‑150 (MQ=255)
   tatGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCaa                                                                                                                                               >  2:1451916/1‑150 (MQ=255)
    atGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAg                                                                                                                                              >  1:805287/1‑150 (MQ=255)
    atGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAg                                                                                                                                              >  2:2762530/1‑150 (MQ=255)
        ggATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTg                                                                                                                                          >  2:2627641/1‑150 (MQ=255)
         gATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGa                                                                                                                                         >  1:1234259/1‑150 (MQ=255)
          aTGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAg                                                                                                                                        <  1:1613717/150‑1 (MQ=255)
          aTGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAg                                                                                                                                        >  2:4387822/1‑150 (MQ=255)
            gAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAgtg                                                                                                                                      <  1:2226230/150‑1 (MQ=255)
             aaaCGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAgtgt                                                                                                                                     >  2:3779212/1‑150 (MQ=255)
                 ggACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGaaa                                                                                                                                 >  1:367709/1‑150 (MQ=255)
                 ggACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGaaa                                                                                                                                 <  2:2189574/150‑1 (MQ=255)
                 ggACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGaaa                                                                                                                                 >  1:4589150/1‑150 (MQ=255)
                  gACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAAc                                                                                                                                <  2:608013/150‑1 (MQ=255)
                   aCGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAAcc                                                                                                                               >  1:1046051/1‑150 (MQ=255)
                   aCGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAAcc                                                                                                                               >  2:4209103/1‑150 (MQ=255)
                    cGAGTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCa                                                                                                                              <  1:415286/150‑1 (MQ=255)
                            ttGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACg                                                                                                                      <  2:3465165/150‑1 (MQ=255)
                              gAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGaa                                                                                                                    >  2:4458160/1‑150 (MQ=255)
                               aaaaCAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAAt                                                                                                                   >  1:3605833/1‑150 (MQ=255)
                                aaaCAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAAtc                                                                                                                  >  1:4233684/1‑150 (MQ=255)
                                 aacaacCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAAtcc                                                                                                                 <  2:2921379/150‑1 (MQ=255)
                                 aacaacCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAAtcc                                                                                                                 <  2:4485313/150‑1 (MQ=255)
                                   caacCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAAtccgt                                                                                                               >  2:4718359/1‑150 (MQ=255)
                                    aacCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAAtccgtc                                                                                                              >  2:1557856/1‑150 (MQ=255)
                                     acCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAAtccgtcc                                                                                                             <  1:3318610/150‑1 (MQ=255)
                                       cGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCgg                                                                                                           >  1:2949950/1‑150 (MQ=255)
                                       cGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCgg                                                                                                           >  1:559278/1‑150 (MQ=255)
                                          tGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTa                                                                                                        >  1:1160528/1‑150 (MQ=255)
                                            ctggCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTaaa                                                                                                      >  1:2808592/3‑150 (MQ=255)
                                                cAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAagag                                                                                                  >  2:3926335/1‑150 (MQ=255)
                                                cAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAagag                                                                                                  >  2:2744637/1‑150 (MQ=255)
                                                    cTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAg                                                                                              >  1:883627/1‑150 (MQ=255)
                                                     tgtgAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAgc                                                                                             >  2:3337501/1‑150 (MQ=255)
                                                       tgAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAgcgc                                                                                           >  2:2594646/1‑150 (MQ=255)
                                                       tgAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAgcgc                                                                                           >  1:3175447/1‑150 (MQ=255)
                                                        gAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGcc                                                                                          >  1:3667455/1‑150 (MQ=255)
                                                         aaCGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCt                                                                                         <  2:2512026/150‑1 (MQ=255)
                                                         aaCGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCt                                                                                         >  1:1589688/1‑150 (MQ=255)
                                                          aCGGAGCATATTCGCCAGCTGGATTAAGCCGCGAGGCGGCTGACGCATCCCTAACAACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCtt                                                                                        >  1:4288688/1‑150 (MQ=255)
                                                             gAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTtctc                                                                                     <  1:2503315/150‑1 (MQ=255)
                                                               gCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCtt                                                                                   <  2:867462/150‑1 (MQ=255)
                                                                 atatTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTtc                                                                                 >  2:3225026/1‑150 (MQ=255)
                                                                   atTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTtctc                                                                               <  1:1136877/150‑1 (MQ=255)
                                                                    ttCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCg                                                                              <  2:3003844/150‑1 (MQ=255)
                                                                      cGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGgc                                                                            >  1:1749666/1‑150 (MQ=255)
                                                                      cGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGgc                                                                            >  2:1002023/1‑150 (MQ=255)
                                                                       gCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGgcg                                                                           <  2:1940164/150‑1 (MQ=255)
                                                                        ccAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGgcgc                                                                          >  2:2485340/1‑150 (MQ=255)
                                                                         cAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCt                                                                         >  2:4762684/1‑150 (MQ=255)
                                                                           gCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCCGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCa                                                                       <  2:205733/150‑1 (MQ=255)
                                                                           gCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCa                                                                       <  2:1600759/150‑1 (MQ=255)
                                                                            cTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCAc                                                                      <  2:4778340/150‑1 (MQ=255)
                                                                            cTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCAc                                                                      >  1:4020177/1‑150 (MQ=255)
                                                                             tGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACt                                                                     <  1:2756233/150‑1 (MQ=255)
                                                                             tGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACt                                                                     <  2:1690847/150‑1 (MQ=255)
                                                                              ggATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACtt                                                                    <  1:4302568/150‑1 (MQ=255)
                                                                              ggATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACtt                                                                    >  2:2266905/1‑150 (MQ=255)
                                                                                  tAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGcc                                                                >  2:4066873/1‑150 (MQ=255)
                                                                                   aaGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCt                                                               <  2:805287/150‑1 (MQ=255)
                                                                                   aaGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCt                                                               <  1:4075396/150‑1 (MQ=255)
                                                                                    aGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCAATTCTCCTa                                                              >  1:210399/1‑150 (MQ=255)
                                                                                      ccGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTaaa                                                            >  1:4790308/1‑150 (MQ=255)
                                                                                      ccGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTaaa                                                            <  2:2686536/150‑1 (MQ=255)
                                                                                       cgcgTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTaaaa                                                           <  2:1249737/150‑1 (MQ=255)
                                                                                         cgTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAAt                                                         <  1:3782390/150‑1 (MQ=255)
                                                                                             gcggcTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCa                                                     <  1:4433794/150‑1 (MQ=255)
                                                                                              cggcTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCag                                                    >  2:3575857/1‑150 (MQ=255)
                                                                                               ggcTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCaga                                                   >  1:2682531/1‑150 (MQ=255)
                                                                                               ggcTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCaga                                                   >  1:1878123/1‑150 (MQ=255)
                                                                                                  tGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTc                                                <  1:2360218/150‑1 (MQ=255)
                                                                                                  tGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTc                                                >  1:217376/1‑150 (MQ=255)
                                                                                                  tGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTc                                                >  1:939191/1‑150 (MQ=255)
                                                                                                  tGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTc                                                >  1:4545003/1‑150 (MQ=255)
                                                                                                     cGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGa                                             <  2:1831710/150‑1 (MQ=255)
                                                                                                      gCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGAt                                            <  2:1175383/150‑1 (MQ=255)
                                                                                                          cccTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCaa                                        <  2:2788746/150‑1 (MQ=255)
                                                                                                           ccTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCaaa                                       >  1:3486086/1‑150 (MQ=255)
                                                                                                           ccTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCaaa                                       >  2:3092944/1‑150 (MQ=255)
                                                                                                            cTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAAt                                      <  2:2330919/150‑1 (MQ=255)
                                                                                                              aaCCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTa                                    <  1:2981076/150‑1 (MQ=255)
                                                                                                                    accacacGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATcggcg                              >  1:3867755/1‑150 (MQ=255)
                                                                                                                       acacGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGgcggcg                           >  2:1470201/1‑150 (MQ=255)
                                                                                                                       acacGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGgcggcg                           <  2:3136952/150‑1 (MQ=255)
                                                                                                                          cgcgAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGAcc                        >  2:1371456/1‑150 (MQ=255)
                                                                                                                           gcgAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCg                       <  2:2949950/150‑1 (MQ=255)
                                                                                                                               aCGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACg                   <  2:571691/150‑1 (MQ=255)
                                                                                                                                cGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGt                  >  2:4335967/1‑150 (MQ=255)
                                                                                                                                 gCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAACCCAGGTTGACGCATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTc                 <  2:4703830/150‑1 (MQ=255)
                                                                                                                                   aGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTccc               >  1:949559/1‑150 (MQ=255)
                                                                                                                                   aGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTccc               <  2:4541424/150‑1 (MQ=255)
                                                                                                                                    gCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCg              >  1:186740/1‑150 (MQ=255)
                                                                                                                                    gCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCg              <  1:620402/150‑1 (MQ=255)
                                                                                                                                     ccATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGc             >  1:168049/1‑150 (MQ=255)
                                                                                                                                        ttGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTc          >  1:1545441/1‑150 (MQ=255)
                                                                                                                                         tGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCa         <  1:3105689/150‑1 (MQ=255)
                                                                                                                                                cGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGa  <  2:4592904/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGA  >  CP009273/158140‑158433

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: ‑