Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NZ_CP022686 1,633,846 T→G 24.8% N253T (AAC→ACC)  yqeB ← selenium‑dependent molybdenum cofactor biosynthesis protein YqeB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0226861,633,8460TG24.8% 312.8 / 59.7 242N253T (AAC→ACC) yqeBselenium‑dependent molybdenum cofactor biosynthesis protein YqeB
Reads supporting (aligned to +/- strand):  ref base T (49/129);  new base G (59/5);  total (108/134)
Fisher's exact test for biased strand distribution p-value = 2.24e-20
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 245 total aligned reads displayed.

CACTTCCAACATGATCACTTTAAAACCCGCATGATACAGACGTAGCGCTACACCACTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAGCACGCTAATGGCGATCTCCTGCGGCGTTTCCGCACCTATGTTGTAACCAACGGGCGCATGTAAACGGGCAATATGTTCTTCAGCCACGCCTTTCTCACGCAATTGGCGCAG  >  NZ_CP022686/1633707‑1633985
                                                                                                                                           |                                                                                                                                            
    tCCAACATGATCACTTTAAAACCCGCATGATACAGACGTAGCGCTACACCACTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGCAGGCTTCATCATCAACCCACCCGGCGCGTTAtt                                                                                                                                          >  2:953864/1‑140 (MQ=255)
       aaCATTATCACTTTAAAAACCACATGATACAGACGTAGCGCTACACCACTGGCGATATCCCCCGCACCGAGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTAttttt                                                                                                                                       <  7:105360/140‑1 (MQ=255)
       aaCATGATCACTTTAAAACCCGCATGATACAGACGTAGCGCTACACCACTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACACGGCGCGTTAttt                                                                                                                                         <  4:1027049/138‑1 (MQ=255)
       aaCATGATCACATTAAAACCCGCATGATACAGACGTAGCTCTCCACCACTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTAttttt                                                                                                                                       <  7:75359/140‑1 (MQ=255)
           tGATCACTTTAAAACCCGCATGATACAGACGTAGCGCTACACCACTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCAcc                                                                                                                                   <  6:134993/140‑1 (MQ=255)
           tGATCACTTTAAAAACCGCATGATACAGTCGTAGCGCTACACCACTGGCGCTATCCCCCTCACCGCGAACCACCACCAGCTTGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCAcc                                                                                                                                   <  5:363267/140‑1 (MQ=255)
               cACTTTAAAACCCGCATGATACAGACGTAGCGCTACACCACTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGGTAttgttgaacgg                                                                                                                                 >  3:847655/1‑132 (MQ=255)
                         cccGCATGATACAGACGTAGCGCTACACCACTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGGTGTTTTTAACCTGTAATATTTCt                                                                                                                       >  3:1043623/1‑138 (MQ=255)
                            gCATGATACAGACGTAGCGCTACACCACTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTc                                                                                                                                      <  7:330608/120‑1 (MQ=255)
                            gCATGATACAGACGTAGCGCTACACCACTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTc                                                                                                                                      >  8:330608/1‑120 (MQ=255)
                            gCATGATACAGACGTAGCGCTACACCACTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGcc                                                                                                                    <  3:566932/138‑1 (MQ=255)
                            gCATGATACAGACGTAGCGCTACACCACTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAg                                                                                                                  <  2:789713/140‑1 (MQ=255)
                            gCATGATACAGACGTAGCGCTACACCACTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAg                                                                                                                  <  5:226054/140‑1 (MQ=255)
                            gCATGATACAGACGTAGCGCTACACCACTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAg                                                                                                                  <  5:1120147/140‑1 (MQ=255)
                            gCATGATACAGACGTAGCGCTACACCACTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAg                                                                                                                  <  2:355530/140‑1 (MQ=255)
                            gCATGATACAGACGTAGCGCTACACCACTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTAACCTGTAATATTTATGCCac                                                                                                                  >  6:876049/1‑139 (MQ=255)
                            gCATGATACAGACGTAGCGCTACACCACTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGGTATTTTTAAACTTTTAAATTTATGCCAg                                                                                                                  >  7:861377/1‑140 (MQ=255)
                            gCATGATACAGACGTAGCGCTACACCACTGGCGATATCCCCAGCACCGCGAATCACCACCAGCTGGTGTCAGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAg                                                                                                                  <  8:828413/140‑1 (MQ=255)
                            cagtGATACAGACGTAGCGCTACACCACTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAg                                                                                                                  <  5:145695/137‑1 (MQ=255)
                              aTGATACAGACGTAGCGCTACACCACTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAg                                                                                                                  <  3:915769/138‑1 (MQ=255)
                              aTGATACAGACGTAGCGCTACACCACTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAGCa                                                                                                                <  8:1016708/140‑1 (MQ=255)
                                 aTACAGACGTAGCGCTACACCACTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAGCAc                                                                                                               <  4:348374/138‑1 (MQ=255)
                                  tACAGACGTAGCGCTACACCACTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAGCACGCt                                                                                                           <  6:1077127/140‑1 (MQ=255)
                                    cAGACGTAGCGCTACACCACTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAGCACGCTaa                                                                                                         <  7:131563/140‑1 (MQ=255)
                                    cAGACGTAGCGCTACACCACTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAGCACGCTaa                                                                                                         <  7:57005/140‑1 (MQ=255)
                                    cAGACGTAGCGCTACACCACTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGGTATTTTTCACCGGTAATATTTATGGCAGCACGCTaa                                                                                                         >  6:828752/1‑140 (MQ=255)
                                       aCGTAGCGCTACACCACTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGGTAATTTTAACCTGTAATATTTATGGCAGCACGCGAATgg                                                                                                      >  1:1312498/1‑140 (MQ=255)
                                       aCGTAGCGCTACACCACTGGCGATATCCCCCGAACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCCCCCGGCGCGTTTTGTTTGTCCAGGAATATTTGTGACAGCACGCTaaaag                                                                                                      >  2:851031/1‑137 (MQ=255)
                                         gTAGCGCTACACCACTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTGAATTTTATCCAGGAAAATTTCTTACACCACCCCAATgg                                                                                                      >  3:1177568/1‑138 (MQ=255)
                                         gTAGCGCTACACCACTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGGTATTTTTAGCAGGTATTATTTTTGCCAGCCCGCTAAttg                                                                                                      >  4:194624/1‑136 (MQ=255)
                                          tAGCGCTACACCACTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGGTATTTTTGACCTGTAATATTTCTGCCAGCACGCTAATGGCgg                                                                                                   >  2:981650/1‑139 (MQ=255)
                                          tAGCGCTACACCACTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGGGTTTTTTTAACAGTAATATTTATTCCAGCACCCTAAttggtt                                                                                                   >  6:145695/1‑135 (MQ=255)
                                           aGCGCTACACCACTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAGCACGCTAATGGCGAt                                                                                                  <  2:563263/140‑1 (MQ=255)
                                           aGCGCTACACCACTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAGCACGCTAATGGCGAt                                                                                                  <  1:882186/140‑1 (MQ=255)
                                            gcgcTACACCACTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCACCCCACCCGGCGCGGGATTTTTGCCCGGTAATATTTGAGCCAGCCCGCAATTggggctc                                                                                                 >  7:81355/1‑135 (MQ=255)
                                            gcgcTACACCACTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCCCCCGGCGCGTGATTTTTGACCGGTTATATTTATGACCGCCCGCGAATgggggtc                                                                                                 >  7:370344/1‑135 (MQ=255)
                                            gcgcTACACCACTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTTTTTTTGACCTGTTATATTTATGCCAGCACGCGAAAGGCGAtc                                                                                                 >  7:482099/1‑140 (MQ=255)
                                            gcgcTACACCACTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTTTTTTTCACCTGTAATATTTCTGCCAGCAAGCTAATGGGGAtc                                                                                                 >  8:802459/1‑140 (MQ=255)
                                            gcgcTACACCACTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGGTATTTTTGACCAGTTATAATTCTTCCAGACCGCATAATgggaga                                                                                                  >  5:89583/1‑136 (MQ=255)
                                             cgcTACACCACTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGGTATTTTTCACCTGTAATATTTATGCCAGCACGCTAATGGGGaa                                                                                                  >  4:1039222/1‑137 (MQ=255)
                                              gcTACACCACTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGGGTTTTTTAAACTGGTATATTTATGCCCGCCCGCTAATGGGGAta                                                                                                 >  3:877791/1‑137 (MQ=255)
                                                  caccacTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTTTTTTTAACCTGTAATATTTATGGCAGCCCCCGAATGGCGATCACCAgc                                                                                           >  6:107361/1‑140 (MQ=255)
                                                    ccacTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTATGCCAGCACGCTAATGGCGAGATCCTgc                                                                                           >  4:558424/1‑138 (MQ=255)
                                                       cTGGCGATATCCCCCGCACCGCGAATCCCCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGGGTTTTTTTACCGGGAATAATTTTGACACCACCGGACTAGCGAGATCCTgcggggg                                                                                      >  1:1083914/1‑137 (MQ=255)
                                                       cTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTGACCTGTAATATTTAGGCCCGCACGCTAATGGCGATCTCCAgcggcgg                                                                                      >  2:1092782/1‑139 (MQ=255)
                                                       cTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGGTTTTTTTAACCTGTAATATTTATTCCAGCACCCGAAATGCGATCACCAgcggcgg                                                                                      >  8:57591/1‑139 (MQ=255)
                                                       cTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGGTATTTTTAACCTGTAATATTTATGCCAGCACGCTAAATGCGATCACCAGCGGCGt                                                                                      >  6:999231/1‑140 (MQ=255)
                                                       cTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCAACCCGGCGCGGTATTTTTGAACTGGAATATTTATGACAGCCCGCTAATGGCGCTCACCGgaggcgg                                                                                      >  2:544406/1‑139 (MQ=255)
                                                         ggCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTGACCTGGTATATTTCTTCCAGCCCGCTAATGGCTATCTCCAgcggcgg                                                                                      >  4:506742/1‑137 (MQ=255)
                                                         ggCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCACGTTATTTTTCACCTGTAATATTTCTGCCAGCACGCTAATGGCGATCTCCTGCGGCGt                                                                                      >  4:805872/1‑138 (MQ=255)
                                                          gCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAGCACGCTAATGGCGATCTCCTGCGGCGTTTc                                                                                   <  1:699205/140‑1 (MQ=255)
                                                          gCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAGCACGCTAATGGCGATCTCCTGCGGCGTTTc                                                                                   <  2:1376621/140‑1 (MQ=255)
                                                          gCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAGCACGCTAATGGCGATCTCCTGCGGCGTTTc                                                                                   <  6:431983/140‑1 (MQ=255)
                                                          gCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAGCACGCTAATGGCGATCTCCTGCGGCGTTTc                                                                                   <  1:558865/140‑1 (MQ=255)
                                                          gCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACACGGAGCGTTATTTTTCACCTGTAATATTTCTGCCAGCACGCTAATGGCGATCTCCTGCGGCGTTTc                                                                                   <  5:521426/140‑1 (MQ=255)
                                                          gCGATATCCCCCGCAACGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAGCACGCTAATGGCGATCTCCTGCGGCGtt                                                                                     <  3:503755/138‑1 (MQ=255)
                                                            gATCTCCCCCGCACCGCGCATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCCGCACGCTAATGGCGATCTCCTGCGGCGTTTCCg                                                                                 <  2:1231445/140‑1 (MQ=255)
                                                            gATATCCCCCGCACCGCGAATCACCCCCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAGCACGCTAATGGCGATCTCCTGCGGCGTTTCCg                                                                                 <  8:26216/140‑1 (MQ=255)
                                                            gATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGACGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTGATTTTTCACCTGTAATCTTTCTGCCAGCACGCTAATGGCGATCTCCTGCGGCGTTTCCg                                                                                 <  8:941276/140‑1 (MQ=255)
                                                            gATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGCAGGCTTCATCATCAACCCACCAGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAGCACGCTAATGTCGATCTCCTGCGGCGTTTc                                                                                   <  4:25298/138‑1 (MQ=255)
                                                            gATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCTGCACGCTAATGGCGATCTCCTGCGGCGTTTCCg                                                                                 <  6:962677/140‑1 (MQ=255)
                                                            gATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAGCACGCTAATGGCGATCTCCTGCGGCGTTTc                                                                                   <  3:1132360/138‑1 (MQ=255)
                                                            gATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAGCACGCTAATGGCGATCTCCTGCGGCGTTTCCg                                                                                 <  1:639723/140‑1 (MQ=255)
                                                            gATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTGCATCATCAACCCACACGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAGCACGCTAATGGCGATCTCCTGCGGCGTTTCCg                                                                                 <  6:719589/140‑1 (MQ=255)
                                                               atCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTTTTTTTAACCTGTAATATTTCTGGCCGCACGCTAATGGCGATCTCCTGCGGCGGTTACGccc                                                                              >  8:895899/1‑138 (MQ=255)
                                                               atCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGGTATTTTTAACCGGGAATATTTGTGCCAGCCCGCAAATAGCGCTCAACGGCGGCGTTTACGGAc                                                                              >  6:1125544/1‑140 (MQ=255)
                                                                 cccccGCACCGCGAATCCCCACCAGCTGGTGTCCGGAAGGATGAGCAGGCATCAGCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATATCTGCCAGCACGCTAATGGCGATCTCCTGCGGCGTTTCCGCACCt                                                                            <  8:1143435/140‑1 (MQ=255)
                                                                 cccccGCACCGCGAAGCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAGCACGCTAATGGCGATCTCCTGCGGCGTTTCCGCACCt                                                                            <  1:98155/140‑1 (MQ=255)
                                                                 cccccGCACCGCGAAGCACCACCAGCTGGAGGCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAGCACGCTAATGGCGATCTCCTGCGGCGTTTCCGCACCt                                                                            <  2:1216270/140‑1 (MQ=255)
                                                                 cccccGCAACGCGTATCCCCAACAGCTTTGTTGCGGTAGGATGTGAATGAGTCATAATCAACCCAACCGGCGCGTTATTTTTCAACTTTAATTTTTCTGCCAGCACGCTAATGGCGATCTCCTGCGGCGTTTCCGCAGCt                                                                            <  2:315257/140‑1 (MQ=255)
                                                                    ccGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAGCACGCTAATGGCGATCTCCTGCGGCGTTTCCGCACCTa                                                                           <  4:598509/138‑1 (MQ=255)
                                                                    ccGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAGCACGCTAATGGCGATCTCCTGCGGCGTTTCCGCACCTAtg                                                                         <  5:497170/140‑1 (MQ=255)
                                                                         ccGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAGCACGCTAATGGCGATCTCCTGCGGCGTTTCCGCACCTATGTTGTa                                                                    <  7:266563/140‑1 (MQ=255)
                                                                         ccGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAGCACGCTAATGGCGATCTCCTGCGGCGTTTCCGCACCTATGTTGTa                                                                    <  8:1179164/140‑1 (MQ=255)
                                                                            cgAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTAACCTGTAATATTTATGCCAGCACGCGAATTGCGATCTCCTGCGGGGTTTACGCCCCTATGTTGTaa                                                                   >  4:1293228/1‑138 (MQ=255)
                                                                            cgAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGGTTTTTTTTAACAGTAATATTTATGCCAGCCCGCTAATGGGGGTATCCAGCGGCGTTTACGCCCCTAtgttttta                                                                   >  3:678431/1‑134 (MQ=255)
                                                                              aaTCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTAACCAGTAATATTTATGCCAGCACGCTAAAGGGGATCACCAGCGGCGGTTACGCACCCATTTTTTaaac                                                                 >  4:551615/1‑136 (MQ=255)
                                                                                  accaccAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAGCACGCTAATGGCGATCTCCTGCGGCGTTTCCGCACCTATGTTGTAACCAACggg                                                           <  7:910500/140‑1 (MQ=255)
                                                                                    caccaGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAGCACGCTAATGGCGATCTCCTGCGGCGTTTCCGCACCTATGTTGTAACCAACGGgcg                                                         <  8:418261/140‑1 (MQ=255)
                                                                                    caccaGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTAACCTGTAATATTTATGCCAGCCCCCTAATGGCGATCTCCAGCGGCGGTTACCCACCAATTTTGTAACCAACcgggg                                                         >  1:31370/1‑138 (MQ=255)
                                                                                           tGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGGTATTTTTAACCAGTAATATTTATGGCAGCACGCGTATGGCGATCACCAGCGGCCGTTCCGCCCCAATGTTTGAAACCACGGGCGCATGt                                                    >  3:1288873/1‑138 (MQ=255)
                                                                                                tCCGGAAGGATGAGAAGTCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAGCACGCTAATGGCGATCTCCTGCGGCGTTTCCGCACCTATGTTGTAACCAACGGGCGCATGTAAACggg                                             <  2:1369020/140‑1 (MQ=255)
                                                                                                tCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTTATAATTCTGGCCGCACCCTAATGGCGATCTCCTGCGGCCTTTCCGCCCCTATGTTGTAACAAAAGGGCGCAAGTAAACg                                               >  3:738311/1‑138 (MQ=255)
                                                                                                tCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATGTTTCTCCTGTAATATTTCTGCCAGCTCGCTAATGGCGATCTCCTGCGGCGTTTCCGCACCTATGTTGTAACCAACGGGCGCATGTAAACggg                                             <  5:635308/140‑1 (MQ=255)
                                                                                                  cGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAGCACGCTAATGGCGATCTCCTGCGGCGTTTCCGCACCTATGTTGTAACCAACGGGCGCATGTAAACggg                                             <  4:144960/138‑1 (MQ=255)
                                                                                                  cGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAGCACGCTAATGGCGATCTCCTGCGGCGTTTCCGCACCTATGTTGTAACCAACGGGCGCATGTAAACggg                                             <  3:281217/138‑1 (MQ=255)
                                                                                                    gAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAGCACGCTAATGGCGATCTCCTGCGGCGTTTCCGCACCTATGTTGTAACCAACGGGCGCATGTAAACGGGCa                                           <  3:1332221/138‑1 (MQ=255)
                                                                                                     aaGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAGCACGCTAATGGCGATCTCCTGCGGCGTTTCCGCACCTATGTTGTAACCAACGGGCGCATGTAAACGGGCAata                                        <  6:1210409/140‑1 (MQ=255)
                                                                                                       ggATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAGCACGCTAATGGCGATCTCCTGCGGCGTTTCCGCACCTATGTTGTAACCAACGGGCGCATGTAAACGGGCAATATg                                      <  1:142178/140‑1 (MQ=255)
                                                                                                       ggATGAGAAGGCTTCATCATCAACCCACCCGGCGAGTTATTTTTCACCTGTAATATTTCTGCCAGCACGCTAATGGCGATCTCCTGCGGCGTTTCCGCACCTATGTTGTAACCAACGGGCGCATGTAAACGGGCAATATg                                      <  6:296193/140‑1 (MQ=255)
                                                                                                          tGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAGCACGCTAATGGCGATCTCCTGCGGCGTTTCCGCACCTATGTTGTAACCAACGGGCGCATGTAAACGGGCAATATGt                                     <  4:690262/138‑1 (MQ=255)
                                                                                                          tGAGAAGGCTTCATCATCAACCCACCCGGCGCGGTATTTTTGACCGGTAATATTTATGCCCCCACGCCAAAGGCGCTGTCCTGCGGCGTTTACGCCCCAATTTTTTTACCAACGGGCGCATGTTAACGGGCAATAttttt                                   >  1:35795/1‑139 (MQ=255)
                                                                                                           gagaAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAGCACGCTAATGGCGATCTCCTGCGGCGTTTCCGCACCTATGTTGTAACCAACGGGCGCATGTAAACGGGCAATATGttct                                  <  2:401546/140‑1 (MQ=255)
                                                                                                            agaAGGCTTCATCATCAACCCACCCGGCGCGGTATTTTTAACCAGTAATATTTCTTCCACCACCCAAAATGCGCAATCCTTCGGCCTTTGCCCACCAATGTTGTAAACAAAGGGCGCGAGTAAAACGGCAATAtttta                                   >  4:886893/1‑137 (MQ=255)
                                                                                                                ggCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAGCACGCTAATGGCGATCTCCTGCGGCGTTTCCGCACCTATGTTGTAACCAACGGGCGCATGTAAACGGGCAATATGTTCTTCa                               <  3:399277/138‑1 (MQ=255)
                                                                                                                ggCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAGCACGCTAATGGCGATCTCCTGCGGCGTTTCCGCACCTATGTTGTAACCAACGGGCGCATGTAAACGGGCAATATGTTCTTCAGc                             <  1:926358/140‑1 (MQ=255)
                                                                                                                   ttCATCATCAACCCACCCGGCGCGTTATTTTTCACCAGTAATATTTATTCCAGCACGCAAATGGCGATCTACTGCCGCGTTTCCGCACCTATGTTGTAACCAACCGGCGCAAGTAAACGGGGAATATGTTCTTCATCCAc                          >  2:1114907/1‑140 (MQ=255)
                                                                                                                              cccacccGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAGCACGCTAATGGCGATCTCCTGCGGCGTTTCCGCACCTATGTTGTAACCAACGGGCGCATGTAAACGGGCAATATGTTCTTCAGCCACGCCTTTCTCAc               <  1:540861/140‑1 (MQ=255)
                                                                                                                                       cgcgTTATTTTTCACCTGTAATATTTCTGCCAGCACGCTAATGGCGATCTCCTGCGGCGTTTCCGCACCTATGTTGTAACCAACGGGCGCATGTAAACGGGCAATATGTGCTTCAGCCACGCCTTTCTCACGCAATTg        <  4:947895/138‑1 (MQ=255)
                                                                                                                                           ttATTTTTCACCTGTAATATTTCTGCCAGCACGCTAATGGCGATCTCCTGCGGCGTTTCCGCACCTATGTTGTAACCAACGGGCGCATGTAAACGGGCAATATGTTCTTCAGCCACGCCTTTCTCACGCAATTGGCGCAg  <  5:856139/140‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                            
CACTTCCAACATGATCACTTTAAAACCCGCATGATACAGACGTAGCGCTACACCACTGGCGATATCCCCCGCACCGCGAATCACCACCAGCTGGTGTCCGGAAGGATGAGAAGGCTTCATCATCAACCCACCCGGCGCGTTATTTTTCACCTGTAATATTTCTGCCAGCACGCTAATGGCGATCTCCTGCGGCGTTTCCGCACCTATGTTGTAACCAACGGGCGCATGTAAACGGGCAATATGTTCTTCAGCCACGCCTTTCTCACGCAATTGGCGCAG  >  NZ_CP022686/1633707‑1633985

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: