Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099743,847,5810TC71.2% 250.4 / 189.4 272intergenic (+94/‑246)RPPX_17150/RPPX_17155tRNA‑dihydrouridine synthase C/heat‑shock protein
Reads supporting (aligned to +/- strand):  ref base T (45/32);  new base C (70/120);  total (115/157)
Fisher's exact test for biased strand distribution p-value = 1.64e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 286 total aligned reads displayed.

GATCACTCGGTTGCTGGGGGTAGAGGCTCAAGCCTTGGAGGCTTGCGTCGCCTGAACAGGCCTCTGGTACTGCACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCTGCTACCGTATTCGCGGGTAAACCCGCTCCTACGCCGGATCGAGGATGGACACAGGTCGGTGGCGGGGCACAGGCCGTCGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCATCAAACCCCGTGCAAAAAACAAATTTTTTCTAACTACCCCT  >  CP009974/3847433‑3847732
                                                                                                                                                      |                                                                                                                                                       
    aCTCGGTTGCTGGGGGTAGAGGCTCAAGCCTTGGAGGCTTGCGTCGCCTGAACAGGCCTCTGGTACTGCACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCCCCACCCGCCGTTCACCGATGCCGCCGCTAcc                                                                                                                                                   <  2:1815413/151‑1 (MQ=25)
        ggTTGCTGGGGGTAGAGGCTCAAGCCTTGGAGGCTTGCGTCGCCTGAACAGGCCTCTGGTACTGCACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCCCCACCCGCCGTTCACCGATGCCGCCGCTACCGTAt                                                                                                                                               <  1:4592424/151‑1 (MQ=255)
        ggTTGCTGGGGGTAGAGGCTCAAGCCTTGGAGGCTTGCGTCGCCTGAACAGGCCTCTGGTACTGCACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCCCCACCCGCCGTTCACCGATGCCGCCGCTACCGTAt                                                                                                                                               <  1:1296668/151‑1 (MQ=25)
         gtgggtgggggtgggggCTCACGCCTTGGAGGCTTGCGTCGCCTGAACAGGCCTCTGGTACTGCACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCCCCACCCGCCGTTCACCGATGCCGCCGCTACCGTAtt                                                                                                                                              <  2:1042528/146‑1 (MQ=21)
         gTTGCTGGGGGTAGAGGCTCAAGCCTTGGAGGCTTGCGTCGCCTGAACAGGCCTCTGGTACTGCACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCCCCACCCGCCGTTCACCGTTGCCGCCGCTGCCGTAtt                                                                                                                                              >  2:4395776/1‑151 (MQ=255)
          ttGCTGGGGGTAGAGGCTCAAGCCTTGGAGGCTTGCGTCGCCTGAACAGGCCTCTGGTACTGCACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCCCCACCCGCCGTTCACCGATGCCGCCGCTACCGTATTc                                                                                                                                             <  2:1702163/151‑1 (MQ=25)
          ttGCTGGGGGTAGAGGCTCAAGCCTTGGAGGCTTGCGTCGCCTGAACAGGCCTCTGGTACTGCACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCCCCACCCGCCGTTCACCGATGCCGCCGCTACCGTATTc                                                                                                                                             <  2:2023799/151‑1 (MQ=25)
          gtGCTGGGGGTAGAGGCTCAAGCCTTGGAGGCTTGCGTCGCCTGAACAGGCCTCTGGTACTGCACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCCCCACCCGCCGTTCACCGATGCCGCCGCTACCGTATTc                                                                                                                                             <  2:775419/150‑1 (MQ=255)
              tGGGGGTAGAGGCTCAAGCCTTGGAGGCTTGCGTCGCCTGAACAGGCCTCTGGTACTGCACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCCCCACCCGCCGTTCACCGATGCCGCCGCTACCGTATTCGCgg                                                                                                                                         <  2:5293570/151‑1 (MQ=21)
              ggggggTAGAGGCGCAAGCGTTGGAGGGTTGCGTCGCGTGGGCAGGCCGCTGGTACTGCACCTGTAGGAGCGGGTTTGCCCGCGAAGGGGCCGGCCCTGCTGGCCAATTCTCACCCGCCGTTCACTGTTGCCGCCGCTGCCGTATTCGCgg                                                                                                                                         <  2:3726375/150‑1 (MQ=17)
                  ggTAGAGGCTCAACCCTTGGAGGCTTGCGTCGCCTGAACATGCCTCTGGTACTGCACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCCCCACCCGCCGTTCACCGATGCCGCCGCTACCGTATTCGCGGGTaa                                                                                                                                     <  2:6158646/151‑1 (MQ=18)
                    tAGAGGCTCAAGCCTTGGAGGCTTGCGTCGCCTGAACAGGCCTCTGGTACTGCACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCCCCACCCGCCGTTCACCGATGCCGCCGCTACCGTATTCGCGGGTAAAc                                                                                                                                   <  2:3398459/151‑1 (MQ=255)
                    tAGAGGCTCAAGCCTTGGAGGCTTGCGTCGCCTGAACAGGCCTCTGGTACTGCACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCCCCACCCGCCGTTCACCGATGCCGCCGCTACCGTATTCGCGGGTAAAc                                                                                                                                   <  2:3287702/151‑1 (MQ=255)
                     agagGCTCAAGCCTTGGAGGCTTGCGTCGCCTGAACAGGCCTCTGGTACTGCACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCCCCACCCGCCGTTCACCGATGCCGCCGCTACCGTATTCGCGGGTAAAcc                                                                                                                                  <  2:1122864/151‑1 (MQ=21)
                       agGCTCAAGCCTTGGAGGCTTGCGTCGCCTGAACAGGCCTCTGGTACTGCACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCCCCAGCCGCCGTTCACCGATGCCGCCGCTACCGTATTCGCGGGTAAACCCg                                                                                                                                >  2:5707200/1‑151 (MQ=21)
                        ggCTCAAGCCTTGGAGGCTTGCGTCGCCTGAACAGGCCTCTGGTACTGCACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCCCCACCCGCCGTTCACCGATGCCGCCGCTACCGTATTCGCGGGTAAACCCGc                                                                                                                               <  1:3651334/151‑1 (MQ=18)
                           tCAAGCCTTGGAGGCTTGCGTCGCCTGAACAGGCCTCTGGTACTGCACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCCCCACCCGCCGTTCACCGATGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTcc                                                                                                                            <  2:568692/151‑1 (MQ=18)
                            cAAGCCTTGGAGGCTGGCGTCGCCTGAACAGGCCTCTGGTACTGCACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCCCCACCCGCCGTTCACCGATGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCt                                                                                                                           <  2:6044247/151‑1 (MQ=18)
                             aaGCCTTGGAGGCTGGCGTCGCCTGACCGGGCCTCTGGTACTGCACCTGTAGGGGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCGGGCCAATCCCCGCCCGCCGTTCACCGATGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTa                                                                                                                          <  2:27880/151‑1 (MQ=17)
                               gCCTTGGAGGCTTGCGTCGCCTGAACAGGCCTCTGGTACTGCACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCCCCACCCGCCGTTCACCGATGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTAcc                                                                                                                        >  2:3548427/1‑150 (MQ=18)
                               gCCTTGGAGGCTTGCGTCGCCTGAACAGGCCTCTGGTACTGCACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCCCCACCCGCCGTTCACCGATGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTAca                                                                                                                        <  2:2081596/151‑2 (MQ=255)
                               cggatgtaGGCTGGGGGCGCCTGAGCGGGCCTCTGGGACTGGACCTGTACGAGCGGGTTTACCCCCGAAGGGGCCGGCCCTGCTGGCCAATCCCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACg                                                                                                                        <  2:6263845/144‑1 (MQ=14)
                                   tGGAGGCTTGCGTCGCCTGAACAGGCCTCTGGTACTGCACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCCCCACCCGCCGTTCACCGATGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTAcacaga                                                                                                                    >  1:1810452/1‑146 (MQ=17)
                                    ggAGGCTTGCGTCGCCTGAACAGGCCTCTGGTACTGCACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCCCCACCCGCCGTTCACCGATGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTAcacagat                                                                                                                   >  1:1439937/1‑145 (MQ=17)
                                     gAGGCTTGCGTCGCCTGAACAGGCCTCTGGTACTGCACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCCCCACCCGCCGTTCACCGATGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTAcacagatt                                                                                                                  <  1:3328805/151‑8 (MQ=17)
                                     gAGGCTTGCGTCGCCTGAACAGGCCTCTGGTACTGCACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCCCCACCCGCCGTTCACCGATGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTAcacagatt                                                                                                                  <  1:5130705/151‑8 (MQ=17)
                                      aGGCTTGCGTCGCCTGAACAGGCCTCTGGTACTGCACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCCCCACCCGCCGTTCACCGATGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACGcaaacgg                                                                                                                 >  1:6263845/1‑145 (MQ=18)
                                       ggCTTGCGTCGCCTGAACAGGCCTCTGGTACTGCACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCCCCACCCGCCGTTCACCGATGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACCCCgcttgg                                                                                                                >  2:4312307/1‑146 (MQ=17)
                                         cTTGCGTCGCCTGAACAGGCCTCTGGTACTGCACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCCCCACCCGCCGTTCACCGATGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACCGGAtggtg                                                                                                              >  1:4282531/1‑147 (MQ=17)
                                                   cTGAACAGGCCTCTGGTACTGCACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCCCCACCCGCCGTTCACCGATGCCGCCGCTACCGTATTCGCGGGTAAAccc                                                                                                                                 >  1:6173812/1‑122 (MQ=17)
                                                                     gccgtcgtGTAGGAGCGGGAGGGCCCGGAAAGAGGCCGGCCATGCTGGCCAATTTTCACCCGCCGTTCACTGTTGCCGCCGCTGCCGTATTCGCGGGTAAACCCGCTCCTACGCCGGATCGAGGATGGACACAGGTCGGTGGCGGGGcaca                                                                                  <  1:914532/144‑1 (MQ=11)
                                                                       cgtcgtGTAGGAGCGGGTTTGCCCGCGAAGAGGCCGGCCCTGCTGGCCAATTTCCACCCGCCGTTCACTGTTGCCGCCGCTGCCGTATTCGCGGGTAAACCCGCTCCTACGCCGGATCGAGGATGGACACAGGTCGGTGGCGGGGCACAgg                                                                                <  1:3211278/146‑1 (MQ=11)
                                                                        ccgtcgTAGGAGCGGGTTCACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGTCGTTCAACGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACGCCGGATCGAGGATGGACACAGGTCGGTGGCGGGGCACAggg                                                                               <  1:3918747/146‑2 (MQ=14)
                                                                         tcgtGTAGGAGCGGGTTTGCCCGCGAAGAGGCCGGCCCTGCTGGCCAATTTTCACCCGCCGTTCACTGTTGCCGCCGCTGCCGTATTCGCGGGTAAACCCGCTCCTACGCCGGATCGAGGATGGACACAGGTCGGTGGCGGGGCACAGGcc                                                                              <  1:5707200/148‑1 (MQ=14)
                                                                         gcggggaGGAGCGCGGTTGCCCGCGAAGAGGCCGGCCCTGCTGGCCAATTTTCACCCGCCGTTCACTGTTGCCGCCGCTGCCGTATTCGCGGGTAAACCCGCTCCTACGCCGGATCGAGGATGGACACAGGTCGGTGGCGGGGCACAGGcc                                                                              <  1:4047271/145‑1 (MQ=14)
                                                                         cgtcgTAGGAGCGGGTTCACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGTCGTTCAACGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACGCCGGATCGAGGATGGACACAGGTCGGTGGCGGGGCACAgggc                                                                              <  2:5197477/147‑3 (MQ=14)
                                                                         ccgtGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATTACCACCCGCCGTTCACCGTTGCCGCTGCTACCGTATTCGCGGGTAAACCCGCTCCTACGCCTGATCGAGGATGGACACAGGTCGGTGGCGGGCGCAGGCcg                                                                             >  2:679575/4‑151 (MQ=14)
                                                                           tcgTAGGAGCGGGTTCACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGTCGTTCAACGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACGCCGGATCGAGGATGGACACAGGTCGGTGGCGGGGCACAGGGcgc                                                                            >  1:5517190/3‑150 (MQ=14)
                                                                           tcgTAGGAGCGGGTGCACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGTCGTTCAACGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACGCCGGATCGAGGATGGACACAGGTCGGTGGCGGGGCACAGGGcgt                                                                            <  2:1431668/149‑1 (MQ=17)
                                                                           gtGTAGGAGCGGGTTGACCCGCGAAGGGGCGGGCCCTGCTGGCTAATTACCACCCGCCGTTCACCGTTGCCGCTGCTACCGTATTCGCGGGTAAACCCGCTCCTACGCCTGATCGAGGATGGACACAGGTCGGTGGCGGGCGCAGGCcgtc                                                                           <  2:154902/150‑1 (MQ=11)
                                                                             gTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCACCACCCGCCGGTCACCGATGCCGCCGCTACCGTATTCGCGGGGAAACCCGCGCCTACGCATAATCGAGGATGGACACAGGTCGGGGGCTGGCGCAGGCcctcgg                                                                         >  2:2259809/1‑150 (MQ=11)
                                                                             gTAGGAGCGGGTTCACCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGTCGTTCAACGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACGCCGGATCGAGGATGGACACAGGTCGGTGGCGGGGCACAGGGcgtcg                                                                          <  1:2956811/151‑1 (MQ=17)
                                                                               aGGAGCGGGTTTGCCCGCGAAGAGGCCGGCCCTGCTGGCCAATTTTCACCCGCCGTTCACTGTTGCCGCCGCTGCCGTATTCGCGGGTAAACCCGCTCCTACGCCGGATCGAGGATGGACACAGGTCGGTGGCGGGGCACAGGCCGTCgtg                                                                        >  2:4915444/1‑151 (MQ=14)
                                                                                ggAGCGGGTTTGCCCGCGAAGAGGCCGGCCCTGCTGGCCAATTTTCACCCGCCGTTCACTGTTGCCGCCGCTGCCGTATTCGCGGGTAAACCCGCTCCTACGCCGGATCGAGGATGGACACAGGTCGGTGGCGGGGCGCCGGCCGTCgtgt                                                                       >  1:3081763/1‑151 (MQ=11)
                                                                                    cGGGTTTGCCCGCGAAGAGGCCGGCCCTGCTGGCCAATTTTCACCCGCCGTTCACTGTTGCCGCCGCTGCCGTATTCGCGGGTAAACCCGCTCCTACGCCGGATCGAGGATGGACACAGGTCGGTGGCgggc                                                                                      >  1:2650991/1‑131 (MQ=11)
                                                                                    cGGGTTTACCCGCGAAGAGGCCGGCCCTACTGGCCAATTACCACCCGCCGTTCACCGTTGCCGCTGCTACCGTATTCGCGGGTAAACCCGCTCCTACGCCTGATCGAGGATGGACACAGGTCGGTGGCGGGCGCAGGCCGTCGTGTAGGAg                                                                  >  2:3416705/1‑151 (MQ=11)
                                                                                     gggTTTGCCCGCGAAGAGGCCGGCCCTGCTGGCCAATTATCACCCGCCGTTCACCGTTGCCGCTGCTACCGTATTCGCGGGTAAACCCGCTCCTAGGCCTGATCGAGGATGGACACAGGTCGGTGGCGGGCGCAGGCcgt                                                                            >  1:3836004/1‑140 (MQ=11)
                                                                                     gggTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATTACCACCCGCCGTTCACCGTTGCCgctg                                                                                                                                                        <  1:1180600/65‑1 (MQ=16)
                                                                                     gggTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATTACCACCCGCCGTTCACCGTTGCCgctg                                                                                                                                                        >  2:1180600/1‑65 (MQ=16)
                                                                                      ggTTTGCCCGCGAAGAGGCCGGCCCTGCTGGCCAATTTTCACCCGCCGTTCACTGTTGCCGCCGCTGCCGTATTCGCGGGTAAACCCGCTACTACGCCGGATCGAGGATGGACACAGGTCGGTGGCGGGGCACAGGCCGTCGTGTAGGAGc                                                                 <  1:3130038/151‑1 (MQ=14)
                                                                                        tttgcccGCGAAGAGGCCGGCCCTGCTGGCCAATTTTCACCCGCCGTTCACTGTTGCCGCCGCTGCCGTATTCGCGGGTAAACCCGCTCCTACGCCGGATCGAGGATGGACACAGGTCGGTGGCGGGGCACAGGCCGTCGTGTAGGAGCgg                                                               <  1:6301729/147‑1 (MQ=14)
                                                                                        ttcaCCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGTCGTTCAACGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACGCCGGATCGAGGATGGACACAGGTCGGTGGCGGGGCACAGGGCGTCGTGTAGGAGCgg                                                               >  1:1642554/4‑151 (MQ=17)
                                                                                         ttgcccGCGAAGAGGCCGGCCCTGCTGGCCAATTTTCACCCGCCGTTCACTGTTGCCGCCGCTGCCGTATTCGCGGGTAAACCCGCTCCTACGCCGGATCGAGGATGGACACAGGTCGGTGGCGGGGCACAGGCCGTCGTGTAGGAGCggg                                                              <  1:4411479/148‑1 (MQ=14)
                                                                                          tgcccGCGAAGAGGCCGGCCCTGCTGGCCAATTTTCACCCGCCGTTCACTGTTGCCGCCGCTGCCGTATTCGCGGGTAAACCCGCTCCTACGCCGGATCGAGGATGGACACAGGTCGGTGGCGGGGCACAGGCCGTCGTGTAGGAGCGGGt                                                             <  1:104370/149‑1 (MQ=14)
                                                                                                  gaagaggccggccCTGCTGGCCAATTTTCACCCGCCGTTCACTGTTGCCGCCGCTGCCGTATTCGCGGGTAAACCCGCTCCTACGCCGGATCGAGGATGGACACAGGTCGGTGGCGGGGCACAGGCCGTCGTGTAGGAGCGGGTTTAcccc                                                      >  1:2427614/6‑150 (MQ=14)
                                                                                                  aaaGGGCCGGCCCTGCTGGCCAATCTCCACCCGTCGTTCAACGTTGCCGCCACCACCGTATTCGCGGGTACACCCGCTCCTACGCCGGATCGAGGAAGGACACAGGTCGGTGGCTGGGCACAGGGCGTCGTGTAGGAGCGGGTTTACCcgc                                                     >  2:2695891/1‑151 (MQ=14)
                                                                                                    ggggCGGGGCCGGGTGGCCAATCTCCGCCCGTCGTTCAACGTTGCCGCCACCACCGTATTCGCGGGTGAGCCCGCGCCGACGCCGGATCGAGGATGGACACAGGTCGGTGGCGGGGCACAGGGCGTCGTGTAGGAGCGGGTTTACCCGCGa                                                   <  2:5334430/150‑1 (MQ=14)
                                                                                                    agcggcccgccCCGCCGGGCAACTTACCCCCGCCCTTCCACCGTTGCCGCGGCAACCGTATCCGGGGGTAACCCCCCCCCACGGCCGAACCGAGGATGGACACAGTTCGGTAGCGGGCGCAGGCCGTCGTGTAGGAGCGGGTTTACCcgcg                                                    <  2:4594566/148‑1 (MQ=2)
                                                                                                    agaggccggccCTGCTGGCCAATTTTCACCCGCCGTTCACTGTTGCCGCCGCTGCCGTATTCGCGGGTAAACCCGCTCCTACGCCGGATCGCGGGTGGACACAGGTCGGTGGCGGGGCACAGGCCGCCGGGTAGGAGCGGGTTTACCcgcg                                                    >  2:4668767/4‑151 (MQ=11)
                                                                                                    agaggccggccCTGCTGGCCAATTTTCACCCGCCGTTCACTGTTGCCGCCGCTGCCGTATTCGCGGGTAAACCCGCTCCTACGCCGGATCGAGGATGGACACAGGTCGGTGGCGGGGCACAGGCCGTCGTGTAGGAGCGGGTTTACCcgcg                                                    <  1:5304405/148‑1 (MQ=17)
                                                                                                    aGGGCCGGCCCTGCTGGCCAATCTCCACCCGTCGTTCAACGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACGCCGGATCGAGGATGGACACAGGTCGGTGGCGGGGCACAGGGCGTCGTGTAGGAGCGGGTTTACCCGCGa                                                   <  1:987590/151‑1 (MQ=17)
                                                                                                     gaggccggccCTGCTGGCCAATTTTCACCCGCCGTTCACTGTTGCCGCCGCTGCCGTATTCGCGGGTAAACCCGCTCCTACGCCGGATCGAGGATGGACACAGGTCGGTGGCGGGGCACAGGCCGTCGTGTAGGAGCGGGTTTACCCGCGa                                                   >  1:1972003/3‑151 (MQ=14)
                                                                                                     gggCCGGACCTGTTGGCCTATCTCCGCCCGTCGTTCAACGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACGCCGGATCGAGAATGGACACAGGTCGGTGGCGGGGCACAGGGCGTCGTGTAGGAGCGGGTTTACCCGCGaa                                                  <  2:2925928/151‑1 (MQ=17)
                                                                                                       ggccggccCTGCTGGCCAATCTCCACCCGTCGTTCAACGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACGCCGGATCGAGGATGGACACAGGTCGGTGGCgg                                                                                        <  1:2073234/112‑1 (MQ=14)
                                                                                                        gccggcgCTGCTGGCCAAGTTTCACCCGGCGTTCACTGTGGCCGCCGCTGCTGTATTGGCGGGGAGACCCCCGCGAGCGCCGGACCGAGGCTGCACACAGGTCGGTGGCGGGGCACAGGCCGTCGTGTAGGAGCGGGTTTACCCGCGAaga                                                <  2:6229763/151‑1 (MQ=9)
                                                                                                        gccggccCTGCTGGCCAATTACCACCCGCCGTTCACCGTTGCCGCTGCCACCGTATTCGCGGGTAAACCCACTCCTACGCCGGATTGCTGCTCGACCCAGTTCGGGGGCTGGGCGCCGGCCGTGGTGTGGGACCGGGTGTACCCACGAagc                                                >  2:2459169/1‑150 (MQ=2)
                                                                                                         ctggccCTGCTGGCCAGCCGCCACCCGGCGTTCATCGTTCGCGCCGCCACGGTAATCGCGGGTAAACCCGCTCCGACGCCGGATCGAGGATGGACACAGGTCGGTGGCGGGGCACCGGGCGTCGTGTAGGGGCGGGTTTACCCGCGAagag                                               <  2:6451617/149‑1 (MQ=17)
                                                                                                           ggccCTGGTGGCCAATCTCCACCCGTCGTTCAACGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACGCCGGATCGAGGATGGACACAGGTCGGTGGCGGGGCACAGGGCGTCGTGTAGGAGCGGGTTTACCCGCGAAGAGGc                                             <  2:5887234/151‑1 (MQ=17)
                                                                                                           ggccCTGCTGGCCAATTACCACCCGCCGTTCACCGTTGCCGCTGCTACCGTATTCGCGGGTAAACCCGCTCCTACGCCTGATCGAGGATGGACACAGGTCGGTGGCGGGCGCAGGCCGTCGTGTAGGAGCGGGTTTACCCGCGAAGAGGcc                                            >  2:640038/1‑151 (MQ=14)
                                                                                                            gccCTGCTGGCCAATCTCCACCCGTCGTTCAACGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACGCCGGATCGAGGATGGACACAGGTCGGTGGCGGGGCACAGGGCGTCGTGTAGGAGCGGGTTTACCCGCGAAGAGGcc                                            <  2:2203279/151‑1 (MQ=17)
                                                                                                             cccTGCTGGCCAATTTTCACCCGCCGTTCACTGTTGCCGCCGCTGCCGTATTCGCGGGTAAACCCGCTCCTACGCCGGATCGAGGATGGACACAGGTCGGTGGCGGGGCACAGGCCGTCGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCa                                           >  2:1203316/1‑151 (MQ=17)
                                                                                                              ccTGCGGGCCAATTTTCACCCGCCGTTCACTGTTGCCGCCGCTGCCGGATTCGCGGGTAAACCCGCTCCTACGCCGGATCGAGGATGGACACAGGTCGGTGGCGGGGCACAGGCCGTCGTGTAGGAGCGGGTTTACCCGCGAAGAGGccgg                                          <  1:3732105/151‑3 (MQ=14)
                                                                                                                tgctgGCCAATTACCACCCGCCGTTCACCGTTGCCGCTGCTACCGTATTCGCGGGTAAACCCGCTCCTACGCCTGATCGAGGATGGACACAGGTCGGTGGCGGGCGCAGGCCGTCGTGTAGGAGCGGGGTTACCCGCGAAGAGGCCATCaa                                       >  2:3163406/1‑151 (MQ=17)
                                                                                                                tgctgGCCAATCTCCACCCGTCGTTCAACGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACGCCGGATCGAGGATGGACACAGGTCGGTGGCGGGGCACAGGGCGTCGTGTAGGAGCGGGTTTACCCGCGAAGAGGccggcc                                        >  1:6345228/1‑147 (MQ=17)
                                                                                                                 gctgGCCAATTTTCACCCGCCGTTCACCGTTGCCGCTGCTACCGTATTCGCGGGTAAACCCGCTCCTACGCCTGATCGAGGATGGACACAGGTCGGTGGCGGGCGCAGGCCGTCGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCATCaaa                                      >  1:2570781/1‑151 (MQ=17)
                                                                                                                   tgGCCAATTTTCACCCGCCGTTCACTGTTGCCGCCGCTGCCGTATTCGCGGGTAAACCCGCTCCTACGCCGGATCGAGGATGGACACAGGTCGGTGGCGGGGCACAGGCCGTCGTGTAGGAGCGGGTTTACCCGCGAAGAGGccggcccta                                     >  1:2487396/1‑144 (MQ=14)
                                                                                                                   tgGCCAATTACCACCCGCCGTTCACCGTTGCCGCTGCTACCGTATTCGCGGGTAACCCCGCTCCTACGCCTGACCGAGCATGGACACAGGTCGCTGGCGGGCGGAGGCCGTCGTGTAGGCGCGGGTTTACCCGCGAGGAGGCGATCACAcc                                    >  2:4451853/1‑151 (MQ=14)
                                                                                                                     gCCAATTCCCACCCGCCGTTCACCGTTGCCGCTGCTACCGTATTCGCGGGTAAACCCGCTCCTACGCCTGATCGAGGATGGACACAGGTCGGTGGCGGGCGCAGGCCGTCGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCATCAAAcccc                                  <  1:5830004/151‑1 (MQ=255)
                                                                                                                     gCCAATTACCACCCGCCGTTCACCGTTGCCGCTGCTACCGTATTCGCGGGTAAACCCGCTCCTACGCCTGATCGAGGATGGACACAGGTCGGTGGCGGGCGCAGGCCGTCGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCATCAAAcccc                                  <  1:5053370/151‑1 (MQ=255)
                                                                                                                     ctcacgttCCCCCCGCCGTTCACTGTTGCCGCCGCTGCCGTATCCGCGGGTAAACCCGCTCCTACGCCTGATCGAGGAGGGACACAGGTCGGTGGCGGGCGCAGGCCGTCGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCATCAAAcccc                                  <  1:2687633/145‑1 (MQ=17)
                                                                                                                      ccAATTACCACCCGCCGTTCACCGTTGCCGCTGCTACCGTATTCGCGGGTAAACCCGCTCCTACGCCTGATCGAGGATGGACACAGGTCGGTGGCGGGCGCAGGCCGTCGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCATCAAACCCCg                                 >  1:1822678/1‑151 (MQ=17)
                                                                                                                      ccAATTACCACCCGCCGTTCACCGTTGCCGCTGCTACCGTATTCGCGGGTAAACCCGCTCCTACGCCTGATCGAGGATGGACACAGGTCGGTGGCGGGCGCAGGCCGTCGTGTAGGAGCAGGTTTACCCGCGAAGAGGCCATCGAACCCCg                                 >  2:3938293/1‑151 (MQ=17)
                                                                                                                       cAATTACCACCCGCCGTTCACCGTTGCCGCTGCTACCGTATTCGCGGGTAAACCCGCTCCTACGCCTGATCGAGGATGGACACAGGTCGGTGGCGGGCGCAGGCCGTCGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCATCAAACCCCGt                                <  1:640517/151‑1 (MQ=17)
                                                                                                                       cAATTACCACCCGCCGTTCACCGTTGCCGCTGCTACCGTATTCGCGGGTAAACCCGCTCCTACGCCTGATCGAGGATGGACACAGGTCGGTGGCGGGCGCAGGCCGTCGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCATCAAACCCCGt                                <  1:3366670/151‑1 (MQ=17)
                                                                                                                         atctCCACCCGTCGTTCAACGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACGCCGGATCGAGGATGGACACAGGTCGGTGGCGGGGcac                                                                                   >  1:1428145/4‑99 (MQ=17)
                                                                                                                           taccACCCGCCGTTCACCGTTGCCGCTGCTACCGTATTCGCGGGTAAACCCGCTCCTACGCCTGATCGAGGATGGACACAGGTCGGTGGCGGGCGCAGGCCGTCGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCATCACACCCCGTGCaa                            >  2:3826082/3‑151 (MQ=17)
                                                                                                                             ccccacgccgGTCACCTTTGCCGCCGCTACCGCAAGCGCGGGTAATCCCGCTCCTACGCCGGATCGAGGATGGACACAGTTCGGTGGCGGGCGCAGGCCGTCGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCATCAAACCCCGTGCaaaa                          <  1:1034141/146‑1 (MQ=18)
                                                                                                                              cACCCGCCGTTCACCGTTGCCGCTGCTACCGTATTCGCGGGTAAACCCGCTCCTACGCCTGATCGAGGATGGACACAGGTCGGTGGCGGGCGCAGGCCGTCGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCATCAAACCCCGTGCaaaaa                         >  1:2771815/1‑151 (MQ=255)
                                                                                                                              cACCCGCCGTTCACCGTTGCCGCTGCTACCGTATTCGCGGGTAAACCCGCTCCTACGCCTGATCGAGGATGGACACAGGTCGGTGGCGGGCGCAGGCCGTCGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCATCAAACCCCGTGCaaaaa                         <  1:2259809/151‑1 (MQ=255)
                                                                                                                              cACCCGCCGTTCACCGTTGCCGCTGCTACCGTATTCGCGGGTAAACCCGCTCCTACGCCTGATCGAGGATGGACACAGGTCGGTGGCGGGCGCAGGCCGTCGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCATCAAACCCCGTGCaaaaa                         >  1:4699923/1‑151 (MQ=255)
                                                                                                                               ccccGCCGTTCACCGTTGCCGCTGCTACCGTATTCGCGGGTAAACCCGCTCCTACGCCTGATCGAGGATGGACACAGGTCGGTGGCGGGCGCAGGCCGTCGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCATCAAACCCCGTGCaaaaaa                        <  1:2368582/150‑1 (MQ=255)
                                                                                                                                cccGCCGTTCACCGTTGCCGCTGCTACCGTATTCGCGGGTAAACCCGCTCCTACGCCTGATCGAGGATGGACACAGGTCGGTGGCGGGCGCAGGCCGTCGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCATCAAACCCCGTGCAAAAAAc                       <  2:2570781/151‑1 (MQ=255)
                                                                                                                                      gTTCACCGTTGCCGCTGCTACCGTATTCGCGGGTAAACCCGCTCCTACGCCTGATCGAGGATGGACACAGGTCGGTGGCGGGCGCAGGCCGTCGTGTAGGAGCGGGTTTACCCGCGAAGAGTCCATCAAACCCCGTGCAAAAAACAAAttt                 >  1:4514227/1‑151 (MQ=255)
                                                                                                                                        tCACCGTTGCCGCTGCTACCGTATTCGCGGGTAAACCCGCTCCTACGCCTGATCGAGGACGGACACAGGTCGGTGGCGGGCGCAGGCCGTCGTGTAGGCTCGTGTTTACACGCGAAGAGGACATCAAACTCCGTGCAAAAAACAAAtattt               >  2:4517041/1‑151 (MQ=18)
                                                                                                                                           cccGTTGCCCCTGCTACCTTATCCGCGGGTAACCCCGCTCCACCGCCGGATGGAGGTTGAACACAGGCGGGTGGCGGGCGCAGGCCGTCGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCATCAAACCCCGTGCAAAAAACAAATTTTTTc             <  1:4105115/150‑1 (MQ=18)
                                                                                                                                              gTTGCCGCTGCTACCGTATTCGCGGGTAAACCCGCTCCTACGCCTGATCGAGGATGGACACAGGTCGGTGGCGGGCGCAGGCCGTCGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCATCAAACCCCGTGCAAAAAACAAATTTTTTCTaa          <  1:3571845/151‑1 (MQ=255)
                                                                                                                                               ttGCCGCTGCTACCGTATTCGCGGGTAAACCCGCTCCTACGCCTGATCGAGGATGGACACAGGTCGGTGGCGGGCGCAGGCCGTCGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCATCAAACCCCGTGCAAAAAACAAATTTTTTCTAAc         <  1:5745608/151‑1 (MQ=255)
                                                                                                                                                agCCGCTGCGAACTAATTCCCGGCTAAACACGATCCTACTCCTGATCGAGGAGGGACACAGGTAGGTGGCGGGCGCAGGCCGTCGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCATCAAACCCCGTGCAAAAAACAAATTTTTTCTAACt        <  2:6543501/150‑1 (MQ=21)
                                                                                                                                                   cGCTGCTACCGTATTCGCGGGTAAACCCGCTCCTACGCCTGATCGAGGATGGACACAGGTCGGTGGCGGGCGCAGGCCGTCGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCATCAAACCCCGTGCAAAAAACAAATTTTTTCTAACTAcc     >  2:404991/1‑151 (MQ=255)
                                                                                                                                                     ctgctACCGTATTCGCGGGTAAACCCGCTCCTACGCCTGATCGAGGATGGACACAGGTCGGTGGCGGGCGCAGGCCGTCGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCATCAAACCCCGTGCAAAAAACAAATTTTTTCTAACTAcccc   >  1:4793589/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                       
GATCACTCGGTTGCTGGGGGTAGAGGCTCAAGCCTTGGAGGCTTGCGTCGCCTGAACAGGCCTCTGGTACTGCACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCTGCTACCGTATTCGCGGGTAAACCCGCTCCTACGCCGGATCGAGGATGGACACAGGTCGGTGGCGGGGCACAGGCCGTCGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCATCAAACCCCGTGCAAAAAACAAATTTTTTCTAACTACCCCT  >  CP009974/3847433‑3847732

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: