Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 158,287 C→A 100% 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,2870CA95.8% 866.3 / 25.7 265D78Y (GAC→TAC) ligT2'‑5' RNA ligase
Reads supporting (aligned to +/- strand):  ref base C (6/5);  new base A (122/132);  total (128/137)
Fisher's exact test for biased strand distribution p-value = 7.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 270 total aligned reads displayed.

AATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGG  >  CP009273/158141‑158432
                                                                                                                                                  |                                                                                                                                                 
aaTATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGtaatc                                                                                                                                                >  1:1127917/1‑150 (MQ=255)
aaTATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGtaatc                                                                                                                                                >  2:16638/1‑150 (MQ=255)
  tatGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCaa                                                                                                                                              >  1:657911/1‑150 (MQ=255)
   atGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAg                                                                                                                                             <  2:1250982/150‑1 (MQ=255)
    tgtgGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGc                                                                                                                                            <  1:3667452/150‑1 (MQ=255)
    tgtgGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGc                                                                                                                                            <  1:3565117/150‑1 (MQ=255)
      tgGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGt                                                                                                                                          >  1:327684/1‑150 (MQ=255)
      tgGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGt                                                                                                                                          >  1:2971147/1‑150 (MQ=255)
         aTGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAg                                                                                                                                       <  1:602796/150‑1 (MQ=255)
         aTGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAg                                                                                                                                       >  2:2915366/1‑150 (MQ=255)
              aCGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGa                                                                                                                                  >  1:241474/1‑150 (MQ=255)
               cGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGaa                                                                                                                                 <  2:309530/150‑1 (MQ=255)
                 gACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAAc                                                                                                                               >  2:1536845/1‑150 (MQ=255)
                 gACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGAGAACGGAGCATATTCGCCAGCTGGATAAAGCCGAGAGGCGGCTGAAGCATCACTAACCACACCACACGCGAACGCAGCCATAGTCCCGCGTAATCAAGCGTGAGTGTGaaca                                                                                                                               >  2:1536841/1‑148 (MQ=255)
                  aCGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAAcc                                                                                                                              <  2:3555607/150‑1 (MQ=255)
                    gATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAg                                                                                                                            <  2:3000538/150‑1 (MQ=255)
                        gCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTg                                                                                                                        <  1:1663917/150‑1 (MQ=255)
                          tttGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGAc                                                                                                                      >  1:1688534/1‑150 (MQ=255)
                           ttGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACg                                                                                                                     >  2:4376193/1‑150 (MQ=255)
                              aaaaCAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAAt                                                                                                                  <  1:516116/150‑1 (MQ=255)
                                aacaacCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAAtcc                                                                                                                <  2:1164556/150‑1 (MQ=255)
                                 acaacCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAAtccg                                                                                                               >  1:42225/1‑150 (MQ=255)
                                   aacCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAAtccgtc                                                                                                             >  1:3732689/1‑150 (MQ=255)
                                   aacCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAAtccgtc                                                                                                             >  1:3980247/1‑150 (MQ=255)
                                    acCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAAtccgtcc                                                                                                            <  1:4563528/150‑1 (MQ=255)
                                    acCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAAtccgtcc                                                                                                            >  1:2337913/1‑150 (MQ=255)
                                      cGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCgg                                                                                                          <  1:708270/150‑1 (MQ=255)
                                      cGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCgg                                                                                                          <  1:4081953/150‑1 (MQ=255)
                                        cTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAGGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCt                                                                                                        <  1:83489/150‑1 (MQ=255)
                                         tGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTa                                                                                                       <  2:4270448/150‑1 (MQ=255)
                                           cGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTaaa                                                                                                     >  2:874878/1‑150 (MQ=255)
                                              gCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAaga                                                                                                  <  1:1368610/150‑1 (MQ=255)
                                              gCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAaga                                                                                                  >  1:4636124/1‑150 (MQ=255)
                                               cAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAagag                                                                                                 <  2:2135682/150‑1 (MQ=255)
                                               cAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAagag                                                                                                 <  2:4463137/150‑1 (MQ=255)
                                                   cTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAg                                                                                             <  2:4195536/150‑1 (MQ=255)
                                                   cTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAg                                                                                             <  1:234944/150‑1 (MQ=255)
                                                    tgtgAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGATAAAgc                                                                                            >  1:1885239/1‑150 (MQ=255)
                                                    tgtgAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAgc                                                                                            >  2:3363745/1‑150 (MQ=255)
                                                      tgAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAgcgc                                                                                          >  2:3967380/1‑150 (MQ=255)
                                                       gAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGcc                                                                                         >  2:18482/1‑150 (MQ=255)
                                                        aaCGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCt                                                                                        <  2:443745/150‑1 (MQ=255)
                                                        aaCGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCt                                                                                        <  1:2046582/150‑1 (MQ=255)
                                                        aaCGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCt                                                                                        >  2:2087418/1‑150 (MQ=255)
                                                         aCGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCtt                                                                                       <  2:1648808/150‑1 (MQ=255)
                                                           ggAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTtct                                                                                     <  1:1965677/150‑1 (MQ=255)
                                                             aGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTctct                                                                                   <  1:3745537/150‑1 (MQ=255)
                                                               cATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCttt                                                                                 <  2:971147/150‑1 (MQ=255)
                                                                  atTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTtctc                                                                              >  1:2906079/1‑150 (MQ=255)
                                                                  agtCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTtctc                                                                              <  2:3922992/148‑1 (MQ=255)
                                                                     cGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGgc                                                                           <  2:386645/150‑1 (MQ=255)
                                                                     cGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGgc                                                                           <  1:2650161/150‑1 (MQ=255)
                                                                         aGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTc                                                                       <  1:2064550/150‑1 (MQ=255)
                                                                          gCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCa                                                                      <  2:670999/150‑1 (MQ=255)
                                                                           cTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCAc                                                                     >  1:3025833/1‑150 (MQ=255)
                                                                             ggATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACtt                                                                   >  2:4213822/1‑150 (MQ=255)
                                                                             ggATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACtt                                                                   >  2:2251482/1‑150 (MQ=255)
                                                                              gATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTc                                                                  >  2:739867/1‑150 (MQ=255)
                                                                               aTTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCg                                                                 <  1:3474135/150‑1 (MQ=255)
                                                                               aTTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCg                                                                 >  1:2510716/1‑150 (MQ=255)
                                                                               aTTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCg                                                                 <  1:3506117/150‑1 (MQ=255)
                                                                                   aGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTa                                                             >  1:3101138/1‑150 (MQ=255)
                                                                                     ccGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTaaa                                                           >  2:132353/1‑150 (MQ=255)
                                                                                      cgcgTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTaaaa                                                          >  2:2532421/1‑150 (MQ=255)
                                                                                        cgTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAAt                                                        >  1:2218077/1‑150 (MQ=255)
                                                                                        cgTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAAt                                                        <  1:3905464/150‑1 (MQ=255)
                                                                                        cgTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAAt                                                        <  1:3905516/150‑1 (MQ=255)
                                                                                        cgTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAAt                                                        <  1:4265918/150‑1 (MQ=255)
                                                                                         gTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATg                                                       <  1:4434544/150‑1 (MQ=255)
                                                                                             cggcTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCag                                                   >  2:1450332/1‑150 (MQ=255)
                                                                                              ggcTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCaga                                                  >  1:2535171/1‑150 (MQ=255)
                                                                                              ggcTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCaga                                                  <  1:4501996/150‑1 (MQ=255)
                                                                                                cTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGt                                                >  2:253565/1‑150 (MQ=255)
                                                                                                 tGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTc                                               >  1:4307900/1‑150 (MQ=255)
                                                                                                 tGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTc                                               >  1:4308065/1‑150 (MQ=255)
                                                                                                   aCGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAg                                             <  1:149822/150‑1 (MQ=255)
                                                                                                   aCGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAg                                             <  1:3152537/150‑1 (MQ=255)
                                                                                                    cGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGa                                            >  1:2722994/1‑150 (MQ=255)
                                                                                                         cccTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCaa                                       >  2:3075353/1‑150 (MQ=255)
                                                                                                          ccTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCaaa                                      <  1:3498879/150‑1 (MQ=255)
                                                                                                           cTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAAt                                     <  1:3383791/150‑1 (MQ=255)
                                                                                                            tAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAAtt                                    <  1:3366747/150‑1 (MQ=255)
                                                                                                              accacaccacACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTAt                                  >  2:4049021/1‑150 (MQ=255)
                                                                                                              accacaccacACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTAt                                  >  2:3840266/1‑150 (MQ=255)
                                                                                                               ccacaccacaCGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATc                                 >  2:385156/1‑150 (MQ=255)
                                                                                                                cacaccacacGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATcg                                <  2:2655303/150‑1 (MQ=255)
                                                                                                                cacaccacacGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATcg                                >  1:1006963/1‑150 (MQ=255)
                                                                                                                    ccacacGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATcggcgg                            <  2:2338575/150‑1 (MQ=255)
                                                                                                                     cacacGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCggcggc                           <  2:2863797/150‑1 (MQ=255)
                                                                                                                      acacGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGgcggcg                          <  1:4022553/150‑1 (MQ=255)
                                                                                                                            gAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGa                    >  2:2772954/1‑150 (MQ=255)
                                                                                                                              aCGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACg                  >  2:3310917/1‑150 (MQ=255)
                                                                                                                                gCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTc                >  2:321458/1‑150 (MQ=255)
                                                                                                                                  aGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTccc              <  1:989156/150‑1 (MQ=255)
                                                                                                                                     cATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGcc           >  1:3189303/1‑150 (MQ=255)
                                                                                                                                      aTTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGTGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCt          >  2:2361277/1‑150 (MQ=255)
                                                                                                                                        tGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCa        >  1:1815310/1‑150 (MQ=255)
                                                                                                                                          tCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAgg      <  1:1779753/150‑1 (MQ=255)
                                                                                                                                            cGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTg    <  2:549306/150‑1 (MQ=255)
                                                                                                                                             ggCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTgg   <  1:750014/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
AATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGG  >  CP009273/158141‑158432

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: