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.9% 5106.9 / NA 1461M1M (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 (739/721);  minor base G (0/1);  total (739/722)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 1462 total aligned reads displayed.

TAATACATACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATC  >  NZ_CP025268/4579735‑4580006
                                                                                                                                         |                                                                                                                                      
  atacatacCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATg                                                                                                                                    >  2:854458/1‑140 (MQ=255)
   tacatacCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGa                                                                                                                                   >  1:10213083/1‑140 (MQ=255)
    acagacCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAg                                                                                                                                  <  2:6617198/140‑1 (MQ=255)
      atacCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTc                                                                                                                                <  2:6887148/140‑1 (MQ=255)
        acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAg                                                                                                                                  <  2:5434033/136‑1 (MQ=255)
        acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTccc                                                                                                                              >  2:10697322/1‑140 (MQ=255)
        acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTccc                                                                                                                              >  1:1671689/1‑140 (MQ=255)
        acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTccc                                                                                                                              >  2:5121700/1‑140 (MQ=255)
        acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTccc                                                                                                                              >  1:7190695/1‑140 (MQ=255)
        acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTccc                                                                                                                              >  2:7835910/1‑140 (MQ=255)
        acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTccc                                                                                                                              >  2:2659735/1‑140 (MQ=255)
        acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTccc                                                                                                                              >  2:3790993/1‑140 (MQ=255)
        acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAAGAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTccc                                                                                                                              <  1:10007268/140‑1 (MQ=255)
        acCCCAGGATATCAAGAAGATTGTTACCGGGTAAGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTccc                                                                                                                              >  2:5887969/1‑140 (MQ=255)
         ccccAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTcccc                                                                                                                             <  2:3424287/140‑1 (MQ=255)
             aGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCtaat                                                                                                                         >  1:4049629/1‑140 (MQ=255)
                  atCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAAttt                                                                                                                    >  1:11380905/1‑140 (MQ=255)
                   tCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATttg                                                                                                                   >  2:10968171/1‑140 (MQ=255)
                    cAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATttgt                                                                                                                  >  1:5853772/1‑140 (MQ=255)
                    cAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATttgt                                                                                                                  >  2:7429201/1‑140 (MQ=255)
                    cAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATttgt                                                                                                                  >  1:855658/1‑140 (MQ=255)
                      agaagaTTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATttgttg                                                                                                                >  1:5114642/1‑140 (MQ=255)
                      agaagaTTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATAACGGGCTGTTCCACTATGAGTCCCCTAATAATttgttg                                                                                                                >  2:2113157/1‑140 (MQ=255)
                         agaTTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTc                                                                                                             >  2:87394/1‑140 (MQ=255)
                           aTTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCa                                                                                                           >  2:7266237/1‑140 (MQ=255)
                           aTTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCa                                                                                                           >  2:2828203/1‑140 (MQ=255)
                             tGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCAtt                                                                                                         >  2:7053330/1‑140 (MQ=255)
                             tGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCAtt                                                                                                         >  2:11156818/1‑140 (MQ=255)
                              gTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCAttt                                                                                                        <  2:10351270/140‑1 (MQ=255)
                              gTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCAttt                                                                                                        <  2:8089261/140‑1 (MQ=255)
                                tACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCt                                                                                                      >  1:8908965/1‑140 (MQ=255)
                                tACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCt                                                                                                      <  1:5927089/140‑1 (MQ=255)
                                 aCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtg                                                                                                     <  1:655925/140‑1 (MQ=255)
                                 aCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtg                                                                                                     <  2:1381380/140‑1 (MQ=255)
                                 aCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtg                                                                                                     <  2:2538759/140‑1 (MQ=255)
                                  ccGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgt                                                                                                    <  1:10449366/140‑1 (MQ=255)
                                   cGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt                                                                                                   <  1:4564827/140‑1 (MQ=255)
                                   cGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt                                                                                                   <  2:8908964/140‑1 (MQ=255)
                                   cGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt                                                                                                   <  1:2659736/140‑1 (MQ=255)
                                     ggTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt                                                                                                   >  2:10525336/1‑138 (MQ=255)
                                       tATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGttt                                                                                               >  1:3255292/1‑140 (MQ=255)
                                         tGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTtat                                                                                             >  1:6814091/1‑140 (MQ=255)
                                         tGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTtat                                                                                             >  2:9504973/1‑140 (MQ=255)
                                         tGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCACCTAATAATTTGTTGGTCCATTTCTGTTGTTtat                                                                                             >  1:9699667/1‑140 (MQ=255)
                                          ggCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAg                                                                                                                                  <  2:9080421/102‑1 (MQ=255)
                                          ggCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAg                                                                                                                                  >  1:9080422/1‑102 (MQ=255)
                                          ggCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTtata                                                                                            <  1:4338674/140‑1 (MQ=255)
                                                aaGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGtc                                                                                      >  2:9685597/1‑140 (MQ=255)
                                                    tttGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGtcatca                                                                                  <  1:11208575/140‑1 (MQ=255)
                                                    tttGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGtcatca                                                                                  <  1:2237475/140‑1 (MQ=255)
                                                    tttGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGtcatca                                                                                  <  2:6098678/140‑1 (MQ=255)
                                                     ttGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCaa                                                                                 >  1:4711613/1‑140 (MQ=255)
                                                                aTTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATAttct                                                                      <  1:3082911/140‑1 (MQ=255)
                                                                aTTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATAttct                                                                      <  1:3815962/140‑1 (MQ=255)
                                                                    ccTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGa                                                                  >  2:11469967/1‑140 (MQ=255)
                                                                       gTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTa                                                               <  2:3211036/140‑1 (MQ=255)
                                                                         aaaTAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTaaa                                                             <  2:3022893/140‑1 (MQ=255)
                                                                             aaCCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGGCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAgggg                                                         <  2:9218792/140‑1 (MQ=255)
                                                                               ccccATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCg                                                       >  1:3546874/1‑140 (MQ=255)
                                                                                  cctGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGata                                                    <  1:7819100/138‑1 (MQ=255)
                                                                                  cATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGata                                                    >  1:1265742/1‑140 (MQ=255)
                                                                                    tGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTgg                                                                                                               <  1:6975806/79‑1 (MQ=255)
                                                                                        tAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGtc                                              <  1:5288896/140‑1 (MQ=255)
                                                                                          ggTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTCTAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGtcat                                            <  1:10609381/140‑1 (MQ=255)
                                                                                          ggTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGtcat                                            <  1:9331809/140‑1 (MQ=255)
                                                                                          ggTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGtcat                                            <  2:3611410/140‑1 (MQ=255)
                                                                                              aGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATcaca                                        >  2:2360230/1‑140 (MQ=255)
                                                                                              aGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATcaca                                        >  2:4995290/1‑140 (MQ=255)
                                                                                              aGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATcaca                                        >  2:4022371/1‑140 (MQ=255)
                                                                                                cATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAAt                                      >  2:7690144/1‑140 (MQ=255)
                                                                                                cATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAAt                                      <  1:10889503/140‑1 (MQ=255)
                                                                                                cATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAAt                                      <  1:9483592/140‑1 (MQ=255)
                                                                                                  tGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAAttt                                    <  2:4591461/140‑1 (MQ=255)
                                                                                                   gTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAggg                                                                                         <  1:2526329/86‑1 (MQ=255)
                                                                                                     aaTAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt                                                                                                   <  2:10170410/74‑1 (MQ=255)
                                                                                                     aaTAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCAt                                 >  1:1429219/1‑140 (MQ=255)
                                                                                                      aTAGATATTACGGACAGGTATCACGTGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCAtt                                <  2:4776901/140‑1 (MQ=255)
                                                                                                      aTAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCAtt                                >  1:9887408/1‑140 (MQ=255)
                                                                                                      aTAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCAtt                                <  1:2906161/140‑1 (MQ=255)
                                                                                                      aTAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTACATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCAtt                                <  2:1916097/140‑1 (MQ=255)
                                                                                                          atatTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTa                               >  1:8304875/1‑137 (MQ=255)
                                                                                                          atatTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGc                            >  2:5242297/1‑140 (MQ=255)
                                                                                                           tatTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGcc                           <  2:824851/140‑1 (MQ=255)
                                                                                                            atTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCAtt                                <  1:7615503/134‑1 (MQ=255)
                                                                                                              tACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCGCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGGCATCAAAGAATTATTCTTCGAGTAAAGGGGCGATATCCGTCATCACTATTTCAGTAAGCCAtt                        >  1:8703103/1‑140 (MQ=255)
                                                                                                              tACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCAtt                        >  2:1538202/1‑140 (MQ=255)
                                                                                                              tACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCAtt                        >  1:7395521/1‑140 (MQ=255)
                                                                                                              tACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCAtt                        >  1:813177/1‑140 (MQ=255)
                                                                                                               aCGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTg                       <  1:7982921/140‑1 (MQ=255)
                                                                                                               aCGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTg                       <  2:1713134/140‑1 (MQ=255)
                                                                                                               aCGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTATATAGGGTTAATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTg                       >  2:4148267/1‑140 (MQ=255)
                                                                                                                cGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGc                      <  1:1818318/140‑1 (MQ=255)
                                                                                                                          tCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGgag            <  1:3436505/140‑1 (MQ=255)
                                                                                                                                 cTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGcc     >  2:4373350/1‑140 (MQ=255)
                                                                                                                                 cTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGcc     >  1:5606845/1‑140 (MQ=255)
                                                                                                                                  tGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAACTATTCTTCGAGTAAAGGTGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCa    <  2:636761/140‑1 (MQ=255)
                                                                                                                                  tGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCa    >  2:5619291/1‑140 (MQ=255)
                                                                                                                                    ttCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATc  >  1:2023939/1‑140 (MQ=255)
                                                                                                                                    ttCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATc  >  1:499618/1‑140 (MQ=255)
                                                                                                                                    ttCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTAGAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATc  >  1:6241455/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                      
TAATACATACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATC  >  NZ_CP025268/4579735‑4580006

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: