Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 158,287 C→A 79.4% 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,2870CA79.4% 759.8 / 247.8 349D78Y (GAC→TAC) ligT2'‑5' RNA ligase
Reads supporting (aligned to +/- strand):  ref base C (37/35);  new base A (128/149);  total (165/184)
Fisher's exact test for biased strand distribution p-value = 5.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 360 total aligned reads displayed.

GTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGG  >  CP009273/158139‑158431
                                                                                                                                                    |                                                                                                                                                
gTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGtc                                                                                                                                                  >  2:3115912/1‑149 (MQ=255)
    tatGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGtaatca                                                                                                                                              <  2:2486639/149‑1 (MQ=255)
         ggATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGt                                                                                                                                         >  1:1352763/1‑149 (MQ=255)
          gATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTg                                                                                                                                        >  1:1259820/1‑149 (MQ=255)
           aTGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGa                                                                                                                                       >  2:401347/1‑149 (MQ=255)
            tGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAg                                                                                                                                      >  1:1705629/1‑149 (MQ=255)
             gAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAgt                                                                                                                                     <  2:3895971/149‑1 (MQ=255)
             gAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAgt                                                                                                                                     >  2:600219/1‑149 (MQ=255)
              aaaCGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAgtg                                                                                                                                    <  1:3020958/149‑1 (MQ=255)
               aaCGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAgtgt                                                                                                                                   >  1:1304246/1‑149 (MQ=255)
                 cGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGa                                                                                                                                 >  2:4485536/1‑149 (MQ=255)
                     cGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAAcc                                                                                                                             <  1:2071920/149‑1 (MQ=255)
                      gATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGACACCa                                                                                                                            >  2:4560102/1‑149 (MQ=255)
                       aTTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAg                                                                                                                           >  2:4322383/1‑149 (MQ=255)
                       aTTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAg                                                                                                                           >  2:825206/1‑149 (MQ=255)
                        ttGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAgg                                                                                                                          >  2:1256295/1‑149 (MQ=255)
                         tGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGt                                                                                                                         >  2:1171402/1‑149 (MQ=255)
                           cTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTg                                                                                                                       >  2:269772/1‑149 (MQ=255)
                             ttGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGAc                                                                                                                     <  1:4530702/149‑1 (MQ=255)
                             ttGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGAc                                                                                                                     <  2:103983/149‑1 (MQ=255)
                             ttGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGAc                                                                                                                     <  1:345198/149‑1 (MQ=255)
                               gAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGa                                                                                                                   >  1:4369567/1‑149 (MQ=255)
                               gAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGa                                                                                                                   >  2:1518231/1‑149 (MQ=255)
                               gAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGa                                                                                                                   <  2:2539373/149‑1 (MQ=255)
                                aaaaCAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGaa                                                                                                                  >  1:885231/1‑149 (MQ=255)
                                 aaaCAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAAt                                                                                                                 >  1:3159674/1‑149 (MQ=255)
                                  aacaacCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCCCACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAAtc                                                                                                                <  2:48241/149‑1 (MQ=255)
                                  aacaacCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAAtc                                                                                                                <  1:2572070/149‑1 (MQ=255)
                                  aacaacCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAAtc                                                                                                                <  1:1895560/149‑1 (MQ=255)
                                  aacaacCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAAtc                                                                                                                <  1:721860/149‑1 (MQ=255)
                                     aacCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAAtccgt                                                                                                             <  2:3513617/149‑1 (MQ=255)
                                      acCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAAtccgtc                                                                                                            <  1:190432/149‑1 (MQ=255)
                                      acCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAAtccgtc                                                                                                            >  1:632589/1‑149 (MQ=255)
                                       ccGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGATGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAAtccgtcc                                                                                                           >  1:3096097/1‑149 (MQ=255)
                                       ccGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAAtccgtcc                                                                                                           <  1:2646449/149‑1 (MQ=255)
                                       ccGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAAtccgtcc                                                                                                           <  2:2414460/149‑1 (MQ=255)
                                          cTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGc                                                                                                        <  2:4456935/149‑1 (MQ=255)
                                           tGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCt                                                                                                       <  2:3348549/149‑1 (MQ=255)
                                              gggCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTaaa                                                                                                    <  2:1963565/149‑1 (MQ=255)
                                               ggCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTaaaa                                                                                                   >  1:1638867/1‑149 (MQ=255)
                                                     cTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGaaa                                                                                             >  1:4211059/1‑149 (MQ=255)
                                                     cTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGaaa                                                                                             <  2:2057806/149‑1 (MQ=255)
                                                      tgtgAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAg                                                                                            >  1:2038075/1‑149 (MQ=255)
                                                      tgtgAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAg                                                                                            >  1:1187413/1‑149 (MQ=255)
                                                       gtgAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAgc                                                                                           >  2:845124/1‑149 (MQ=255)
                                                           aCGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCt                                                                                       <  1:3791026/149‑1 (MQ=255)
                                                             ggAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTtc                                                                                     >  1:227060/1‑149 (MQ=255)
                                                              gAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTtct                                                                                    >  1:4325011/1‑149 (MQ=255)
                                                               aGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTtctc                                                                                   <  1:3523526/149‑1 (MQ=255)
                                                               aGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTtctc                                                                                   <  1:2346553/149‑1 (MQ=255)
                                                                gCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTctct                                                                                  >  2:1466388/1‑149 (MQ=255)
                                                                gCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTctct                                                                                  >  1:574961/1‑149 (MQ=255)
                                                                  atatTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCttt                                                                                <  2:2610188/149‑1 (MQ=255)
                                                                   tatTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTtc                                                                               >  1:3842897/1‑149 (MQ=255)
                                                                     ttCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTtctc                                                                             <  2:2957354/149‑1 (MQ=255)
                                                                      tCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCg                                                                            <  1:3204785/149‑1 (MQ=255)
                                                                        gCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGgc                                                                          >  2:2971146/1‑149 (MQ=255)
                                                                            gCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTc                                                                      <  1:4670945/149‑1 (MQ=255)
                                                                            gCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTc                                                                      >  2:2776357/1‑149 (MQ=255)
                                                                               ggATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACt                                                                   >  1:1817702/1‑149 (MQ=255)
                                                                               ggATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACt                                                                   >  1:1509514/1‑149 (MQ=255)
                                                                                 aTTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACCCGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTc                                                                 <  1:1725290/149‑1 (MQ=255)
                                                                                 aTTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTc                                                                 <  1:4484086/149‑1 (MQ=255)
                                                                                    aaGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGcc                                                              <  1:2606551/149‑1 (MQ=255)
                                                                                     aGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCt                                                             <  1:139886/149‑1 (MQ=255)
                                                                                        cgcgTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTaaa                                                          <  1:4108806/149‑1 (MQ=255)
                                                                                         gcgTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTaaaa                                                         <  2:628129/149‑1 (MQ=255)
                                                                                           gTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAAt                                                       <  2:1705629/149‑1 (MQ=255)
                                                                                           gTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCTGCGCTCACTTCGCCTAAAAAt                                                       <  2:1269614/149‑1 (MQ=255)
                                                                                           gTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAAt                                                       <  2:4166378/149‑1 (MQ=255)
                                                                                                  cTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCagag                                                >  1:767498/1‑149 (MQ=255)
                                                                                                  cTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCagag                                                <  1:1035996/149‑1 (MQ=255)
                                                                                                  cTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCagag                                                <  2:1647691/149‑1 (MQ=255)
                                                                                                   tGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGt                                               <  1:4575748/149‑1 (MQ=255)
                                                                                                   tGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGt                                               >  1:2212443/1‑149 (MQ=255)
                                                                                                   tGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGt                                               >  1:4056929/1‑149 (MQ=255)
                                                                                                    gACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTc                                              >  2:1804939/1‑149 (MQ=255)
                                                                                                     aCGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCa                                             <  1:2791617/149‑1 (MQ=255)
                                                                                                     aCGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCa                                             >  2:3172223/1‑149 (MQ=255)
                                                                                                         aTCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATg                                         >  2:321316/1‑149 (MQ=255)
                                                                                                         aTCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATg                                         >  1:2563216/1‑149 (MQ=255)
                                                                                                               aaCCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAAtt                                   <  2:2273419/149‑1 (MQ=255)
                                                                                                               aaCCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAAtt                                   <  2:2820202/149‑1 (MQ=255)
                                                                                                               aaCCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAAtt                                   <  2:3745233/149‑1 (MQ=255)
                                                                                                               aaCCACACCACACGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAAtt                                   <  2:4390403/149‑1 (MQ=255)
                                                                                                                   acaccacacGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCTCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATcg                               <  2:3032537/149‑1 (MQ=255)
                                                                                                                     accacacGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATcggc                             <  2:1022495/149‑1 (MQ=255)
                                                                                                                          acGCGAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGa                        >  2:4641832/1‑149 (MQ=255)
                                                                                                                            gcgAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGAcc                      <  1:1013107/149‑1 (MQ=255)
                                                                                                                            gcgAACGCAGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGAcc                      >  2:2366693/1‑149 (MQ=255)
                                                                                                                                    aGCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTcc              <  2:3478034/149‑1 (MQ=255)
                                                                                                                                     gCCATTGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTccc             >  1:2481876/1‑149 (MQ=255)
                                                                                                                                       cATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGCCGTCCCGc           <  1:2722314/149‑1 (MQ=255)
                                                                                                                                         ttGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCt         >  1:2789450/1‑149 (MQ=255)
                                                                                                                                          tGTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTc        <  1:37824/149‑1 (MQ=255)
                                                                                                                                           gTCCGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCa       >  1:2871512/1‑149 (MQ=255)
                                                                                                                                             ccGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAgg     <  2:2883127/149‑1 (MQ=255)
                                                                                                                                             ccGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAgg     >  2:3925131/1‑149 (MQ=255)
                                                                                                                                             ccGGCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAgg     >  2:2581546/1‑149 (MQ=255)
                                                                                                                                               ggCGTAATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTg   >  1:3972546/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGG  >  CP009273/158139‑158431

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: