Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 158,287 C→A 83.5% 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,2870CA83.5% 887.3 / 192.7 352D78Y (GAC→TAC) ligT2'‑5' RNA ligase
Reads supporting (aligned to +/- strand):  ref base C (29/29);  new base A (128/166);  total (157/195)
Fisher's exact test for biased strand distribution p-value = 3.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 358 total aligned reads displayed.

AATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGA  >  CP009273/158141‑158433
                                                                                                                                                  |                                                                                                                                                  
   atGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCaa                                                                                                                                               <  2:1228182/149‑1 (MQ=255)
   atGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCaa                                                                                                                                               <  2:2836757/149‑1 (MQ=255)
   atGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCaa                                                                                                                                               <  1:1714061/149‑1 (MQ=255)
    tgtgGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAg                                                                                                                                              <  1:674626/149‑1 (MQ=255)
    tgtgGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAg                                                                                                                                              <  1:3086122/149‑1 (MQ=255)
     gtgGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGc                                                                                                                                             <  1:3502212/149‑1 (MQ=255)
       ggATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGt                                                                                                                                           >  2:4060081/1‑149 (MQ=255)
        gATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACAACACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTg                                                                                                                                          <  1:820534/149‑1 (MQ=255)
            aaaCGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAgtg                                                                                                                                      >  2:1732840/1‑149 (MQ=255)
             aaCGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAgtgt                                                                                                                                     <  2:3217411/149‑1 (MQ=255)
             aaCGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAgtgt                                                                                                                                     <  2:1675848/149‑1 (MQ=255)
                ggACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGaa                                                                                                                                  >  1:3920775/1‑149 (MQ=255)
                 gACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGaaa                                                                                                                                 >  2:2974886/1‑149 (MQ=255)
                  aCGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAAc                                                                                                                                <  1:4317636/149‑1 (MQ=255)
                   cGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAAcc                                                                                                                               >  2:3205660/1‑149 (MQ=255)
                     aTTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAg                                                                                                                             >  1:475094/1‑149 (MQ=255)
                     aTTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAg                                                                                                                             >  2:313371/1‑149 (MQ=255)
                       tGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGt                                                                                                                           <  1:929146/149‑1 (MQ=255)
                       tGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGt                                                                                                                           <  2:2193030/149‑1 (MQ=255)
                        gCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGtt                                                                                                                          <  2:1878653/149‑1 (MQ=255)
                         cTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTg                                                                                                                         >  1:3787268/1‑149 (MQ=255)
                          tttGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGa                                                                                                                        <  2:2626820/149‑1 (MQ=255)
                          tttGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGa                                                                                                                        <  2:2309734/149‑1 (MQ=255)
                          tttGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGa                                                                                                                        >  2:1861771/1‑149 (MQ=255)
                           ttGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGAc                                                                                                                       <  2:5149619/149‑1 (MQ=255)
                              aaaaCAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGaa                                                                                                                    >  1:951236/1‑149 (MQ=255)
                               aaaCAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAAt                                                                                                                   >  1:2087699/1‑149 (MQ=255)
                                 acaacCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAAtcc                                                                                                                 <  1:1949454/149‑1 (MQ=255)
                                    acCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAAtccgtc                                                                                                              >  1:4894776/1‑149 (MQ=255)
                                     ccGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAAtccgtcc                                                                                                             <  1:1431694/149‑1 (MQ=255)
                                     ccGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAAtccgtcc                                                                                                             <  2:2058976/149‑1 (MQ=255)
                                     ccGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAAtccgtcc                                                                                                             >  1:3158581/1‑149 (MQ=255)
                                     ccGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAAtccgtcc                                                                                                             >  2:1315578/1‑149 (MQ=255)
                                        cTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGc                                                                                                          <  2:4633189/149‑1 (MQ=255)
                                        cTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGc                                                                                                          <  1:4569147/149‑1 (MQ=255)
                                        cTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGc                                                                                                          <  1:4708974/149‑1 (MQ=255)
                                         tGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCt                                                                                                         <  2:821265/149‑1 (MQ=255)
                                         tGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCt                                                                                                         >  2:2716680/1‑149 (MQ=255)
                                             ggCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTaaaa                                                                                                     <  1:4231827/149‑1 (MQ=255)
                                               cAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAaga                                                                                                   <  2:4174982/149‑1 (MQ=255)
                                                 gCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAgaga                                                                                                 <  1:3303211/149‑1 (MQ=255)
                                                 gCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAgaga                                                                                                 >  1:773216/1‑149 (MQ=255)
                                                  ccTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGaa                                                                                                >  2:3006931/1‑149 (MQ=255)
                                                    tgtgAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAg                                                                                              >  1:544594/1‑149 (MQ=255)
                                                         aCGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCt                                                                                         >  2:4365743/1‑149 (MQ=255)
                                                         aCGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCt                                                                                         >  1:1126363/1‑149 (MQ=255)
                                                         aCGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCGGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCt                                                                                         <  1:1368648/149‑1 (MQ=255)
                                                           ggAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTtc                                                                                       <  2:3366732/149‑1 (MQ=255)
                                                            gAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTtct                                                                                      >  1:275550/1‑149 (MQ=255)
                                                             aGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTtctc                                                                                     <  2:1493531/149‑1 (MQ=255)
                                                              gCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTctct                                                                                    <  2:4860450/149‑1 (MQ=255)
                                                              gCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTctct                                                                                    <  2:4815627/149‑1 (MQ=255)
                                                              gCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTctct                                                                                    >  2:1180332/1‑149 (MQ=255)
                                                              gCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTctct                                                                                    <  1:1055300/149‑1 (MQ=255)
                                                               cATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCtt                                                                                   <  1:5087161/149‑1 (MQ=255)
                                                                   ttCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAACGCGCCTTCTCTTtctc                                                                               <  1:4830147/149‑1 (MQ=255)
                                                                   ttCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTtctc                                                                               <  2:4696665/149‑1 (MQ=255)
                                                                   ttCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTtctc                                                                               >  1:3417241/1‑149 (MQ=255)
                                                                    tCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCg                                                                              <  1:46395/149‑1 (MQ=255)
                                                                      gCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGgc                                                                            <  2:4922244/149‑1 (MQ=255)
                                                                          gCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTc                                                                        <  2:349867/149‑1 (MQ=255)
                                                                           cTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCa                                                                       >  2:4536074/1‑149 (MQ=255)
                                                                            tGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCAc                                                                      <  2:808589/149‑1 (MQ=255)
                                                                             ggATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAACGCGCCTTCTCTTTCTCGGCGCTCACt                                                                     <  1:2566362/149‑1 (MQ=255)
                                                                              gATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGACAGCGCCTTCTCTTTCTCGGCGCTCACtt                                                                    <  1:1746116/149‑1 (MQ=255)
                                                                               aTTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTc                                                                   <  1:2177460/149‑1 (MQ=255)
                                                                                ttAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCg                                                                  >  1:1268854/1‑149 (MQ=255)
                                                                                   aGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCt                                                               <  1:1732840/149‑1 (MQ=255)
                                                                                        cgTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTaaaaa                                                          <  1:562165/149‑1 (MQ=255)
                                                                                         gTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAAt                                                         >  2:2474315/1‑149 (MQ=255)
                                                                                           ggcggcTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGc                                                       <  1:1507096/149‑1 (MQ=255)
                                                                                              ggcTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCag                                                    <  2:4159865/149‑1 (MQ=255)
                                                                                              ggcTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCag                                                    >  1:5199027/1‑149 (MQ=255)
                                                                                               gcTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCaga                                                   >  2:4330166/1‑149 (MQ=255)
                                                                                                   aCGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCa                                               >  1:4916322/1‑149 (MQ=255)
                                                                                                        tCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGc                                          >  2:5081856/1‑149 (MQ=255)
                                                                                                            tAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAAt                                      <  1:2670599/149‑1 (MQ=255)
                                                                                                              accacaccacACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTa                                    <  2:1636100/149‑1 (MQ=255)
                                                                                                              accacaccacACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTa                                    >  1:4339016/1‑149 (MQ=255)
                                                                                                               ccacaccacaCGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTAt                                   >  2:4129185/1‑149 (MQ=255)
                                                                                                               ccacaccacaCGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTAt                                   <  1:330431/149‑1 (MQ=255)
                                                                                                                cacaccacacGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATc                                  >  2:679489/1‑149 (MQ=255)
                                                                                                                  caccacacGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATcgg                                <  2:5212098/149‑1 (MQ=255)
                                                                                                                    ccacacGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATcggcg                              >  2:1397028/1‑149 (MQ=255)
                                                                                                                     cacacGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATcggcgg                             <  2:2797762/149‑1 (MQ=255)
                                                                                                                      acacGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCggcggc                            >  2:57832/1‑149 (MQ=255)
                                                                                                                        acGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGa                          >  1:218372/1‑149 (MQ=255)
                                                                                                                         cgcgAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGAc                         <  2:3702948/149‑1 (MQ=255)
                                                                                                                         cgcgAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGAc                         <  2:541155/149‑1 (MQ=255)
                                                                                                                                 cAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTc                 <  2:1985277/149‑1 (MQ=255)
                                                                                                                                  aGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTcc                <  2:1062704/149‑1 (MQ=255)
                                                                                                                                    ccATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCg              >  1:400619/1‑149 (MQ=255)
                                                                                                                                     cATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGc             >  1:4643580/1‑149 (MQ=255)
                                                                                                                                        tGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTc          <  2:85551/149‑1 (MQ=255)
                                                                                                                                        tGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTc          <  2:2608778/149‑1 (MQ=255)
                                                                                                                                          tCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAg        >  2:3890186/1‑149 (MQ=255)
                                                                                                                                           ccGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAgg       >  1:4510346/1‑149 (MQ=255)
                                                                                                                                           ccGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAgg       <  1:3825278/149‑1 (MQ=255)
                                                                                                                                             ggCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTg     >  1:2413581/1‑149 (MQ=255)
                                                                                                                                                gTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGa  >  2:3290123/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGA  >  CP009273/158141‑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: ‑