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,2470CT99.5% 655.1 / ‑4.4 198W206* (TGG→TGA) dusCtRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (85/112);  total (86/112)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 198 total aligned reads displayed.

TTTGCAGCAAAGCAACCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGTAGCGCCAGCCT  >  CP009273/2223102‑2223389
                                                                                                                                                 |                                                                                                                                              
       caaagcaaCCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAg                                                                                                                                       <  2:693880/148‑1 (MQ=255)
       caaagcaaCCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAg                                                                                                                                       >  2:2915641/1‑148 (MQ=255)
        aaagcaaCCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGa                                                                                                                                      <  1:678286/148‑1 (MQ=255)
         aagcaaCCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGAt                                                                                                                                     <  1:254036/148‑1 (MQ=255)
         aagcaaCCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGAt                                                                                                                                     >  2:1313787/1‑148 (MQ=255)
         aagcaaCCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGAt                                                                                                                                     <  1:2666396/148‑1 (MQ=255)
            caaCCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTc                                                                                                                                  >  1:1869095/1‑148 (MQ=255)
             aaCCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCa                                                                                                                                 <  1:334117/148‑1 (MQ=255)
              accaccTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGACTATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCAc                                                                                                                                >  2:2672229/1‑148 (MQ=255)
              accaccTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCAc                                                                                                                                >  2:2671161/1‑148 (MQ=255)
              accaccTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCAc                                                                                                                                >  2:265985/1‑148 (MQ=255)
               ccaccTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCAcc                                                                                                                               >  2:1984302/1‑148 (MQ=255)
               ccaccTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCAcc                                                                                                                               <  1:592989/148‑1 (MQ=255)
                caccTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCg                                                                                                                              <  1:2970185/148‑1 (MQ=255)
                  ccTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGtt                                                                                                                            >  2:428706/1‑148 (MQ=255)
                       ggCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCaa                                                                                                                       >  1:2898371/1‑148 (MQ=255)
                       ggCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCaa                                                                                                                       <  2:944867/148‑1 (MQ=255)
                       ggCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCaa                                                                                                                       <  2:138444/148‑1 (MQ=255)
                        gCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAAt                                                                                                                      <  2:1965162/148‑1 (MQ=255)
                         ccACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATc                                                                                                                     >  1:1683918/1‑148 (MQ=255)
                          cACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCa                                                                                                                    <  1:1818134/148‑1 (MQ=255)
                           aCGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCAc                                                                                                                   <  1:321395/148‑1 (MQ=255)
                               cATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGa                                                                                                               <  2:3239982/148‑1 (MQ=255)
                                aTTCGCGGTTAGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGaa                                                                                                              <  1:3061172/148‑1 (MQ=255)
                                 ttCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAat                                                                                                             >  2:1363344/1‑148 (MQ=255)
                                 ttCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAat                                                                                                             >  2:1200907/1‑148 (MQ=255)
                                 ttCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAat                                                                                                             >  1:964224/1‑148 (MQ=255)
                                 ttCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAat                                                                                                             <  2:1894365/148‑1 (MQ=255)
                                 ttCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAat                                                                                                             >  2:321499/1‑148 (MQ=255)
                                 ttCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAat                                                                                                             >  1:816662/1‑148 (MQ=255)
                                 ttCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATAGCCATGCATTGTTGCGAGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAat                                                                                                             >  2:2501282/1‑148 (MQ=255)
                                  tCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAata                                                                                                            <  1:62414/148‑1 (MQ=255)
                                  tCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAata                                                                                                            >  1:2312240/1‑148 (MQ=255)
                                    gcgGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATAtt                                                                                                          <  1:2878238/148‑1 (MQ=255)
                                       gTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAg                                                                                                       >  1:2137977/1‑148 (MQ=255)
                                       gTTCGTTATAATTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAg                                                                                                       <  2:3206361/148‑1 (MQ=255)
                                        ttCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGc                                                                                                      <  2:2778794/148‑1 (MQ=255)
                                         tCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGcc                                                                                                     >  1:472881/1‑148 (MQ=255)
                                          cGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCg                                                                                                    <  1:1820385/148‑1 (MQ=255)
                                          cGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCg                                                                                                    >  1:123185/1‑148 (MQ=255)
                                           gTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGc                                                                                                   <  1:2563215/148‑1 (MQ=255)
                                           gTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGc                                                                                                   <  2:2312240/148‑1 (MQ=255)
                                           gTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGc                                                                                                   <  1:804703/148‑1 (MQ=255)
                                             tataTTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTg                                                                                                 <  1:3164015/148‑1 (MQ=255)
                                              atatTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTgg                                                                                                <  1:2718230/148‑1 (MQ=255)
                                                atTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCg                                                                                              >  1:1357951/1‑148 (MQ=255)
                                                atTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCg                                                                                              >  1:1532362/1‑148 (MQ=255)
                                                   ttACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAat                                                                                           <  1:2551322/148‑1 (MQ=255)
                                                     accaccCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAatat                                                                                         >  1:841391/1‑148 (MQ=255)
                                                      ccaccCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATc                                                                                        <  2:2125632/148‑1 (MQ=255)
                                                        accCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGc                                                                                      <  2:2066330/148‑1 (MQ=255)
                                                         cccGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGcc                                                                                     >  2:2647750/1‑148 (MQ=255)
                                                         cccGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGcc                                                                                     >  2:871733/1‑148 (MQ=255)
                                                          ccGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCa                                                                                    <  1:3349592/148‑1 (MQ=255)
                                                             gCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGACGCAGCCGCTGATCCCCATGCATAGTTGCGCGCTCTGTCAGTCCCAGATTACACCGTTGGCAATCACCGGAAAATTCAGCCGCTCGCGAATATCGCCAATc                                                                                 >  2:689508/1‑148 (MQ=255)
                                                              cTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCg                                                                                >  1:2108422/1‑148 (MQ=255)
                                                                cAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGcc                                                                              <  2:348227/148‑1 (MQ=255)
                                                                  ggTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTg                                                                            >  1:1132307/1‑148 (MQ=255)
                                                                  ggTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTg                                                                            >  1:1442022/1‑148 (MQ=255)
                                                                  ggTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTg                                                                            >  1:3472023/1‑148 (MQ=255)
                                                                   gTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGc                                                                           >  1:2456778/1‑148 (MQ=255)
                                                                     tGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCa                                                                         <  2:1365829/148‑1 (MQ=255)
                                                                         aaTATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCa                                                                     <  2:773616/148‑1 (MQ=255)
                                                                          atatTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCaa                                                                    <  1:1983899/148‑1 (MQ=255)
                                                                          atatTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCaa                                                                    <  2:3345730/148‑1 (MQ=255)
                                                                          atatTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCaa                                                                    <  2:3115717/148‑1 (MQ=255)
                                                                          atatTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCaa                                                                    <  2:913464/148‑1 (MQ=255)
                                                                          atatTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCaa                                                                    <  1:2239895/148‑1 (MQ=255)
                                                                          atatTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCaa                                                                    <  1:474818/148‑1 (MQ=255)
                                                                            atTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAata                                                                  >  2:799940/1‑148 (MQ=255)
                                                                            atTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAata                                                                  <  1:1748732/148‑1 (MQ=255)
                                                                            atTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAata                                                                  >  2:1624491/1‑148 (MQ=255)
                                                                                   gcCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCcg                                                           >  1:1750950/1‑148 (MQ=255)
                                                                                   gcCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCcg                                                           <  2:319812/148‑1 (MQ=255)
                                                                                           aCCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTaa                                                   <  1:2797316/148‑1 (MQ=255)
                                                                                           aCCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTaa                                                   <  1:1949648/148‑1 (MQ=255)
                                                                                            ccAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAAc                                                  >  2:1758761/1‑148 (MQ=255)
                                                                                            ccAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAAc                                                  >  2:1167088/1‑148 (MQ=255)
                                                                                              aaTCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAAccc                                                >  1:112456/1‑148 (MQ=255)
                                                                                                tcatcaCTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTg                                              >  1:2993387/1‑148 (MQ=255)
                                                                                                tcatcaCTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTg                                              >  1:2529999/1‑148 (MQ=255)
                                                                                                tcatcaCTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTg                                              <  2:158082/148‑1 (MQ=255)
                                                                                                tcatcaCTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTg                                              >  2:1280243/1‑148 (MQ=255)
                                                                                                   tcaCTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCCGTCAATATGCTCCGCGCGGTAACCCTGctc                                           >  2:18739/1‑148 (MQ=255)
                                                                                                   tcaCTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGctc                                           <  1:1624228/148‑1 (MQ=255)
                                                                                                     aCTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCtt                                         <  1:876576/148‑1 (MQ=255)
                                                                                                     aCTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCtt                                         <  1:2028590/148‑1 (MQ=255)
                                                                                                      cTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCttt                                        >  2:1467205/1‑148 (MQ=255)
                                                                                                        gCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCg                                      <  2:364083/148‑1 (MQ=255)
                                                                                                            cGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCg                                  >  2:2565369/1‑148 (MQ=255)
                                                                                                                gCCGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCa                              <  2:2048672/148‑1 (MQ=255)
                                                                                                                  cGCTGATCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATg                            <  2:2482067/148‑1 (MQ=255)
                                                                                                                       aTCGCCATGCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGcacca                       >  2:1422399/1‑148 (MQ=255)
                                                                                                                               gCATTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTcc               >  1:1850997/1‑148 (MQ=255)
                                                                                                                                 aTTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGt             <  1:2075809/148‑1 (MQ=255)
                                                                                                                                 aTTGTTGCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGt             <  1:2501282/148‑1 (MQ=255)
                                                                                                                                   tgttgCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGTAg           <  1:1263740/148‑1 (MQ=255)
                                                                                                                                   tgttgCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGTAg           <  1:735907/148‑1 (MQ=255)
                                                                                                                                      tgCGCGCTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGTAgcgc        >  1:352448/1‑148 (MQ=255)
                                                                                                                                          cgcTCTGTCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGTAGCGCCAGc    <  2:3344628/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
TTTGCAGCAAAGCAACCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGTAGCGCCAGCCT  >  CP009273/2223102‑2223389

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: ‑