Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,252,556 A→T 5.6% L320F (TTA→TTT dusA → tRNA‑dihydrouridine synthase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,252,5560AT 5.6% 836.3 / 28.4 265L320F (TTA→TTTdusAtRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base A (120/130);  new base T (8/7);  total (128/137)
Fisher's exact test for biased strand distribution p-value = 7.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 271 total aligned reads displayed.

CCGATGCCGATCCGGTGGCGGTAGTGCGCGCCATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTAAATCAA  >  CP009273/4252408‑4252705
                                                                                                                                                    |                                                                                                                                                     
ccGATGCCGATCCGGTGGCGGTAGTGCGCGCCATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTaa                                                                                                                                                      <  1:4430778/150‑1 (MQ=255)
   aTGCCGATCCGGTGGCGGTAGTGCGCGCCATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTg                                                                                                                                                   <  2:1836114/150‑1 (MQ=255)
     gCCGATCCGGTGGCGGTAGTGCGCGCCATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGaa                                                                                                                                                 >  2:3181627/1‑150 (MQ=255)
       cGATCCGGTGGCGGTAGTGCGCGCCATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTTAGTGaaaa                                                                                                                                               >  2:1615630/1‑150 (MQ=255)
             ggTGGCGGTAGTGCGCGCCATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCa                                                                                                                                         <  2:1413360/150‑1 (MQ=255)
             ggTGGCGGTAGTGCGCGCCATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGGTATTTAAGTGAAAATGCCCa                                                                                                                                         <  2:3494998/150‑1 (MQ=255)
              gTGGCGGTAGTGCGCGCCATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGTCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCAt                                                                                                                                        <  1:2842677/150‑1 (MQ=255)
              gTGGCGGTAGTGCGCGCCATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCAt                                                                                                                                        <  2:3969684/150‑1 (MQ=255)
              gTGGCGGTAGTGCGCGCCATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCAt                                                                                                                                        >  2:1702595/1‑150 (MQ=255)
                ggCGGTAGTGCGCGCCATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATaa                                                                                                                                      <  2:2394748/150‑1 (MQ=255)
                   ggTAGTGCGCGCCATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGc                                                                                                                                   <  1:3919253/150‑1 (MQ=255)
                     tAGTGCGCGCCATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCgg                                                                                                                                 >  1:2872498/1‑150 (MQ=255)
                           gcgcCATGTATCCGTACATTGAGCGTGAACTCAGCCAGTGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAg                                                                                                                           <  1:4450554/150‑1 (MQ=255)
                             gcCATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGAc                                                                                                                         <  2:4949526/150‑1 (MQ=255)
                              ccATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACa                                                                                                                        >  2:4004183/1‑150 (MQ=255)
                                 tGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTa                                                                                                                     >  2:2463307/1‑150 (MQ=255)
                                     tCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAAtgt                                                                                                                 >  1:3318007/1‑150 (MQ=255)
                                      ccGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAAtgtg                                                                                                                >  1:1357691/1‑150 (MQ=255)
                                       cGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGc                                                                                                               >  2:1299290/1‑150 (MQ=255)
                                        gTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCt                                                                                                              <  1:494087/150‑1 (MQ=255)
                                        gTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCt                                                                                                              <  2:1758466/150‑1 (MQ=255)
                                        gTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCt                                                                                                              <  1:3081878/150‑1 (MQ=255)
                                          aCATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTgg                                                                                                            <  1:1168916/150‑1 (MQ=255)
                                          aCATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTgg                                                                                                            <  1:1052510/150‑1 (MQ=255)
                                            aTTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCAGATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACAGTAATGTGCTGGaa                                                                                                          >  2:2733197/1‑150 (MQ=255)
                                               gAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAacac                                                                                                       >  2:2526927/1‑150 (MQ=255)
                                                aGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACAcg                                                                                                      <  2:624825/150‑1 (MQ=255)
                                                   gTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACgcgc                                                                                                   >  1:583767/1‑150 (MQ=255)
                                                   gTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACgcgc                                                                                                   >  1:583814/1‑150 (MQ=255)
                                                     gAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCAATAAAGAGGGTGCAGACATTAATGTGCTGGAACACGCGCTc                                                                                                 >  2:2189942/1‑150 (MQ=255)
                                                      aaCTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCa                                                                                                <  1:1271727/150‑1 (MQ=255)
                                                        cTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCaaa                                                                                              >  1:4016120/1‑150 (MQ=255)
                                                        cTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGGTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCaaa                                                                                              <  1:622241/150‑1 (MQ=255)
                                                         tCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAAc                                                                                             >  2:4909021/1‑150 (MQ=255)
                                                              cAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACtggtg                                                                                        <  2:1759463/150‑1 (MQ=255)
                                                              cAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACtggtg                                                                                        <  1:4175580/150‑1 (MQ=255)
                                                                ggggACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGc                                                                                      >  1:4747862/1‑150 (MQ=255)
                                                                 gggACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCg                                                                                     <  1:2300274/150‑1 (MQ=255)
                                                                    aCGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGAt                                                                                  <  1:3854884/150‑1 (MQ=255)
                                                                    aCGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGAt                                                                                  <  2:1889558/150‑1 (MQ=255)
                                                                        aTCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGc                                                                              >  2:1789690/1‑150 (MQ=255)
                                                                           tcGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGtt                                                                           <  1:101747/150‑1 (MQ=255)
                                                                            cGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAACATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTa                                                                          <  2:1980851/150‑1 (MQ=255)
                                                                             ggCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTaa                                                                         <  2:2358728/150‑1 (MQ=255)
                                                                             ggCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTaa                                                                         >  1:318814/1‑150 (MQ=255)
                                                                             ggCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTaa                                                                         <  2:4057575/150‑1 (MQ=255)
                                                                              gCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAAc                                                                        <  1:3785564/150‑1 (MQ=255)
                                                                                 atatTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACtta                                                                     <  2:3851460/150‑2 (MQ=255)
                                                                                  tatTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACtttt                                                                    >  1:3146253/1‑150 (MQ=255)
                                                                                  tatTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACtttt                                                                    >  1:3145872/1‑150 (MQ=255)
                                                                                       cccGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCa                                                               <  1:633681/150‑1 (MQ=255)
                                                                                        ccGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGATGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCaa                                                              >  1:1951418/1‑150 (MQ=255)
                                                                                         cGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCaaa                                                             >  1:2380889/1‑150 (MQ=255)
                                                                                         cGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCaaa                                                             >  1:2043355/1‑150 (MQ=255)
                                                                                          ggCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCaaaa                                                            >  2:3226280/1‑150 (MQ=255)
                                                                                           gCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCaaaaa                                                           <  1:1856172/150‑1 (MQ=255)
                                                                                            cATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAag                                                          <  1:5092626/150‑1 (MQ=255)
                                                                                            cATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAag                                                          >  1:1786202/1‑150 (MQ=255)
                                                                                              tatGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTGAACTTTTCACCAAAAagta                                                        >  2:2127710/1‑150 (MQ=255)
                                                                                                 gTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTc                                                     >  1:1419797/1‑150 (MQ=255)
                                                                                                  ttGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTTAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCa                                                    <  2:188307/150‑1 (MQ=255)
                                                                                                  ttGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCa                                                    >  1:2711873/1‑150 (MQ=255)
                                                                                                       ctagTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTTCAGAAATTAATGTGATGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTGTTCACCAAAAAGTAGTCAAATTc                                               >  1:1674579/4‑150 (MQ=255)
                                                                                                         tGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTcac                                             >  2:4992214/1‑150 (MQ=255)
                                                                                                         tGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTcac                                             >  2:2809869/1‑150 (MQ=255)
                                                                                                         tGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGAAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTGAACTTTTCACCAAAAAGTAGTCAAATTcac                                             >  2:662047/1‑150 (MQ=255)
                                                                                                          gTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTcacc                                            >  2:607804/1‑150 (MQ=255)
                                                                                                           ttCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAAGTcacca                                           >  1:2659590/1‑150 (MQ=255)
                                                                                                            tCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTcaccac                                          <  1:3656279/150‑1 (MQ=255)
                                                                                                             ccAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACg                                         >  1:2977473/1‑150 (MQ=255)
                                                                                                              cAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGc                                        >  1:3695116/1‑150 (MQ=255)
                                                                                                              cAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGc                                        >  1:2432939/1‑150 (MQ=255)
                                                                                                               aGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTTAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGcc                                       >  1:1273441/1‑150 (MQ=255)
                                                                                                                 ggTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTTAGTGAAAATGCCCATAAATCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCt                                     >  2:125638/1‑150 (MQ=255)
                                                                                                                 ggTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACAAAAAAGTAGTCAAATTCACCACGCCCt                                     >  1:2641530/1‑150 (MQ=255)
                                                                                                                  gTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTg                                    <  1:416651/150‑1 (MQ=255)
                                                                                                                  gTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTg                                    <  2:769104/150‑1 (MQ=255)
                                                                                                                       ccTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCAc                               >  1:253613/1‑150 (MQ=255)
                                                                                                                       ccTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCAc                               >  1:749186/1‑150 (MQ=255)
                                                                                                                        cTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCAcc                              <  1:4909021/150‑1 (MQ=255)
                                                                                                                         tGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCg                             >  1:4353228/1‑150 (MQ=255)
                                                                                                                           gcgcgcGGCAGTGGCGGCGTTATTTTAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTc                           <  1:1452914/150‑1 (MQ=255)
                                                                                                                              cgcgGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTcgcg                        >  1:161835/1‑150 (MQ=255)
                                                                                                                              cgcgGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTcgcg                        <  2:4278647/150‑1 (MQ=255)
                                                                                                                               gcgGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCgg                       <  1:3544441/150‑1 (MQ=255)
                                                                                                                                cgGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTGAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCggg                      >  1:4488146/1‑150 (MQ=255)
                                                                                                                                  gCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGc                    <  2:2666059/150‑1 (MQ=255)
                                                                                                                                  gCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGc                    <  1:4286588/150‑1 (MQ=255)
                                                                                                                                     gTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGtt                 <  2:3731069/150‑1 (MQ=255)
                                                                                                                                     gTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGtt                 <  1:1000226/150‑1 (MQ=255)
                                                                                                                                       ggcggcGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGtttt               <  1:4626995/150‑1 (MQ=255)
                                                                                                                                          ggcgTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGTTTTGCt            >  2:4420259/1‑150 (MQ=255)
                                                                                                                                           gcgTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGTTTTGCTg           <  2:4553929/150‑1 (MQ=255)
                                                                                                                                            cgTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGt          >  1:4245735/1‑150 (MQ=255)
                                                                                                                                             gTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGtt         >  2:4387845/1‑150 (MQ=255)
                                                                                                                                             gTTATTTAAGTGAAAATGCCCATAAAGCGGGGGCAGACCTTAATGGGCTGGAACACGCGCTCCACCTGGTGGCGGATAAGCGTTAACTTTTCCCCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGtt         <  1:917303/150‑1 (MQ=255)
                                                                                                                                               tATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTaa       <  1:1530387/150‑1 (MQ=255)
                                                                                                                                                aTTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTaaa      >  2:789149/1‑150 (MQ=255)
                                                                                                                                                 tttAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTAAAt     <  2:4062062/150‑1 (MQ=255)
                                                                                                                                                    aaGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTAAATCaa  <  2:3099960/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCGATGCCGATCCGGTGGCGGTAGTGCGCGCCATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTAAATCAA  >  CP009273/4252408‑4252705

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: