Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,223,247 C→T 100% 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,2470CT97.8% 734.3 / 5.0 223W206* (TGG→TGA) dusCtRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base T (86/132);  total (89/134)
Fisher's exact test for biased strand distribution p-value = 3.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 223 total aligned reads displayed.

CAGCAAAGCAACCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGTAGCGCCAGCCTGT  >  CP009273/2223106‑2223391
                                                                                                                                             |                                                                                                                                                
cagcaAAGCAACCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTcc                                                                                                                                            <  1:1558426/148‑1 (MQ=255)
   caaagcaaCCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAg                                                                                                                                         >  1:1585656/1‑148 (MQ=255)
     aagcaaCCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGAt                                                                                                                                       >  1:3077666/1‑148 (MQ=255)
     aagcaaCCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGAt                                                                                                                                       >  1:2272439/1‑148 (MQ=255)
     aagcaaCCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGAt                                                                                                                                       >  1:3535783/1‑148 (MQ=255)
     aagcaaCCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGAt                                                                                                                                       <  2:1565162/148‑1 (MQ=255)
       gcaaCCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGAttt                                                                                                                                     <  1:2212075/148‑1 (MQ=255)
       gcaaCCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGAttt                                                                                                                                     >  1:604410/1‑148 (MQ=255)
          accaccTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCAc                                                                                                                                  <  2:3323289/148‑1 (MQ=255)
           ccaccTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCAcc                                                                                                                                 >  2:621634/1‑148 (MQ=255)
           ccaccTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCAcc                                                                                                                                 >  1:1616940/1‑148 (MQ=255)
             accTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGt                                                                                                                               <  2:801906/148‑1 (MQ=255)
              ccTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGtt                                                                                                                              <  1:2801275/148‑1 (MQ=255)
              ccTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGtt                                                                                                                              >  1:2336186/1‑148 (MQ=255)
              ccTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGtt                                                                                                                              >  1:2481356/1‑148 (MQ=255)
                tCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTgg                                                                                                                            <  1:1647054/148‑1 (MQ=255)
                   ggCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCaa                                                                                                                         <  2:858443/148‑1 (MQ=255)
                    gCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAAt                                                                                                                        >  1:3226066/1‑148 (MQ=255)
                     ccACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATc                                                                                                                       >  1:1931004/1‑148 (MQ=255)
                      cACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCa                                                                                                                      <  1:2136090/148‑1 (MQ=255)
                       aCGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCAc                                                                                                                     <  1:3122148/148‑1 (MQ=255)
                       aCGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAACGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCAc                                                                                                                     >  2:1799612/1‑148 (MQ=255)
                        cGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCAcc                                                                                                                    <  1:3110996/148‑1 (MQ=255)
                        cGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCAcc                                                                                                                    <  1:649079/148‑1 (MQ=255)
                              tCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAata                                                                                                              >  2:243527/1‑148 (MQ=255)
                              tCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAata                                                                                                              >  1:1989959/1‑148 (MQ=255)
                                      cGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCg                                                                                                      <  1:3388218/148‑1 (MQ=255)
                                      cGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCg                                                                                                      >  2:1773933/1‑148 (MQ=255)
                                      cGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCg                                                                                                      <  1:3261240/148‑1 (MQ=255)
                                      cGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCg                                                                                                      >  1:64052/1‑148 (MQ=255)
                                      cGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCg                                                                                                      <  1:1422618/148‑1 (MQ=255)
                                      cGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCg                                                                                                      <  1:194556/148‑1 (MQ=255)
                                       gTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGc                                                                                                     <  2:1280078/148‑1 (MQ=255)
                                        ttATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCt                                                                                                    >  1:2925485/1‑148 (MQ=255)
                                            atTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGACTATTCAGCCGCTGGCg                                                                                                >  2:1881363/1‑148 (MQ=255)
                                             ttttACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGa                                                                                               <  1:2058177/148‑1 (MQ=255)
                                              tttACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGaa                                                                                              >  2:510162/1‑148 (MQ=255)
                                               ttACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAat                                                                                             <  2:1907216/148‑1 (MQ=255)
                                               ttACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAat                                                                                             <  2:3291343/148‑1 (MQ=255)
                                                tACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAata                                                                                            <  1:1619864/148‑1 (MQ=255)
                                                 accaccCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAatat                                                                                           >  2:2882396/1‑148 (MQ=255)
                                                 accaccCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAatat                                                                                           >  2:25574/1‑148 (MQ=255)
                                                 accaccCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAatat                                                                                           <  1:3491977/148‑1 (MQ=255)
                                                    accCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGc                                                                                        <  2:498623/148‑1 (MQ=255)
                                                    accCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGc                                                                                        <  2:245232/148‑1 (MQ=255)
                                                     cccGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGcc                                                                                       >  2:1320517/1‑148 (MQ=255)
                                                          cTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCg                                                                                  >  1:1679202/1‑148 (MQ=255)
                                                            cAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGcc                                                                                <  1:2891664/148‑1 (MQ=255)
                                                              ggTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTg                                                                              >  2:2760566/1‑148 (MQ=255)
                                                               gTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGc                                                                             <  1:3054042/148‑1 (MQ=255)
                                                                ttGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGcc                                                                            >  1:2407943/1‑148 (MQ=255)
                                                                   ggAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGt                                                                         >  2:2456725/1‑148 (MQ=255)
                                                                    gAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTc                                                                        <  2:14596/148‑1 (MQ=255)
                                                                    gAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTc                                                                        <  2:3082415/148‑1 (MQ=255)
                                                                      atatTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCaa                                                                      >  1:3552414/1‑148 (MQ=255)
                                                                       tatTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAat                                                                     >  1:3334015/1‑148 (MQ=255)
                                                                       tatTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAat                                                                     <  2:3196812/148‑1 (MQ=255)
                                                                       tatTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAat                                                                     >  1:2569390/1‑148 (MQ=255)
                                                                        atTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAata                                                                    >  1:1223388/1‑148 (MQ=255)
                                                                          tGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATg                                                                  <  2:1167628/148‑1 (MQ=255)
                                                                            aGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCt                                                                <  2:1660883/148‑1 (MQ=255)
                                                                                     cgACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGt                                                       <  2:405098/148‑1 (MQ=255)
                                                                                        ccAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCTGTAAc                                                    >  1:2024259/1‑148 (MQ=255)
                                                                                        ccAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAAc                                                    >  2:663030/1‑148 (MQ=255)
                                                                                         cAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAAcc                                                   <  2:1707342/148‑1 (MQ=255)
                                                                                          aaTCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAAccc                                                  <  1:1782159/148‑1 (MQ=255)
                                                                                            tcatcaCTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTg                                                >  2:1894105/1‑148 (MQ=255)
                                                                                              atcaCTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGct                                              <  2:1032702/148‑1 (MQ=255)
                                                                                              atcaCTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGct                                              <  1:994556/148‑1 (MQ=255)
                                                                                              atcaCTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGct                                              <  2:2807965/148‑1 (MQ=255)
                                                                                                caCTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGctct                                            <  1:187566/148‑1 (MQ=255)
                                                                                                 aCTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCtt                                           <  1:2771420/148‑1 (MQ=255)
                                                                                                  cTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTACCCCTGCTCttt                                          <  1:794529/148‑1 (MQ=255)
                                                                                                  cTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCttt                                          <  2:197200/148‑1 (MQ=255)
                                                                                                  cTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCttt                                          <  1:3062034/148‑1 (MQ=255)
                                                                                                     cGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGt                                       >  1:396736/1‑148 (MQ=255)
                                                                                                        cGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCg                                    >  2:54503/1‑148 (MQ=255)
                                                                                                          cAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGcc                                  <  1:501588/148‑1 (MQ=255)
                                                                                                          cAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGcc                                  <  1:1042706/148‑1 (MQ=255)
                                                                                                          cAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGcc                                  >  1:3558294/1‑148 (MQ=255)
                                                                                                            gCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCa                                >  2:3434889/1‑148 (MQ=255)
                                                                                                            gCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCa                                >  2:873850/1‑148 (MQ=255)
                                                                                                            gCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCa                                >  2:873807/1‑148 (MQ=255)
                                                                                                                   aTCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGcacca                         <  2:1588596/148‑1 (MQ=255)
                                                                                                                   aTCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGcacca                         <  2:1045035/148‑1 (MQ=255)
                                                                                                                     cGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCaccacc                       <  1:2908641/148‑1 (MQ=255)
                                                                                                                       ccATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAg                     <  1:585104/148‑1 (MQ=255)
                                                                                                                        cATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGc                    <  1:1947834/148‑1 (MQ=255)
                                                                                                                          tGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTc                  >  2:3550438/1‑148 (MQ=255)
                                                                                                                          tGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTc                  >  2:2395220/1‑148 (MQ=255)
                                                                                                                            cATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCg                <  1:1278966/148‑1 (MQ=255)
                                                                                                                              ttgttgCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGTa              >  1:75694/1‑148 (MQ=255)
                                                                                                                              ttgttgCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGTa              >  1:1091446/1‑148 (MQ=255)
                                                                                                                              ttgttgCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCAGGCACCACCAGCTCCGTa              >  1:2865960/1‑148 (MQ=255)
                                                                                                                                  tgCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGTAgcgc          <  1:2995694/148‑1 (MQ=255)
                                                                                                                                   gcgcgcTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGTAGCGcc         <  1:2990573/148‑1 (MQ=255)
                                                                                                                                         tctGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGTAGCGCCAGCCTg   <  2:3245445/148‑1 (MQ=255)
                                                                                                                                          ctGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGTAGCGCCAGCCTGt  >  2:3538301/1‑148 (MQ=255)
                                                                                                                                          ctGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGTAGCGCCAGCCTGt  >  1:90332/1‑148 (MQ=255)
                                                                                                                                          ctGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGTAGCGCCAGCCTGt  <  1:2519571/148‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                
CAGCAAAGCAACCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGTAGCGCCAGCCTGT  >  CP009273/2223106‑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: ‑