Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,404,059 C→A 100% R141S (CGC→AGC)  dusB → tRNA‑dihydrouridine synthase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,404,0590CA98.8% 1103.4 / 2.2 335R141S (CGC→AGC) dusBtRNA‑dihydrouridine synthase B
Reads supporting (aligned to +/- strand):  ref base C (3/1);  new base A (167/163);  total (171/164)
Fisher's exact test for biased strand distribution p-value = 6.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01
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 335 total aligned reads displayed.

CAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAGTT  >  CP009273/3403912‑3404205
                                                                                                                                                   |                                                                                                                                                  
cAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGc                                                                                                                                                  <  1:1152122/150‑1 (MQ=255)
  gATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCAc                                                                                                                                                <  1:1820256/150‑1 (MQ=255)
       ttGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCt                                                                                                                                           <  1:3785667/150‑1 (MQ=255)
         gATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTgg                                                                                                                                         <  1:5040897/150‑1 (MQ=255)
          atatCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTggg                                                                                                                                        >  1:2844204/1‑150 (MQ=255)
          atatCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTggg                                                                                                                                        <  2:1210067/150‑1 (MQ=255)
            atCAAGATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCa                                                                                                                                      <  2:2300444/150‑1 (MQ=255)
              cAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCAcc                                                                                                                                    <  2:443894/150‑1 (MQ=255)
              cAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCAcc                                                                                                                                    >  2:1436245/1‑150 (MQ=255)
                atatGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCgg                                                                                                                                  >  1:1058567/1‑150 (MQ=255)
                atatGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCgg                                                                                                                                  <  1:634060/150‑1 (MQ=255)
                 tatGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGa                                                                                                                                 >  1:4423828/1‑150 (MQ=255)
                   tGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAac                                                                                                                               <  2:4246797/150‑1 (MQ=255)
                   tGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAac                                                                                                                               >  2:783428/1‑150 (MQ=255)
                     ggTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAacac                                                                                                                             <  1:1927087/150‑1 (MQ=255)
                       ttGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACAcca                                                                                                                           <  2:3503839/150‑2 (MQ=255)
                       ttGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCg                                                                                                                           >  2:5013583/1‑150 (MQ=255)
                       ttGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCATCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCg                                                                                                                           <  1:3558854/150‑1 (MQ=255)
                        tGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGt                                                                                                                          >  1:1368592/1‑150 (MQ=255)
                         gCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTa                                                                                                                         >  1:437815/1‑150 (MQ=255)
                             ggCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTg                                                                                                                     >  2:4670527/1‑150 (MQ=255)
                                tAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGa                                                                                                                  <  2:3350850/150‑1 (MQ=255)
                                tAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGa                                                                                                                  <  2:2438481/150‑1 (MQ=255)
                                 aaaaaaGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGaa                                                                                                                 <  1:4265513/150‑1 (MQ=255)
                                 aaaaaaGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGaa                                                                                                                 >  1:3329215/1‑150 (MQ=255)
                                 aaaaaaGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGaa                                                                                                                 <  1:4691792/150‑1 (MQ=255)
                                 aaaaaaGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGaa                                                                                                                 >  1:2233050/1‑150 (MQ=255)
                                   aaaaGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAaga                                                                                                               >  2:2570015/1‑150 (MQ=255)
                                     aaGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAgaga                                                                                                             >  1:4891297/1‑150 (MQ=255)
                                     aaGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAgaga                                                                                                             >  2:130504/1‑150 (MQ=255)
                                      aGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGAt                                                                                                            >  1:3282987/1‑150 (MQ=255)
                                       gTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGAtt                                                                                                           >  1:944446/1‑150 (MQ=255)
                                       gTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGAtt                                                                                                           <  1:2006075/150‑1 (MQ=255)
                                        tGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTg                                                                                                          <  2:3054784/150‑1 (MQ=255)
                                         gAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGc                                                                                                         >  1:2334425/1‑150 (MQ=255)
                                          aaTCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGcc                                                                                                        <  1:1037501/150‑1 (MQ=255)
                                          aaTCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGcc                                                                                                        >  2:3576517/1‑150 (MQ=255)
                                          aaTCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGcc                                                                                                        >  1:4842953/1‑150 (MQ=255)
                                           aTCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGccc                                                                                                       >  2:3306792/1‑150 (MQ=255)
                                             cGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCaa                                                                                                     <  2:4134656/150‑1 (MQ=255)
                                               cAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACt                                                                                                   <  1:4027210/150‑1 (MQ=255)
                                                aaGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTg                                                                                                  >  1:4818913/1‑150 (MQ=255)
                                                aaGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTg                                                                                                  >  1:512674/1‑150 (MQ=255)
                                                    tCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTg                                                                                              <  1:2911967/150‑1 (MQ=255)
                                                     cGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGa                                                                                             <  1:158814/150‑1 (MQ=255)
                                                      gCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGaa                                                                                            >  1:1562020/1‑150 (MQ=255)
                                                        aGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGa                                                                                          >  1:3152545/1‑150 (MQ=255)
                                                        aGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGa                                                                                          >  2:3124967/1‑150 (MQ=255)
                                                            tCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACtgt                                                                                      >  1:4649169/1‑150 (MQ=255)
                                                             cAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTACCTGCGAAGAGATTGCCCAACTGGCTGAAGACtgtg                                                                                     >  2:4380206/1‑150 (MQ=255)
                                                                cccTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCa                                                                                  <  1:1444800/150‑1 (MQ=255)
                                                                 ccTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCAt                                                                                 >  2:1024633/1‑150 (MQ=255)
                                                                  ctctTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCAtt                                                                                <  2:3309087/150‑1 (MQ=255)
                                                                    ctTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCa                                                                              <  1:1999030/150‑1 (MQ=255)
                                                                      tGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAgg                                                                            >  2:2014236/1‑150 (MQ=255)
                                                                           tACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTg                                                                       <  1:4547089/150‑1 (MQ=255)
                                                                            aCCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGa                                                                      >  2:2331663/1‑150 (MQ=255)
                                                                             cccGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGAc                                                                     >  2:4989344/1‑150 (MQ=255)
                                                                               cGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCa                                                                   >  1:1714127/1‑150 (MQ=255)
                                                                               cGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCa                                                                   >  2:1721869/1‑150 (MQ=255)
                                                                                ggATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCAt                                                                  >  1:3846957/1‑150 (MQ=255)
                                                                                    gTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCAt                                                              >  2:1584396/1‑150 (MQ=255)
                                                                                     tCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATg                                                             <  2:40101/150‑1 (MQ=255)
                                                                                      cGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATgg                                                            >  1:2664524/1‑150 (MQ=255)
                                                                                      cGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATgg                                                            <  1:745492/150‑1 (MQ=255)
                                                                                         tAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCg                                                         <  2:4340485/150‑1 (MQ=255)
                                                                                         tAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCg                                                         >  2:2895898/1‑150 (MQ=255)
                                                                                            aTCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTac                                                      <  2:3777476/150‑1 (MQ=255)
                                                                                                 tCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACAcgcg                                                 >  2:348015/1‑150 (MQ=255)
                                                                                                   cTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGcc                                               >  2:861631/1‑150 (MQ=255)
                                                                                                    ttACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGcca                                              <  2:3449095/150‑2 (MQ=255)
                                                                                                    ttACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCt                                              <  1:3981287/150‑1 (MQ=255)
                                                                                                      aCCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCtgt                                            <  2:4496944/150‑1 (MQ=255)
                                                                                                        cGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCtgttt                                          <  1:4198389/150‑1 (MQ=255)
                                                                                                         gAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCtgtttg                                         <  1:1272557/150‑1 (MQ=255)
                                                                                                          aGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCtgtttgt                                        >  2:4739825/1‑150 (MQ=255)
                                                                                                             tcgtcAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCa                                     <  2:4744753/150‑1 (MQ=255)
                                                                                                              cgtcAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCaa                                    >  1:3409724/1‑150 (MQ=255)
                                                                                                               gtcAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAAt                                   >  1:3074929/1‑150 (MQ=255)
                                                                                                                 cAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATgg                                 >  1:3719845/1‑150 (MQ=255)
                                                                                                                       aGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGc                           >  2:1422247/1‑150 (MQ=255)
                                                                                                                         tGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTg                         <  1:791525/150‑1 (MQ=255)
                                                                                                                             cGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTa                     <  1:2941204/150‑1 (MQ=255)
                                                                                                                              gTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTaa                    <  1:3088728/150‑2 (MQ=255)
                                                                                                                              gTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACAATTCATGGCCGTACACGCGCCTGTTTGATCAATGGAGAAGCTGAGTAc                    >  2:1890074/1‑150 (MQ=255)
                                                                                                                               ttCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACg                   >  2:4537140/1‑150 (MQ=255)
                                                                                                                                tCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGa                  >  1:1319471/1‑150 (MQ=255)
                                                                                                                                  cTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACa                >  2:3206560/1‑150 (MQ=255)
                                                                                                                                   tGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAg               <  1:2016040/150‑1 (MQ=255)
                                                                                                                                   tGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAg               <  2:4863279/150‑1 (MQ=255)
                                                                                                                                        cccTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTc          <  2:2806975/150‑1 (MQ=255)
                                                                                                                                        cccTGAAGATTAGCACCGGCTGGGCACAGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTc          >  1:4473257/1‑150 (MQ=255)
                                                                                                                                         ccTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCg         <  2:2969096/150‑1 (MQ=255)
                                                                                                                                          cTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCgg        >  2:1233903/1‑150 (MQ=255)
                                                                                                                                             aaGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCa     <  2:4370185/150‑1 (MQ=255)
                                                                                                                                             aaGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCa     <  2:3933012/150‑1 (MQ=255)
                                                                                                                                             aaGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCa     <  2:3839290/150‑1 (MQ=255)
                                                                                                                                              aGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAg    <  1:2504644/150‑1 (MQ=255)
                                                                                                                                               gATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAGt   >  2:2190641/1‑150 (MQ=255)
                                                                                                                                                aTTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAGtt  <  1:2144326/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAGTT  >  CP009273/3403912‑3404205

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: