Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,223,247 C→T 92.5% 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,2470CT92.5% 621.5 / 33.6 200W206* (TGG→TGAdusCtRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand):  ref base C (7/8);  new base T (76/109);  total (83/117)
Fisher's exact test for biased strand distribution p-value = 7.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

Only 100 of 200 total aligned reads displayed.

TTGCAGCAAAGCAACCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGTAGCGCCAGCCTGT  >  CP009273/2223103‑2223391
                                                                                                                                                |                                                                                                                                                
ttGCAGCAAAGCAACCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAg                                                                                                                                               >  2:610139/1‑148 (MQ=255)
   cagcaAAGCAACCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTcc                                                                                                                                            <  1:2750732/148‑1 (MQ=255)
    agcaaagcaaCCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTccc                                                                                                                                           <  2:961952/148‑1 (MQ=255)
    agcaaagcaaCCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTccc                                                                                                                                           >  2:1512681/1‑148 (MQ=255)
     gcaaagcaaCCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCa                                                                                                                                          >  2:2312857/1‑148 (MQ=255)
        aagcaaCCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGAt                                                                                                                                       >  2:2234075/1‑148 (MQ=255)
        aagcaaCCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGAt                                                                                                                                       >  1:2089977/1‑148 (MQ=255)
         agcaaCCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGAtt                                                                                                                                      >  1:303952/1‑148 (MQ=255)
          gcaaCCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGAttt                                                                                                                                     <  1:1589636/148‑1 (MQ=255)
            aaCCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCa                                                                                                                                   <  1:867785/148‑1 (MQ=255)
             accaccTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCAc                                                                                                                                  <  1:484320/148‑1 (MQ=255)
              ccaccTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCAcc                                                                                                                                 >  2:3211192/1‑148 (MQ=255)
              ccaccTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCCGTCCCAGATTTCAcc                                                                                                                                 >  2:3378922/1‑148 (MQ=255)
               caccTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCg                                                                                                                                <  2:1198875/148‑1 (MQ=255)
                accTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCCTTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGt                                                                                                                               >  1:2991301/1‑148 (MQ=255)
                accTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGt                                                                                                                               >  1:331753/1‑148 (MQ=255)
                 ccTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGtt                                                                                                                              <  1:2583308/148‑1 (MQ=255)
                 ccTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGtt                                                                                                                              <  2:3116966/148‑1 (MQ=255)
                          aCGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCAc                                                                                                                     >  2:1040454/1‑148 (MQ=255)
                           cGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCAcc                                                                                                                    >  1:635728/1‑148 (MQ=255)
                           cGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCAcc                                                                                                                    >  1:3413979/1‑148 (MQ=255)
                              cATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGa                                                                                                                 <  2:1332071/148‑1 (MQ=255)
                                 tCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAata                                                                                                              >  2:1085433/1‑148 (MQ=255)
                                 tCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAata                                                                                                              <  1:1233728/148‑1 (MQ=255)
                                      gTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAg                                                                                                         <  1:1389764/148‑1 (MQ=255)
                                       ttCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGc                                                                                                        <  1:1953080/148‑1 (MQ=255)
                                        tCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGcc                                                                                                       >  1:2212156/1‑148 (MQ=255)
                                         cGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCg                                                                                                      <  2:1576140/148‑1 (MQ=255)
                                         cGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCg                                                                                                      <  2:1576628/148‑1 (MQ=255)
                                               atTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCg                                                                                                <  2:447588/148‑1 (MQ=255)
                                                ttttACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGa                                                                                               >  1:986900/1‑148 (MQ=255)
                                                  ttACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAat                                                                                             <  1:52371/148‑1 (MQ=255)
                                                   tACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAata                                                                                            <  1:3524241/148‑1 (MQ=255)
                                                    accaccCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAatat                                                                                           >  1:652433/1‑148 (MQ=255)
                                                    accaccCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAatat                                                                                           >  1:925497/1‑148 (MQ=255)
                                                     ccaccCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATc                                                                                          >  2:2893725/1‑148 (MQ=255)
                                                      caccCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCg                                                                                         >  2:879034/1‑148 (MQ=255)
                                                         ccGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCa                                                                                      <  2:3530541/148‑1 (MQ=255)
                                                           ggCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAAt                                                                                    <  1:1242274/148‑1 (MQ=255)
                                                           ggCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAAt                                                                                    <  1:3170107/148‑1 (MQ=255)
                                                           ggCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAAt                                                                                    <  2:844907/148‑1 (MQ=255)
                                                            gCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATc                                                                                   >  1:1031678/1‑148 (MQ=255)
                                                            gCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATAGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATATCACCGTTGGCAATCACCGGAATATACAGCCGCTGGCGAATATCGCCAATc                                                                                   >  2:3312605/1‑148 (MQ=255)
                                                             cTCAGGTTGGGAATATAGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCg                                                                                  >  2:632477/1‑148 (MQ=255)
                                                              tCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGc                                                                                 <  2:2847475/148‑1 (MQ=255)
                                                                 ggTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTg                                                                              >  2:2430875/1‑148 (MQ=255)
                                                                 ggTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTg                                                                              >  2:1397950/1‑148 (MQ=255)
                                                                  gTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGc                                                                             <  2:548736/148‑1 (MQ=255)
                                                                  gTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGc                                                                             >  2:806981/1‑148 (MQ=255)
                                                                  gTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGc                                                                             >  1:2935987/1‑148 (MQ=255)
                                                                  gTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGc                                                                             <  1:1699434/148‑1 (MQ=255)
                                                                    tGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCa                                                                           <  2:604718/148‑1 (MQ=255)
                                                                    tGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCa                                                                           <  2:438304/148‑1 (MQ=255)
                                                                       gAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTc                                                                        >  1:1651394/1‑148 (MQ=255)
                                                                        aaTATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCa                                                                       <  2:303952/148‑1 (MQ=255)
                                                                        aaTATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCa                                                                       <  2:2101212/148‑1 (MQ=255)
                                                                          tatTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAat                                                                     <  1:4966/148‑1 (MQ=255)
                                                                          tatTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAat                                                                     >  1:1983143/1‑148 (MQ=255)
                                                                            ttGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAatat                                                                   >  2:1617582/1‑148 (MQ=255)
                                                                             tGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATg                                                                  <  2:648692/148‑1 (MQ=255)
                                                                              gAGCGCCCCGCGACCAATCATCCCTGCGTCGCAGCCGCTGATCGCCAGGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGc                                                                 <  2:2278453/148‑1 (MQ=255)
                                                                                 cgcCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTcc                                                              >  1:2439228/1‑148 (MQ=255)
                                                                                  gcCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCcg                                                             <  2:2922182/148‑1 (MQ=255)
                                                                                     ccGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCcgcgc                                                          >  2:3345126/1‑148 (MQ=255)
                                                                                      cgcgACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCcgcgcg                                                         >  1:819686/1‑148 (MQ=255)
                                                                                             aaTCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAAccc                                                  <  1:1674497/148‑1 (MQ=255)
                                                                                               tcatcaCTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTg                                                <  2:2038858/148‑1 (MQ=255)
                                                                                               tcatcaCTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTg                                                <  2:2038296/148‑1 (MQ=255)
                                                                                               tcatcaCTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTg                                                <  2:769642/148‑1 (MQ=255)
                                                                                               tcatcaCTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTg                                                <  1:2854413/148‑1 (MQ=255)
                                                                                                  tcaCTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGctc                                             <  1:1824284/148‑1 (MQ=255)
                                                                                                  tcaCTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGctc                                             <  1:1115605/148‑1 (MQ=255)
                                                                                                  tcaCTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGctc                                             >  2:1864174/1‑148 (MQ=255)
                                                                                                  tcaCTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGctc                                             >  1:3377767/1‑148 (MQ=255)
                                                                                                   caCTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGctct                                            <  2:2724166/148‑1 (MQ=255)
                                                                                                   caCTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGctct                                            <  1:3024422/148‑1 (MQ=255)
                                                                                                   caCTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGctct                                            <  2:2530991/148‑1 (MQ=255)
                                                                                                    tcTGCGTCGCAGCCGCTGATCTCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCTCTGGAGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACACTGCTCtt                                           <  2:1562786/147‑1 (MQ=255)
                                                                                                    aCTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCtt                                           <  1:279095/148‑1 (MQ=255)
                                                                                                    aCTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCtt                                           <  1:907760/148‑1 (MQ=255)
                                                                                                        cGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGt                                       <  1:879034/148‑1 (MQ=255)
                                                                                                             cAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGcc                                  <  2:2620545/148‑1 (MQ=255)
                                                                                                               gCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCa                                >  2:2833750/1‑148 (MQ=255)
                                                                                                               gCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCa                                <  2:3026572/148‑1 (MQ=255)
                                                                                                               gCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCa                                <  2:961483/148‑1 (MQ=255)
                                                                                                                     gATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGcacc                          >  1:712973/1‑148 (MQ=255)
                                                                                                                     gATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGcacc                          <  1:1512681/148‑1 (MQ=255)
                                                                                                                        cGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCaccacc                       <  1:522153/148‑1 (MQ=255)
                                                                                                                           cATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGc                    <  1:2312864/148‑1 (MQ=255)
                                                                                                                               cATTGTTGCGCGCTCTGTCAGTCCCTGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCg                <  1:2621074/148‑1 (MQ=255)
                                                                                                                                aTTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGt               <  1:454607/148‑1 (MQ=255)
                                                                                                                                aTTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGt               <  1:3211192/148‑1 (MQ=255)
                                                                                                                                 ttgttgCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGTa              >  2:3192110/1‑148 (MQ=255)
                                                                                                                                    ttgCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGTAgcg           <  2:510652/148‑1 (MQ=255)
                                                                                                                                     tgCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGTAgcgc          <  1:1020708/148‑1 (MQ=255)
                                                                                                                                     tgCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGTAgcgc          >  2:2681453/1‑148 (MQ=255)
                                                                                                                                       cgcgcTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGTAGCGCCa        <  2:677663/148‑1 (MQ=255)
                                                                                                                                         cgcTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGTAGCGCCAGc      <  1:2931011/148‑1 (MQ=255)
                                                                                                                                          gcTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGTAGCGCCAGcc     <  1:908471/148‑1 (MQ=255)
                                                                                                                                             ctGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGTAGCGCCAGCCTGt  >  2:737417/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
TTGCAGCAAAGCAACCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGTAGCGCCAGCCTGT  >  CP009273/2223103‑2223391

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: