Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP025268 4,579,872 T→C M1M (ATG→GTG) †
*460* (TAA→TAG
mcrC ←
mcrB ←
5‑methylcytosine‑specific restriction endonuclease system specificity protein McrC
5‑methylcytosine‑specific restriction endonuclease subunit McrB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0252684,579,8720TC99.6% 878.5 / NA 264M1M (ATG→GTG) †
*460* (TAA→TAG
mcrC
mcrB
5‑methylcytosine‑specific restriction endonuclease system specificity protein McrC
5‑methylcytosine‑specific restriction endonuclease subunit McrB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (179/84);  minor base A (0/1);  total (179/85)
Fisher's exact test for biased strand distribution p-value = 3.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 264 total aligned reads displayed.

AATACATACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATC  >  NZ_CP025268/4579736‑4580006
                                                                                                                                        |                                                                                                                                       
   tcatacCCCATGATATCAAGTAGATTGTTACCGGCTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCACATGCATAGGTAAGCATGTAATAGATATTACGGACCCGTCTCACGGGCTGTTCCACTATGAg                                                                                                                                   <  2:3075935/139‑1 (MQ=255)
   acattcACCAGGATATAAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATTACGTGCTGTTCCACTATGAg                                                                                                                                   <  2:2000598/140‑1 (MQ=255)
   acatccCCCATGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAg                                                                                                                                   <  1:3932476/140‑1 (MQ=255)
   acatacCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTGCAAGGTGTGCCTGCTTAATTTCCTGTCAATAACCGCATGCGTAGGTAAGCATGTAATAGATATTACGGGCAGGGATCACGGGCTGTTCCACTATGAg                                                                                                                                   <  2:1373979/140‑1 (MQ=255)
       atctcAGTATATCAATAAGATTGTTACCGGGTATAGCTTCAAGGTTTGCCTGCTTAATTTACTTTAAATAACCCCACGCATAGGTAAGCATGTAATAGTTATTACGGACAGGTATCACTTGCTGTTTCACTATGAGTcac                                                                                                                               <  2:2602196/136‑3 (MQ=255)
       acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTccc                                                                                                                               >  1:2953751/1‑140 (MQ=255)
       acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTccc                                                                                                                               >  1:307826/1‑140 (MQ=255)
       acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTccc                                                                                                                               >  1:2698199/1‑140 (MQ=255)
       acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTccc                                                                                                                               >  2:1475506/1‑140 (MQ=255)
       acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTccc                                                                                                                               <  2:5076550/140‑1 (MQ=255)
       acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTccc                                                                                                                               >  1:2152464/1‑140 (MQ=255)
       acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTccc                                                                                                                               >  1:4151960/1‑140 (MQ=255)
       acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTccc                                                                                                                               >  1:4870156/1‑140 (MQ=255)
       acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTccc                                                                                                                               >  1:712079/1‑140 (MQ=255)
       acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTccc                                                                                                                               >  2:412441/1‑140 (MQ=255)
       acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTccc                                                                                                                               >  1:5310924/1‑140 (MQ=255)
       acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAAGAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTccc                                                                                                                               >  2:2030106/1‑140 (MQ=255)
       acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAAGAGATATTACGGACAGGTATCACGGGCTGTGCCACTATGAGTccc                                                                                                                               >  2:536722/1‑140 (MQ=255)
       acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGGAAGCATGTAATAGATATTACGGACAGGTATCAGGGGCTGTTCCACTAGTAGTccc                                                                                                                               >  2:4274535/1‑140 (MQ=255)
        ccccAGGATATCAAGAATATTGTTACAGGGTATTGCTTAAAGGTTTGCCTGCTTAATTTCATGTAAATACCCCCATGAATAGGTAAGCATGTAATACATATTACGGACAAGTATAACGGTCTGTTCCACTATGACTcccc                                                                                                                              <  2:2218646/140‑1 (MQ=255)
            aGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTATGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCGCCgaat                                                                                                                          >  1:4668178/1‑140 (MQ=255)
            aGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCtaat                                                                                                                          >  2:2262892/1‑140 (MQ=255)
            aGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCtaat                                                                                                                          >  1:4288933/1‑140 (MQ=255)
             ggATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCtaata                                                                                                                         >  2:3729703/1‑140 (MQ=255)
              gATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCtaataa                                                                                                                        >  1:4130056/1‑140 (MQ=255)
                   cAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATttgt                                                                                                                   >  1:3454760/1‑140 (MQ=255)
                   cAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATttgt                                                                                                                   >  1:804320/1‑140 (MQ=255)
                   cAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATttgt                                                                                                                   >  1:4229977/1‑140 (MQ=255)
                   cAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGGTCCACTATGAGTCCCCTAATAATtggt                                                                                                                   >  2:2716734/1‑140 (MQ=255)
                   cAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGATCACGGGCTGTTCCACTATGAGTCCCCTAATAATttgtt                                                                                                                  >  2:313513/1‑140 (MQ=255)
                       aagaTTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGt                                                                                                               >  1:3582367/1‑140 (MQ=255)
                            tGTTACCGTGTACGGCGTCAAGGGTTGCCGGCGTAATTGCCTGTAAATAACGCCATGCGCAGGGAGGCATGGAATAGATATTACGGACAGGGATCAGGGGCTGTTCTACTATGATTCCCCTAATAGTTGGGGTGTCCAtt                                                                                                          >  2:366215/1‑140 (MQ=255)
                            tGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCAtt                                                                                                          >  1:2921968/1‑140 (MQ=255)
                            tGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCAtt                                                                                                          >  1:3652378/1‑140 (MQ=255)
                               tACCGGGTATGTCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGAAAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCt                                                                                                       <  2:49173/140‑1 (MQ=255)
                                aCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtg                                                                                                      <  2:270449/140‑1 (MQ=255)
                                 ccGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCGGTAAATAACCCCATGCAGCGGTAAGCATGTAATAGATATTATGGACAGGTATCACGGGTTGTTCCACTATGAGTCCCCTAATAATTTGTGGGTCCATTTCtgt                                                                                                     >  2:4835951/1‑140 (MQ=255)
                                  cGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt                                                                                                    <  1:1991653/140‑1 (MQ=255)
                                        tGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAATTAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTtat                                                                                              >  2:4797435/1‑140 (MQ=255)
                                        tGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTtat                                                                                              >  1:2642682/1‑140 (MQ=255)
                                              cAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGt                                                                                       >  2:5322258/1‑140 (MQ=255)
                                               aaGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTCATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGtc                                                                                      >  2:846075/1‑140 (MQ=255)
                                                   tttGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACGATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGtcatca                                                                                  <  1:3608572/140‑1 (MQ=255)
                                                     tGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCaaa                                                                                >  2:969034/1‑140 (MQ=255)
                                                     tGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCaaa                                                                                >  1:3037459/1‑140 (MQ=255)
                                                     tGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATATCTGTTGTTTATAGGGGTCATCaaa                                                                                >  2:3706534/1‑140 (MQ=255)
                                                          gCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAat                                                                           <  1:525653/140‑1 (MQ=255)
                                                           cTTAATTTCCTGTAAAGAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGATAGGTATCGCGGGCTGTTTCACTATGAGTCTCCTAATAATTTGTTGGGCCATTTCTGTTGGTTATATGGGTGATCAAGGGAatg                                                                          >  2:4668550/1‑139 (MQ=255)
                                                            ttAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAatat                                                                         <  1:3330365/140‑1 (MQ=255)
                                                                    cTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAg                                                                 <  2:3383458/140‑1 (MQ=255)
                                                                     tGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTGATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGGTCATCAGAGAAATATTCTTCGAg                                                                 >  1:4471072/1‑140 (MQ=255)
                                                                     tGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACAGGCCGGTACACTATGAGTCCCCTAATAATTTGTTGGTCCATTGCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGt                                                                >  2:1305238/1‑140 (MQ=255)
                                                                      gTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTa                                                               >  2:1902855/1‑140 (MQ=255)
                                                                        aaaTAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTaaa                                                             <  1:497018/140‑1 (MQ=255)
                                                                        aaaTAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTaaa                                                             <  2:3504671/140‑1 (MQ=255)
                                                                              ccccATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCg                                                       >  1:3800440/1‑140 (MQ=255)
                                                                              ccccATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATATCTGTTGTGTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCg                                                       >  2:5186804/1‑140 (MQ=255)
                                                                                 cATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGata                                                    >  1:951012/1‑140 (MQ=255)
                                                                                   tGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATc                                                  >  1:2678151/1‑140 (MQ=255)
                                                                                   tGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATc                                                  >  2:4003445/1‑140 (MQ=255)
                                                                                         ggTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGtcat                                            >  1:1380380/1‑140 (MQ=255)
                                                                                         ggTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGtcat                                            >  2:333262/1‑140 (MQ=255)
                                                                                         ggTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGtcat                                            >  1:3519142/1‑140 (MQ=255)
                                                                                          gTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGtcatc                                           >  2:5551773/1‑140 (MQ=255)
                                                                                           tAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGtcatca                                          <  1:2844689/140‑1 (MQ=255)
                                                                                             aGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATcaca                                        >  1:3398358/1‑140 (MQ=255)
                                                                                             aGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGGCCCCTAATAGTTTGTTGGTCCATTTCGGTGGTTGATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATcaca                                        >  2:586117/1‑140 (MQ=255)
                                                                                               cATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTCTAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAAt                                      >  1:3990450/1‑140 (MQ=255)
                                                                                               cATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAAt                                      >  1:5567616/1‑140 (MQ=255)
                                                                                               cATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAAt                                      <  2:804320/140‑1 (MQ=255)
                                                                                               cATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAAt                                      <  2:3586872/140‑1 (MQ=255)
                                                                                                 tGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAAttt                                    >  2:3870303/1‑140 (MQ=255)
                                                                                                   tAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCa                                  <  2:4580686/140‑1 (MQ=255)
                                                                                                    aaTAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAATGGGCGATATCCGTCATCACAATTTCAt                                 >  2:2054230/1‑140 (MQ=255)
                                                                                                    aaTAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCAt                                 >  1:2224199/1‑140 (MQ=255)
                                                                                                    aaTAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCAt                                 >  1:2251859/1‑140 (MQ=255)
                                                                                                    aaTAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCAt                                 >  1:2662567/1‑140 (MQ=255)
                                                                                                    aaTAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGGAAAGGGGCGATATCCGTCATCACAATTTCAt                                 >  1:4108379/1‑140 (MQ=255)
                                                                                                     aTAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTCTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCTTCATCACAATTTCAtt                                <  2:5477109/140‑1 (MQ=255)
                                                                                                         atatTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGaaa                                                                            >  1:1502077/1‑92 (MQ=255)
                                                                                                         atatTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGc                            >  2:4755739/1‑140 (MQ=255)
                                                                                                         atatTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGc                            >  1:1331288/1‑140 (MQ=255)
                                                                                                         atatTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGc                            >  2:24364/1‑140 (MQ=255)
                                                                                                         atatTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGc                            >  1:5468473/1‑140 (MQ=255)
                                                                                                           atTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCa                          >  1:1037538/1‑140 (MQ=255)
                                                                                                            ttACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCAt                         >  2:4218597/1‑140 (MQ=255)
                                                                                                             tACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCAtt                        >  1:2060950/1‑140 (MQ=255)
                                                                                                             tACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCAtt                        >  1:2384914/1‑140 (MQ=255)
                                                                                                             tACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCAtt                        >  2:1595070/1‑140 (MQ=255)
                                                                                                              aCGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTg                       <  2:1073788/140‑1 (MQ=255)
                                                                                                              aCGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTg                       <  1:3748262/140‑1 (MQ=255)
                                                                                                                    aGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATc                 <  1:866383/140‑1 (MQ=255)
                                                                                                                    aGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATc                 >  2:3312486/1‑140 (MQ=255)
                                                                                                                      gTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCg               >  2:3008837/1‑140 (MQ=255)
                                                                                                                            cGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAgg         >  1:3176913/1‑140 (MQ=255)
                                                                                                                              ggCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGtt       >  2:3884899/1‑139 (MQ=255)
                                                                                                                                 tGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCa    <  1:2453448/140‑1 (MQ=255)
                                                                                                                                   ttCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTAGCCGGAGAGGTGCCATc  >  2:3451110/1‑140 (MQ=255)
                                                                                                                                   ttCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGCAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATc  >  2:4942185/1‑140 (MQ=255)
                                                                                                                                   ttCCACTATGAGTCCCCAAATAATTTGTTGGTCCATTTCTGGTGTGTATAGGGGTGATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCAGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCagc  >  2:1224485/1‑138 (MQ=255)
                                                                                                                                        |                                                                                                                                       
AATACATACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATC  >  NZ_CP025268/4579736‑4580006

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