Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,567,114 G→A 10.4% 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,1140GA10.4% 1028.5 / 76.2 318S221F (TCT→TTT) mcrC5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrC
Reads supporting (aligned to +/- strand):  ref base G (144/141);  new base A (20/13);  total (164/154)
Fisher's exact test for biased strand distribution p-value = 3.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

Only 100 of 328 total aligned reads displayed.

TTCCCATTCGTCGTGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATG  >  CP009273/4566965‑4567263
                                                                                                                                                     |                                                                                                                                                     
ttCCCATTCGTCGTGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAag                                                                                                                                                       <  2:100796/150‑1 (MQ=255)
ttCCCATTCGTCGTGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAag                                                                                                                                                       <  2:2414052/150‑1 (MQ=255)
   ccATTCGTCGTGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGc                                                                                                                                                    >  2:347027/1‑150 (MQ=255)
    cATTCGTCGTGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCg                                                                                                                                                   >  2:890065/1‑150 (MQ=255)
    cATTCGTCGTGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCg                                                                                                                                                   <  1:4466793/150‑1 (MQ=255)
    cATTCGTCGTGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCg                                                                                                                                                   <  1:1389264/150‑1 (MQ=255)
     aTTCGTCGTGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCgg                                                                                                                                                  >  1:3823550/1‑150 (MQ=255)
     aTTCGTCGTGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCgg                                                                                                                                                  >  2:1215658/1‑150 (MQ=255)
     aTTCGTCGTGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCgg                                                                                                                                                  >  2:2724461/1‑150 (MQ=255)
        cgtcgtGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGtt                                                                                                                                               >  2:3389932/1‑150 (MQ=255)
         gtcgtGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTtg                                                                                                                                              >  1:3230532/1‑150 (MQ=255)
         gtcgtGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTtg                                                                                                                                              <  1:702634/150‑1 (MQ=255)
              gAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGttt                                                                                                                                         >  2:2754346/1‑150 (MQ=255)
               aaaaaaTGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTg                                                                                                                                        <  2:3233648/150‑1 (MQ=255)
               aaaaaaTGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTg                                                                                                                                        >  1:564152/1‑150 (MQ=255)
               aaaaaaTGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAAAAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTg                                                                                                                                        <  2:204203/150‑1 (MQ=255)
                 aaaaTGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAAAGCGGGTTGTGTTTGCa                                                                                                                                      <  2:1760400/150‑1 (MQ=255)
                        tctTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTa                                                                                                                               <  2:2613768/150‑1 (MQ=255)
                        tctTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTa                                                                                                                               <  2:296975/150‑1 (MQ=255)
                                   ttGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACgg                                                                                                                    >  2:3160578/1‑150 (MQ=255)
                                    tGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGc                                                                                                                   >  2:692825/1‑150 (MQ=255)
                                    tGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAAAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGc                                                                                                                   >  1:297175/1‑150 (MQ=255)
                                        gTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCaaaa                                                                                                               <  1:1942657/150‑1 (MQ=255)
                                        gTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCaaaa                                                                                                               <  1:4272957/150‑1 (MQ=255)
                                         tCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAAt                                                                                                              <  1:1881129/150‑1 (MQ=255)
                                         tCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAAAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAAt                                                                                                              <  1:2664091/150‑1 (MQ=255)
                                           aaCGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTc                                                                                                            >  1:1314899/1‑150 (MQ=255)
                                           aaCGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTc                                                                                                            <  2:3044323/150‑1 (MQ=255)
                                            aCGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCa                                                                                                           >  1:2011832/1‑150 (MQ=255)
                                             cGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCAt                                                                                                          >  1:1034280/1‑150 (MQ=255)
                                             cGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCAt                                                                                                          >  1:2820177/1‑150 (MQ=255)
                                             cGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAAAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCAt                                                                                                          <  2:2225432/150‑1 (MQ=255)
                                               aTAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATaa                                                                                                        <  1:3122348/150‑1 (MQ=255)
                                                tAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATaaa                                                                                                       <  2:1327962/150‑1 (MQ=255)
                                                 aaGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAg                                                                                                      <  2:1638517/150‑1 (MQ=255)
                                                 aaGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAg                                                                                                      <  1:2572714/150‑1 (MQ=255)
                                                 aaGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAg                                                                                                      <  2:2848587/150‑1 (MQ=255)
                                                     aTTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAAc                                                                                                  >  1:4046681/1‑150 (MQ=255)
                                                       tttttCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACtt                                                                                                <  2:2491459/150‑1 (MQ=255)
                                                       tttttCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACtt                                                                                                <  1:4184355/150‑1 (MQ=255)
                                                         tttCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACtttt                                                                                              <  1:2726467/150‑1 (MQ=255)
                                                            cTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGAt                                                                                           <  2:3381432/150‑1 (MQ=255)
                                                             tgatgaGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATa                                                                                          <  2:2149561/150‑1 (MQ=255)
                                                             tgatgaGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATa                                                                                          <  1:2256467/150‑1 (MQ=255)
                                                              gatgaGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATaa                                                                                         <  1:2853861/150‑1 (MQ=255)
                                                               atgaGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATaaa                                                                                        >  2:4456070/1‑150 (MQ=255)
                                                                tgaGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAAAGCGGGTTGTGTATGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAg                                                                                       >  1:3537407/1‑150 (MQ=255)
                                                                    cGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAAt                                                                                   <  1:1174933/150‑1 (MQ=255)
                                                                     gAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGATACTTTTGATAAAGTAATg                                                                                  <  2:616360/150‑1 (MQ=255)
                                                                     gAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATg                                                                                  <  2:4181723/150‑1 (MQ=255)
                                                                     gAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATg                                                                                  >  1:1724619/1‑150 (MQ=255)
                                                                       aTGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGAc                                                                                <  1:1020272/150‑1 (MQ=255)
                                                                             aTGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATctct                                                                          <  1:1767880/150‑1 (MQ=255)
                                                                                     ttCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGtt                                                                  >  1:3450223/1‑150 (MQ=255)
                                                                                        cATTCGAGGTAACAAATTCAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGtttct                                                               <  2:1907657/150‑1 (MQ=255)
                                                                                         aTTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGtttctt                                                              >  1:783728/1‑150 (MQ=255)
                                                                                          ttCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGtttcttt                                                             <  1:3573023/150‑1 (MQ=255)
                                                                                          ttCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGtttcttt                                                             <  1:3319817/150‑1 (MQ=255)
                                                                                          ttCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGtttcttt                                                             <  1:1642572/150‑1 (MQ=255)
                                                                                             gAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCaa                                                          >  1:3123889/1‑150 (MQ=255)
                                                                                              aGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCaaa                                                         >  1:2482840/1‑150 (MQ=255)
                                                                                              aGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCaaa                                                         <  1:2228764/150‑1 (MQ=255)
                                                                                                gTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAAt                                                       >  2:2386768/1‑150 (MQ=255)
                                                                                                  aaCAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCa                                                     <  2:3899735/150‑1 (MQ=255)
                                                                                                  aaCAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCa                                                     >  2:3690121/1‑150 (MQ=255)
                                                                                                     aaaTTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAg                                                  >  2:56258/1‑150 (MQ=255)
                                                                                                         tAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAAc                                              <  1:2102177/150‑1 (MQ=255)
                                                                                                          aaGTGACTGATCCGATATACTCGATGCATCCGATTTTAAATAAAAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACg                                             <  2:2951558/150‑1 (MQ=255)
                                                                                                          aaGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAAAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACg                                             >  1:2343594/1‑150 (MQ=255)
                                                                                                             tGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTa                                          >  1:336815/1‑150 (MQ=255)
                                                                                                               aCTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAgt                                        >  2:2276796/1‑150 (MQ=255)
                                                                                                               aCTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAgt                                        >  2:2693785/1‑150 (MQ=255)
                                                                                                               aCTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAgt                                        >  2:1457800/1‑150 (MQ=255)
                                                                                                                 tGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAgtgt                                      >  2:38708/1‑150 (MQ=255)
                                                                                                                   aTCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTcc                                    >  2:2293731/1‑150 (MQ=255)
                                                                                                                     ccGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCtt                                  >  1:3902661/1‑150 (MQ=255)
                                                                                                                      cGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCttt                                 <  2:3638534/150‑1 (MQ=255)
                                                                                                                      cGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCttt                                 >  2:1863284/1‑150 (MQ=255)
                                                                                                                       gATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTg                                >  2:4590549/1‑150 (MQ=255)
                                                                                                                       gATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAAAGGTAGTGTCCTTTg                                <  2:3122939/150‑1 (MQ=255)
                                                                                                                        atatACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGt                               >  1:4634721/1‑150 (MQ=255)
                                                                                                                           taCTCGATGCATCCCATTTTAAATAAAAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGtttt                            >  1:4616415/1‑150 (MQ=255)
                                                                                                                            aCTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTg                           >  1:140854/1‑150 (MQ=255)
                                                                                                                              tCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGAc                         >  2:1805901/1‑150 (MQ=255)
                                                                                                                              tCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGAc                         >  1:1395046/1‑150 (MQ=255)
                                                                                                                              tCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGAc                         <  1:4342261/150‑1 (MQ=255)
                                                                                                                                   gCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGa                    <  1:3988707/150‑1 (MQ=255)
                                                                                                                                   gCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGa                    >  1:3941179/1‑150 (MQ=255)
                                                                                                                                     aTCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAAt                  <  1:647439/150‑1 (MQ=255)
                                                                                                                                      tCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATa                 <  2:283055/150‑1 (MQ=255)
                                                                                                                                        ccATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTATTTTTCGTTTCTTTCAAAATCATATAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGa               >  1:1142600/1‑150 (MQ=255)
                                                                                                                                               aaaTAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTg        <  1:1614789/150‑1 (MQ=255)
                                                                                                                                               aaaTAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTg        >  2:635622/1‑150 (MQ=255)
                                                                                                                                                aaTAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGa       >  2:1386157/1‑150 (MQ=255)
                                                                                                                                                 aTAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGAc      <  1:4564672/150‑1 (MQ=255)
                                                                                                                                                  tAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACg     <  1:1679178/150‑1 (MQ=255)
                                                                                                                                                   aaGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGa    >  1:271259/1‑150 (MQ=255)
                                                                                                                                                   aaGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGa    >  2:4217120/1‑150 (MQ=255)
                                                                                                                                                   aaGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGa    >  1:2126834/1‑150 (MQ=255)
                                                                                                                                                     gagCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATg  <  2:4041742/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTCCCATTCGTCGTGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATG  >  CP009273/4566965‑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: