Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,223,247 C→T 41.8% W206* (TGG→TGA dusC ← tRNA‑dihydrouridine synthase C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,223,2470CT41.8% 223.8 / inf 328W206* (TGG→TGAdusCtRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand):  ref base C (98/92);  new base T (79/58);  total (178/150)
Fisher's exact test for biased strand distribution p-value = 3.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 331 total aligned reads displayed.

TATTTTTGCAGCAAAGCAACCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGTAGCGCCAGCCTGTTGAAC  >  CP009273/2223098‑2223396
                                                                                                                                                     |                                                                                                                                                     
 atTTTTGCAGCAAAGCAACCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGcc                                                                                                                                                      >  1:4790107/1‑150 (MQ=255)
   ttttGCAGCAAAGCAACCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAg                                                                                                                                                    <  2:3579945/150‑1 (MQ=255)
    tttGCAGCAAAGCAACCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGt                                                                                                                                                   <  1:5087779/150‑1 (MQ=255)
     ttGCAGCAAAGCAACCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTc                                                                                                                                                  <  2:3855500/150‑1 (MQ=255)
      tGCAGCAAAGCAACCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTcc                                                                                                                                                 <  2:2692021/150‑1 (MQ=255)
      tGCAGCAAAGCAACCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTcc                                                                                                                                                 <  2:3751823/150‑1 (MQ=255)
      tGCAGCAAAGCAACCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTcc                                                                                                                                                 <  1:5055155/150‑1 (MQ=255)
          gcaaagcaaCCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGa                                                                                                                                             >  1:414054/1‑150 (MQ=255)
          gcaaagcaaCCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGa                                                                                                                                             <  1:3745298/150‑1 (MQ=255)
          gcaaagcaaCCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGa                                                                                                                                             >  1:1044804/1‑150 (MQ=255)
           caaagcaaCCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGAt                                                                                                                                            >  1:3980138/1‑150 (MQ=255)
            aaagcaaCCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGAtt                                                                                                                                           <  1:3365997/150‑1 (MQ=255)
             aagcaaACACCTCCTGCCACGGCTTTCGCTGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCTCTCTGTCAGTCCCAGAttt                                                                                                                                          <  2:3692508/150‑1 (MQ=255)
               gcaaCCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCa                                                                                                                                        >  2:251552/1‑150 (MQ=255)
                caaCCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCAc                                                                                                                                       <  1:3057110/150‑1 (MQ=255)
                   ccaccTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGt                                                                                                                                    >  2:1873547/1‑150 (MQ=255)
                       cTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGc                                                                                                                                >  1:662622/1‑150 (MQ=255)
                         ccGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCaa                                                                                                                              <  2:3320254/150‑1 (MQ=255)
                           ggCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATc                                                                                                                            <  1:5009356/150‑1 (MQ=255)
                            gCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGACTATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCAGTGTTGCGCGCTCTGCCAGTCCCAGAGTTCACCGGTGGCACTcc                                                                                                                           >  2:3173453/1‑149 (MQ=255)
                                 ggCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGa                                                                                                                      <  2:1093013/150‑1 (MQ=255)
                                 ggCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGa                                                                                                                      >  1:4502611/1‑150 (MQ=255)
                                   cATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAat                                                                                                                    >  2:4999168/1‑150 (MQ=255)
                                      tCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATAtt                                                                                                                 >  1:2061491/1‑150 (MQ=255)
                                       cgcgGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTc                                                                                                                >  2:3430048/1‑150 (MQ=255)
                                        gcgGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCa                                                                                                               <  2:4277662/150‑1 (MQ=255)
                                          ggTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGc                                                                                                             >  1:2592620/1‑150 (MQ=255)
                                           gTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGcc                                                                                                            <  1:3707374/150‑1 (MQ=255)
                                               gTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTg                                                                                                        <  1:2207532/150‑1 (MQ=255)
                                               gTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTg                                                                                                        <  1:719404/150‑1 (MQ=255)
                                                 tataTTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGc                                                                                                      >  1:4392986/1‑150 (MQ=255)
                                                 tataTTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGc                                                                                                      >  1:1595321/1‑150 (MQ=255)
                                                  atatTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCg                                                                                                     >  1:5065803/1‑150 (MQ=255)
                                                  atatTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCg                                                                                                     >  1:2097551/1‑150 (MQ=255)
                                                  atatTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCg                                                                                                     >  1:480766/1‑150 (MQ=255)
                                                   tatTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTTCGACGCTGCCGCTCATCGCCAGGCCTTGATGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCCATCACCGGACTAGTCAGCCGCTGGCGa                                                                                                    >  1:2529198/1‑150 (MQ=255)
                                                    atTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGaa                                                                                                   <  2:20403/150‑1 (MQ=255)
                                                       ttACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAatat                                                                                                >  2:395113/1‑150 (MQ=255)
                                                        tACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATc                                                                                               >  1:825239/1‑150 (MQ=255)
                                                        tACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATc                                                                                               >  2:585678/1‑150 (MQ=255)
                                                        tACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATGTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATc                                                                                               >  2:3270049/1‑150 (MQ=255)
                                                           caccCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGcc                                                                                            >  2:793829/1‑150 (MQ=255)
                                                            accCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCa                                                                                           >  2:892149/1‑150 (MQ=255)
                                                              ccGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGTCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAAt                                                                                         <  1:3460845/150‑1 (MQ=255)
                                                               cGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCTCCGGAATATTCAGCCGCTGGCTAATATCGCCAATc                                                                                        <  1:273591/150‑1 (MQ=255)
                                                                   tCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCCTCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCNt                                                                                    <  1:4945991/150‑1 (MQ=255)
                                                                    cAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTg                                                                                   >  1:1450418/1‑150 (MQ=255)
                                                                    cAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTg                                                                                   <  1:4245952/150‑1 (MQ=255)
                                                                     aGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGc                                                                                  >  2:3461965/1‑150 (MQ=255)
                                                                       gTTGGGAATATTGAGCGCCCCGCGACCACTCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCa                                                                                <  2:1946611/150‑1 (MQ=255)
                                                                        ttGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAg                                                                               >  1:1123075/1‑150 (MQ=255)
                                                                        ttGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAg                                                                               >  2:615152/1‑150 (MQ=255)
                                                                            gAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCaa                                                                           >  2:378888/1‑150 (MQ=255)
                                                                             aaTATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAat                                                                          >  2:4004730/1‑150 (MQ=255)
                                                                              atatTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAata                                                                         <  2:4202908/150‑1 (MQ=255)
                                                                                atTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATg                                                                       >  2:3654027/1‑150 (MQ=255)
                                                                                atTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATg                                                                       >  2:1561623/1‑150 (MQ=255)
                                                                                atTGAGCGCCCCGAGACCAATCAACAATGAGTCGCAGCCTCTGATCGCCATGCATTGTTTCGCGCTCTGACAGTCCCATATTACACAGTTGGCAATCACCGGAATATTAAGCCGCTGGCGAATATCGCAAATCTCCTGCCAGTCAATATg                                                                       >  2:285314/1‑150 (MQ=255)
                                                                                    aGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTcc                                                                   <  2:830851/150‑1 (MQ=255)
                                                                                     gcgcCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCcg                                                                  >  2:179506/1‑150 (MQ=255)
                                                                                     gcgcCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCcg                                                                  <  1:4429468/150‑1 (MQ=255)
                                                                                     gcgcCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGCCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCcg                                                                  <  1:1722103/150‑1 (MQ=255)
                                                                                        ccccGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCcgcgc                                                               <  2:2270737/150‑1 (MQ=255)
                                                                                         cccGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCcgcgcg                                                              >  2:2134836/1‑150 (MQ=255)
                                                                                          ccGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCgg                                                             <  2:4135041/150‑1 (MQ=255)
                                                                                           cgcgACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGt                                                            >  1:3288498/1‑150 (MQ=255)
                                                                                           cgcgACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGt                                                            >  1:2378592/1‑150 (MQ=255)
                                                                                            gcgACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTa                                                           <  1:2914869/150‑1 (MQ=255)
                                                                                             cgACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTaa                                                          <  2:4823297/150‑1 (MQ=255)
                                                                                              gACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAAc                                                         <  1:1383754/150‑1 (MQ=255)
                                                                                                ccAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAAccc                                                       <  2:4385267/150‑1 (MQ=255)
                                                                                                  aaTCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTg                                                     <  1:1088864/150‑1 (MQ=255)
                                                                                                    tcatcaCTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGct                                                   <  2:3020614/150‑1 (MQ=255)
                                                                                                      atcaCTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGctct                                                 <  2:1346350/150‑1 (MQ=255)
                                                                                                        caCTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCttt                                               <  1:1526920/150‑1 (MQ=255)
                                                                                                        caCTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCAAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCttt                                               >  2:2464960/1‑150 (MQ=255)
                                                                                                             cGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTcc                                          <  1:5027861/150‑1 (MQ=255)
                                                                                                              gTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCg                                         <  2:296978/150‑1 (MQ=255)
                                                                                                              gTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCACTATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCg                                         >  2:4212325/1‑150 (MQ=255)
                                                                                                                    gCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATg                                   <  2:2577330/150‑1 (MQ=255)
                                                                                                                      cGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGTTTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAAACCTGCTCTTTCGTCCGCCCATGca                                 >  2:4719020/1‑150 (MQ=255)
                                                                                                                      cGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGca                                 >  1:4847902/1‑150 (MQ=255)
                                                                                                                        cTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGcacc                               <  2:1164238/150‑1 (MQ=255)
                                                                                                                               ccATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCt                        >  1:3583055/1‑150 (MQ=255)
                                                                                                                               ccATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCt                        >  1:854180/1‑150 (MQ=255)
                                                                                                                               ccATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCt                        <  2:4743925/150‑1 (MQ=255)
                                                                                                                                cATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTc                       >  1:630887/1‑150 (MQ=255)
                                                                                                                                  tGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCg                     <  1:2937346/150‑1 (MQ=255)
                                                                                                                                   gCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGt                    >  1:1174481/1‑150 (MQ=255)
                                                                                                                                   gCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGt                    >  2:3235101/1‑150 (MQ=255)
                                                                                                                                   gCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGt                    <  1:3691385/150‑1 (MQ=255)
                                                                                                                                    cATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGTa                   >  1:2333584/1‑150 (MQ=255)
                                                                                                                                       tgttgCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGTAgcg                >  1:543851/1‑150 (MQ=255)
                                                                                                                                        gttgCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGTAgcgc               >  1:4471960/1‑150 (MQ=255)
                                                                                                                                        gttgCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGTAgcgc               >  1:1919346/1‑150 (MQ=255)
                                                                                                                                              cgcTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGTAGCGCCAGCCt         <  2:5065761/150‑1 (MQ=255)
                                                                                                                                              cgcTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGTAGCGCCAGCCt         >  1:2813440/1‑150 (MQ=255)
                                                                                                                                                ctctGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATTCACCACCAGATCCGTAGCGCCAGCCTGt       <  1:911885/150‑1 (MQ=255)
                                                                                                                                                 tctGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGTAGCGCCAGCCTGtt      <  2:662622/150‑1 (MQ=255)
                                                                                                                                                 tctGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGTAGCGCCAGCCTGtt      <  2:1716123/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TATTTTTGCAGCAAAGCAACCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGTAGCGCCAGCCTGTTGAAC  >  CP009273/2223098‑2223396

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: