Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,567,114 G→A 13.8% S221F (TCT→TTT)  mcrC ← 5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,567,1140GA13.8% 681.2 / 77.1 231S221F (TCT→TTT) mcrC5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrC
Reads supporting (aligned to +/- strand):  ref base G (107/92);  new base A (13/19);  total (120/111)
Fisher's exact test for biased strand distribution p-value = 1.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 236 total aligned reads displayed.

TCCCATTCGTCGTGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATG  >  CP009273/4566966‑4567263
                                                                                                                                                    |                                                                                                                                                     
tCCCATTCGTCGTGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAaga                                                                                                                                                      <  1:588433/150‑1 (MQ=255)
tCCCATTCGTCGTGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAaga                                                                                                                                                      <  1:2946422/150‑1 (MQ=255)
tCCCATTCGTCGTGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAaga                                                                                                                                                      <  2:2851153/150‑1 (MQ=255)
   cATTCGTCGTGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCg                                                                                                                                                   >  2:101761/1‑150 (MQ=255)
   cATTCGTCGTGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCg                                                                                                                                                   >  2:825716/1‑150 (MQ=255)
   cATTCGTCGTGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCg                                                                                                                                                   >  2:101913/1‑150 (MQ=255)
       cgtcgtGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGtt                                                                                                                                               <  2:563358/150‑1 (MQ=255)
       cgtcgtGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGtt                                                                                                                                               >  2:2278350/1‑150 (MQ=255)
         tcgtGAAAAAATTCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAAAGCGGGTtgt                                                                                                                                             <  2:3502342/150‑1 (MQ=255)
               aaaaaTGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGc                                                                                                                                       <  2:476888/150‑1 (MQ=255)
                 aaaTGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAg                                                                                                                                     >  1:1832578/1‑150 (MQ=255)
                  aaTGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGa                                                                                                                                    >  2:1196835/1‑150 (MQ=255)
                  aaTGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGa                                                                                                                                    <  1:4409539/150‑1 (MQ=255)
                  aaTGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGa                                                                                                                                    >  2:1868989/1‑150 (MQ=255)
                       tctTAGAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTa                                                                                                                               <  2:2972848/150‑1 (MQ=255)
                         ttATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAAAGCGGGTTGTGTTTGCAGACGTTAAt                                                                                                                             >  2:4073732/1‑150 (MQ=255)
                           ataGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTTAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTc                                                                                                                           <  2:1341413/150‑1 (MQ=255)
                            taGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCa                                                                                                                          >  2:253120/1‑150 (MQ=255)
                             aGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATTTACTCGATGCATCACATTTTAAATAAGAGCTGGTTGTGTTTGCAGACGTTAATTCac                                                                                                                         <  2:1358254/150‑1 (MQ=255)
                              gTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCacg                                                                                                                        >  1:2830079/1‑150 (MQ=255)
                               tATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCacga                                                                                                                       >  2:1413532/1‑150 (MQ=255)
                                aTTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCacgac                                                                                                                      >  1:2684216/1‑150 (MQ=255)
                                aTTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCacgac                                                                                                                      <  2:2789777/150‑1 (MQ=255)
                                 tttGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCacgacg                                                                                                                     <  1:1012995/150‑1 (MQ=255)
                                 tttGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCacgacg                                                                                                                     <  1:19803/150‑1 (MQ=255)
                                  ttGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACgg                                                                                                                    >  2:1461046/1‑150 (MQ=255)
                                   tGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAAAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGc                                                                                                                   >  1:3607136/1‑150 (MQ=255)
                                     gCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCaa                                                                                                                 <  1:1793895/150‑1 (MQ=255)
                                          aaCGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTc                                                                                                            >  2:3536671/1‑150 (MQ=255)
                                              aTAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATaa                                                                                                        <  1:891157/150‑1 (MQ=255)
                                              aTAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATaa                                                                                                        >  2:3586728/1‑150 (MQ=255)
                                                   tATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGaaa                                                                                                   >  1:1420370/1‑150 (MQ=255)
                                                     ttttttCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAAAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACt                                                                                                 <  2:1600506/150‑1 (MQ=255)
                                                     ttttttCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAAAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACt                                                                                                 <  2:310777/150‑1 (MQ=255)
                                                           cTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGAt                                                                                           <  2:4448878/150‑1 (MQ=255)
                                                            tgatgaGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATa                                                                                          <  2:2614522/150‑1 (MQ=255)
                                                            tgatgaGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATa                                                                                          <  1:1553918/150‑1 (MQ=255)
                                                            tgatgaGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATa                                                                                          <  1:2799760/150‑1 (MQ=255)
                                                            tgatgaGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAAAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATa                                                                                          >  1:2510308/1‑150 (MQ=255)
                                                            tgatgaGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAAAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATa                                                                                          >  1:724007/1‑150 (MQ=255)
                                                             gatgaGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATaa                                                                                         >  2:3232829/1‑150 (MQ=255)
                                                              atgaGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATaaa                                                                                        >  2:1385824/1‑150 (MQ=255)
                                                               tgaGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAg                                                                                       >  1:653980/1‑150 (MQ=255)
                                                                   cGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAAt                                                                                   <  1:3857643/150‑1 (MQ=255)
                                                                    gAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATg                                                                                  >  1:4020966/1‑150 (MQ=255)
                                                                    gAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATg                                                                                  <  2:2451736/150‑1 (MQ=255)
                                                                    gAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATg                                                                                  <  1:2435935/150‑1 (MQ=255)
                                                                    gAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATg                                                                                  <  1:2435717/150‑1 (MQ=255)
                                                                     aaTGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGa                                                                                 >  2:186015/1‑150 (MQ=255)
                                                                       tGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACa                                                                               <  1:2900289/150‑1 (MQ=255)
                                                                       tGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAAAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACa                                                                               <  2:4304635/150‑1 (MQ=255)
                                                                          tGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACAtct                                                                            <  2:733246/150‑1 (MQ=255)
                                                                          tGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACAtct                                                                            <  1:1167389/150‑1 (MQ=255)
                                                                             tGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCtt                                                                         >  2:977224/1‑150 (MQ=255)
                                                                              gTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCttt                                                                        >  2:3901680/1‑150 (MQ=255)
                                                                               tCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCtttt                                                                       >  2:2210991/1‑150 (MQ=255)
                                                                                  gTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCg                                                                    >  2:1751611/1‑150 (MQ=255)
                                                                                  gTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCg                                                                    <  1:564308/150‑1 (MQ=255)
                                                                                  gTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAAAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCg                                                                    <  1:1629210/150‑1 (MQ=255)
                                                                                     tCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGttt                                                                 >  2:2957955/1‑150 (MQ=255)
                                                                                       cATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGtttct                                                               <  1:2724803/150‑1 (MQ=255)
                                                                                           cGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCa                                                           >  2:860650/1‑150 (MQ=255)
                                                                                           cGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAAAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCa                                                           >  2:4232125/1‑150 (MQ=255)
                                                                                            gAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCaa                                                          >  2:2387750/1‑150 (MQ=255)
                                                                                              ggTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCaaaa                                                        >  2:3818526/1‑150 (MQ=255)
                                                                                              ggTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCaaaa                                                        <  2:951985/150‑1 (MQ=255)
                                                                                               gTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAAt                                                       >  1:4241894/1‑150 (MQ=255)
                                                                                                tAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATc                                                      >  2:3445181/1‑150 (MQ=255)
                                                                                                 aaCAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCa                                                     >  2:527334/1‑150 (MQ=255)
                                                                                                  aCAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCAt                                                    <  1:901313/150‑1 (MQ=255)
                                                                                                  aCAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCAt                                                    <  2:3727188/150‑1 (MQ=255)
                                                                                                  aCAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAAAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCAt                                                    <  2:167718/150‑1 (MQ=255)
                                                                                                   cAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAAAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATa                                                   >  1:4479765/1‑150 (MQ=255)
                                                                                                     aaTTAAGTGACTGATCCGATATTCTCGATGCATCCAATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACTTCTCTTTTTCGTTTCTTTCAAAATCATTGa                                                 <  1:4286638/150‑1 (MQ=255)
                                                                                                        tAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAAAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAAc                                              >  1:3425821/1‑150 (MQ=255)
                                                                                                         aaGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACg                                             >  2:1576114/1‑150 (MQ=255)
                                                                                                          aGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACgg                                            >  1:2025213/1‑150 (MQ=255)
                                                                                                            tGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTa                                          >  1:1032696/1‑150 (MQ=255)
                                                                                                            tGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTa                                          >  1:3058792/1‑150 (MQ=255)
                                                                                                               cTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAgtg                                       >  2:2634500/1‑150 (MQ=255)
                                                                                                                      gATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTg                                >  1:1362245/1‑150 (MQ=255)
                                                                                                                         ataCTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGttt                             >  2:4585827/1‑150 (MQ=255)
                                                                                                                          taCTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGtttt                            <  1:3499238/150‑1 (MQ=255)
                                                                                                                          taCTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGtttt                            <  1:3499288/150‑1 (MQ=255)
                                                                                                                           aCTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGCCGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTg                           >  1:2239771/1‑150 (MQ=255)
                                                                                                                              cGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGAcc                        >  1:3718150/1‑150 (MQ=255)
                                                                                                                               gATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCt                       <  2:4404387/150‑1 (MQ=255)
                                                                                                                                  gCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGa                    >  2:1255449/1‑150 (MQ=255)
                                                                                                                                  gCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTGTGTTTTGACCTGGa                    >  2:291848/1‑150 (MQ=255)
                                                                                                                                   cATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGaa                   >  1:4077568/1‑150 (MQ=255)
                                                                                                                                   cATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGaa                   >  1:2014764/1‑150 (MQ=255)
                                                                                                                                    aTCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAAt                  <  1:3404092/150‑1 (MQ=255)
                                                                                                                                          ttttAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAatt            <  1:2123732/150‑1 (MQ=255)
                                                                                                                                             tAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAattatt         >  2:1168995/1‑150 (MQ=255)
                                                                                                                                              aaaTAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTg        <  2:3419292/150‑1 (MQ=255)
                                                                                                                                              aaaTAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTg        >  2:2812637/1‑150 (MQ=255)
                                                                                                                                              aaaTAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTg        <  2:3036728/150‑1 (MQ=255)
                                                                                                                                                aTAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTACTTTGTTTTGACCTGGAATAGAATTATTGAc      >  1:238665/1‑150 (MQ=255)
                                                                                                                                                aTAAAAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGAc      <  1:3299205/150‑1 (MQ=255)
                                                                                                                                                   agagCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGAt   >  1:4058789/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TCCCATTCGTCGTGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATG  >  CP009273/4566966‑4567263

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: