Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,572,800 A→T 42.1% V22V (GTT→GTA)  hsdM ← DNA methyltransferase M

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,572,8000AT42.1% 128.6 / 311.0 221V22V (GTT→GTA) hsdMDNA methyltransferase M
Reads supporting (aligned to +/- strand):  ref base A (62/65);  new base T (54/39);  total (116/105)
Fisher's exact test for biased strand distribution p-value = 2.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

Only 100 of 229 total aligned reads displayed.

GGTAGAACTGCAACTGCTCCTGGCCGATGCGGGATTTCAGGTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCAT  >  CP009273/4572652‑4572949
                                                                                                                                                    |                                                                                                                                                     
ggTAGAACTGCAACTGCTCCTGGCCGATGCGGGATTTCAGGTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGaaa                                                                                                                                                      >  2:2850702/1‑150 (MQ=255)
 gTAGAACTGCAACTGCTCCTGGCCGATGCGGGATTTCAGGTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATAc                                                                                                                                                     <  2:1637059/150‑1 (MQ=255)
 gTAGAACTGCAACTGCTCCTGGCCGATGCGGGATTTCAGGTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATAc                                                                                                                                                     <  2:3228112/150‑1 (MQ=255)
 gTAGAACTGCAACTGCTCCTGGCCGATGCGGGATTTCAGGTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATAc                                                                                                                                                     <  1:2185643/150‑1 (MQ=255)
   aGAACTGCAACTGCTCCTGGCCGATGCGGGATTTCAGGTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATAcgc                                                                                                                                                   <  2:2414915/150‑1 (MQ=255)
    gAACTGCAACTGCTCCTGGCCGATGCGGGATTTCAGGTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAAcgcc                                                                                                                                                  <  1:2748659/150‑1 (MQ=255)
         gcaactgcTCCTGGCCGATGCGGGATTTCAGGTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCAt                                                                                                                                             <  2:1628652/150‑1 (MQ=255)
                cTCCTGGCCGATGCGGGATTTCAGGTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAg                                                                                                                                      >  2:3655603/1‑150 (MQ=255)
                cTCCTGGCCGATGCGGGATTTCAGGTCACCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAg                                                                                                                                      <  1:1501514/150‑1 (MQ=255)
                 tCCTGGCCGATGCGGGATTTCAGGTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAgg                                                                                                                                     <  1:1286314/150‑1 (MQ=255)
                 tCCTGGCCGATGCGGGATTTCAGGTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAgg                                                                                                                                     >  1:581525/1‑150 (MQ=255)
                        cGATGCGGGATTTCAGGTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGc                                                                                                                              >  2:1045708/1‑150 (MQ=255)
                         gATGCGGGATTTCAGGTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGca                                                                                                                             >  2:2066113/1‑150 (MQ=255)
                          aTGCGGGATTTCAGGTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGcac                                                                                                                            <  2:2721958/150‑1 (MQ=255)
                                gATTTCAGGTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGcacagcttc                                                                                                                      >  2:3250463/1‑150 (MQ=255)
                                 aTTTCAGGTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGcacagcttcc                                                                                                                     <  1:248303/150‑1 (MQ=255)
                                   ttCAGGTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGcacagcttccac                                                                                                                   <  1:611996/150‑1 (MQ=255)
                                        gTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGcacagcttccacagctt                                                                                                              <  2:2547712/150‑1 (MQ=255)
                                         tCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGcacagcttccacagcttc                                                                                                             <  2:3993290/150‑1 (MQ=255)
                                         tCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGcacagcttccacagcttc                                                                                                             <  2:3537786/150‑1 (MQ=255)
                                          cATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTcg                                                                                                            <  1:3726264/150‑1 (MQ=255)
                                          cATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTcg                                                                                                            <  2:3386506/150‑1 (MQ=255)
                                            tCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTcgcg                                                                                                          <  1:895413/150‑1 (MQ=255)
                                             cccAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGa                                                                                                         >  1:3506964/1‑150 (MQ=255)
                                             cccAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGa                                                                                                         >  2:3324147/1‑150 (MQ=255)
                                              ccAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGAc                                                                                                        <  1:1327403/150‑1 (MQ=255)
                                              ccAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGAc                                                                                                        <  2:1385668/150‑1 (MQ=255)
                                               cAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGAcc                                                                                                       >  2:793548/1‑150 (MQ=255)
                                                  cGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGa                                                                                                    >  2:1031461/1‑150 (MQ=255)
                                                   ggTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGAt                                                                                                   >  2:3873743/1‑150 (MQ=255)
                                                   ggTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGAt                                                                                                   >  2:2254702/1‑150 (MQ=255)
                                                    gTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATc                                                                                                  <  1:785348/150‑1 (MQ=255)
                                                     tAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCg                                                                                                 >  1:3335565/1‑150 (MQ=255)
                                                     tAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCg                                                                                                 >  2:1823174/1‑150 (MQ=255)
                                                          ttCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTAtt                                                                                            <  2:1520896/150‑1 (MQ=255)
                                                            ccGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGt                                                                                          >  1:2069883/1‑150 (MQ=255)
                                                             cGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGtt                                                                                         >  2:293043/1‑150 (MQ=255)
                                                             cGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGtt                                                                                         >  2:2782238/1‑150 (MQ=255)
                                                              ggCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTc                                                                                        <  1:2013810/150‑1 (MQ=255)
                                                              ggCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTc                                                                                        >  2:3624191/1‑150 (MQ=255)
                                                               gCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCa                                                                                       <  1:3395825/150‑1 (MQ=255)
                                                                cAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCAt                                                                                      >  1:2734091/1‑150 (MQ=255)
                                                                 aGGTATTCCGTTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCAtt                                                                                     >  2:1090549/1‑150 (MQ=255)
                                                                 aGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCAtt                                                                                     <  1:2220419/150‑1 (MQ=255)
                                                                  ggTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATtg                                                                                    >  2:3797397/1‑150 (MQ=255)
                                                                         cGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTc                                                                             <  1:1704244/150‑1 (MQ=255)
                                                                           cTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCg                                                                           >  1:4066556/1‑150 (MQ=255)
                                                                            ttCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGt                                                                          <  1:3879785/150‑1 (MQ=255)
                                                                                tGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAAt                                                                      <  1:2051509/150‑1 (MQ=255)
                                                                                  aCCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTa                                                                    <  1:417859/150‑1 (MQ=255)
                                                                                    cGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAag                                                                  >  1:2352829/1‑150 (MQ=255)
                                                                                     ggTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAagc                                                                 <  2:848847/150‑1 (MQ=255)
                                                                                     ggTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAagc                                                                 >  1:1771661/1‑150 (MQ=255)
                                                                                       tctcTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAagcag                                                               <  2:86226/150‑1 (MQ=255)
                                                                                       tctcTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAagcag                                                               >  2:2448859/1‑150 (MQ=255)
                                                                                        ctctTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAagcagc                                                              <  2:1844370/150‑1 (MQ=255)
                                                                                        ctctTTAAACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGGACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAagcagc                                                              <  2:725664/150‑1 (MQ=255)
                                                                                          ctTTCCACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCgg                                                            <  1:3011729/150‑1 (MQ=255)
                                                                                          ctTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCgg                                                            >  2:510511/1‑150 (MQ=255)
                                                                                                 caTTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAAt                                                     >  2:3685438/1‑150 (MQ=255)
                                                                                                 caTTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAAt                                                     >  2:496916/1‑150 (MQ=255)
                                                                                                  aTTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAAtt                                                    >  2:1137895/1‑150 (MQ=255)
                                                                                                    tttCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCa                                                  >  2:3304573/1‑150 (MQ=255)
                                                                                                    tttCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCa                                                  >  2:2459671/1‑150 (MQ=255)
                                                                                                      tCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATc                                                <  1:2539799/150‑1 (MQ=255)
                                                                                                       aaaaaaCAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCg                                               <  2:1025535/149‑1 (MQ=255)
                                                                                                        aaaaaCAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGa                                              >  1:3653469/1‑150 (MQ=255)
                                                                                                        aaaaaCAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGa                                              >  2:2164402/1‑150 (MQ=255)
                                                                                                        aaaaaCAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGa                                              >  1:1894961/1‑150 (MQ=255)
                                                                                                           aaCAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGcc                                           >  2:273766/1‑150 (MQ=255)
                                                                                                           aaCAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGcc                                           <  1:1693071/150‑1 (MQ=255)
                                                                                                             cagcagCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGc                                         >  1:3186505/1‑150 (MQ=255)
                                                                                                              agcagcGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCa                                        <  2:426682/150‑1 (MQ=255)
                                                                                                                 agcGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCagaa                                     <  2:4071412/150‑1 (MQ=255)
                                                                                                                 agcGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCagaa                                     >  2:621476/1‑150 (MQ=255)
                                                                                                                  gcgaggcgagTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCagaag                                    <  2:4164452/150‑1 (MQ=255)
                                                                                                                     aggcgagTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCagaagaag                                 >  2:2764934/1‑150 (MQ=255)
                                                                                                                      ggcgagTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACATCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCagaagatga                                <  1:438397/150‑1 (MQ=255)
                                                                                                                       gcgagTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAgaagaagaa                               >  2:1456501/1‑150 (MQ=255)
                                                                                                                          agTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAAtt                            <  2:3157090/150‑1 (MQ=255)
                                                                                                                          agTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAAtt                            <  1:2604763/150‑1 (MQ=255)
                                                                                                                             tCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGcc                         <  1:328319/150‑1 (MQ=255)
                                                                                                                             tCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGcc                         >  2:1810045/1‑150 (MQ=255)
                                                                                                                                  gACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAggg                    >  1:1007032/1‑150 (MQ=255)
                                                                                                                                   aCGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGt                   <  1:1470795/150‑1 (MQ=255)
                                                                                                                                    cGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTa                  <  1:3413427/150‑1 (MQ=255)
                                                                                                                                       aGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAata               >  1:3075423/1‑150 (MQ=255)
                                                                                                                                          tttGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATac            <  2:3501652/150‑1 (MQ=255)
                                                                                                                                           ttGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATaca           >  1:2895885/1‑150 (MQ=255)
                                                                                                                                             gATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATAcaca         >  1:1402903/1‑150 (MQ=255)
                                                                                                                                              aTAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACaa        >  2:3007876/1‑150 (MQ=255)
                                                                                                                                               tAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACaaa       >  1:1595287/1‑150 (MQ=255)
                                                                                                                                               tAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACaaa       >  1:2954734/1‑150 (MQ=255)
                                                                                                                                                aGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACaaaa      >  2:2240167/1‑150 (MQ=255)
                                                                                                                                                aGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACaaaa      >  2:1943428/1‑150 (MQ=255)
                                                                                                                                                  gAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATc    <  2:4128454/150‑1 (MQ=255)
                                                                                                                                                  gAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATc    <  2:510184/150‑1 (MQ=255)
                                                                                                                                                  gAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATc    >  1:2674861/1‑150 (MQ=255)
                                                                                                                                                   aaaCGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCa   <  1:674762/150‑1 (MQ=255)
                                                                                                                                                    aaCGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCAt  <  2:2821135/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGTAGAACTGCAACTGCTCCTGGCCGATGCGGGATTTCAGGTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCAT  >  CP009273/4572652‑4572949

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: ‑