Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,567,278 G→T 100% I166I (ATC→ATA 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,2780GT99.2% 776.0 / ‑1.7 237I166I (ATC→ATAmcrC5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrC
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (128/107);  total (128/109)
Fisher's exact test for biased strand distribution p-value = 2.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 239 total aligned reads displayed.

GACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGCT  >  CP009273/4567132‑4567423
                                                                                                                                                  |                                                                                                                                                 
 aCGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTtata                                                                                                                                                <  2:1009804/149‑1 (MQ=255)
  cGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATaa                                                                                                                                               <  2:139522/149‑1 (MQ=255)
    ttAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACg                                                                                                                                             >  2:2570232/1‑149 (MQ=255)
     tAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGa                                                                                                                                            <  1:438516/149‑1 (MQ=255)
     tAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGa                                                                                                                                            <  1:2407228/149‑1 (MQ=255)
      aaTTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGaa                                                                                                                                           >  1:2109164/1‑149 (MQ=255)
          cACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTta                                                                                                                                       <  1:604255/149‑1 (MQ=255)
            cgacgGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTtata                                                                                                                                     >  2:682168/1‑149 (MQ=255)
             gacgGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTataa                                                                                                                                    >  1:576796/1‑149 (MQ=255)
             gacgGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTataa                                                                                                                                    >  1:443598/1‑149 (MQ=255)
             gacgGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTataa                                                                                                                                    >  1:2455052/1‑149 (MQ=255)
              acgGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTataat                                                                                                                                   <  2:4108302/149‑1 (MQ=255)
               cgGCAAAATTCATTAATAAACTTTTGATAAAGTAATGACTTCTCTTTTTCGTTTCTTTCAAATTCATAGAAACGGTAGTGTCCTGTGTTTTGACCTGGAATAGAATTATTTACGATGAATTTGCATACAATTATAACGAATTTataatt                                                                                                                                  <  2:4062942/149‑2 (MQ=255)
               cgGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTataata                                                                                                                                  <  1:3556513/149‑1 (MQ=255)
               atgaaataTTCATAAAGAAACTTTTGATAAAGTATTGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACTGTAGTGTCCTTTGTTTTGACCTGGATTAGAATTATTGACGATTAATTTGCAGAAACTTATAACGAATTTataata                                                                                                                                  <  1:3955164/142‑1 (MQ=255)
                ggCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTAtaataa                                                                                                                                 <  1:3788390/149‑1 (MQ=255)
                ggCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTAtaataa                                                                                                                                 <  2:3913459/149‑1 (MQ=255)
                ggCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTAtaataa                                                                                                                                 <  2:3913833/149‑1 (MQ=255)
                 gCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAAc                                                                                                                                <  1:940138/149‑1 (MQ=255)
                   aaaaTTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAAcgc                                                                                                                              <  1:827683/149‑1 (MQ=255)
                       ttaaTAAAGAAACTTTTGATAAATTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTAt                                                                                                                          <  2:3469643/146‑1 (MQ=255)
                       ttCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTAt                                                                                                                          >  2:3611566/1‑149 (MQ=255)
                          aTAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTAtttt                                                                                                                       >  1:3001774/1‑149 (MQ=255)
                          aTAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTAtttt                                                                                                                       <  2:78902/149‑1 (MQ=255)
                          aTAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTAtttt                                                                                                                       >  1:718670/1‑149 (MQ=255)
                           tAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTAttttt                                                                                                                      >  1:3265565/1‑149 (MQ=255)
                             aaGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGTAATAGAATTATTGACGATGAATTTGCAGAAACTTATAACGAATTTATAATAACGCGTATTTTTTc                                                                                                                    >  2:832132/1‑149 (MQ=255)
                                   cTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCAt                                                                                                              >  1:2285201/1‑149 (MQ=255)
                                    ttttGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTCTAACGAATTTATAATAACGCGTATTTTTTCCGCCAtt                                                                                                             <  1:2080912/149‑1 (MQ=255)
                                    ttttGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCAtt                                                                                                             >  2:2163218/1‑149 (MQ=255)
                                         aTAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGAt                                                                                                        >  2:3876741/1‑149 (MQ=255)
                                         aTAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGAt                                                                                                        >  1:2299277/1‑149 (MQ=255)
                                          tAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGTCCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATa                                                                                                       >  2:2821518/1‑149 (MQ=255)
                                               tAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATATAATTATTTACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGa                                                                                                  >  1:1393363/1‑149 (MQ=255)
                                               tAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGa                                                                                                  >  1:2796386/1‑149 (MQ=255)
                                                   gACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAAtg                                                                                              >  1:2896485/1‑149 (MQ=255)
                                                    aCATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAAtgc                                                                                             >  1:1925254/1‑149 (MQ=255)
                                                    aCATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAAtgc                                                                                             <  1:2202745/149‑1 (MQ=255)
                                                         tctTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCgg                                                                                        >  1:1302959/1‑149 (MQ=255)
                                                             tttCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGtt                                                                                    <  2:310230/149‑1 (MQ=255)
                                                                 gTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAAt                                                                                >  1:2731020/1‑149 (MQ=255)
                                                                  tttctttcAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATg                                                                               >  2:3359027/1‑149 (MQ=255)
                                                                   ttctttcAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGa                                                                              >  2:235191/1‑149 (MQ=255)
                                                                       ttcAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCTCCATTCAGATATCTGAAATGCTGCGGAGTTAAATGAagag                                                                          >  2:110188/1‑149 (MQ=255)
                                                                       ttcAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATATCTGAAATGCTGCGGAGTTAAATGAagag                                                                          >  2:2143462/1‑149 (MQ=255)
                                                                       ttcAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAagag                                                                          <  1:1125225/149‑1 (MQ=255)
                                                                        tcAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGt                                                                         <  2:2413898/149‑1 (MQ=255)
                                                                           aaaTCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCt                                                                      >  2:2313128/1‑149 (MQ=255)
                                                                             aTCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTaa                                                                    >  2:1575871/1‑149 (MQ=255)
                                                                              tCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAAt                                                                   <  1:482539/149‑1 (MQ=255)
                                                                               cATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATg                                                                  >  2:808636/1‑149 (MQ=255)
                                                                                 tAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGcc                                                                >  1:3812642/1‑149 (MQ=255)
                                                                                    aaaCGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCgg                                                             >  1:3920168/1‑149 (MQ=255)
                                                                                     aaCGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGt                                                            >  1:192082/1‑149 (MQ=255)
                                                                                      aCGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTa                                                           >  1:261585/1‑149 (MQ=255)
                                                                                         gTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAAtt                                                        <  2:213585/149‑1 (MQ=255)
                                                                                         gTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAAtt                                                        <  2:2290986/149‑1 (MQ=255)
                                                                                         gTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAAtt                                                        <  2:3216359/149‑1 (MQ=255)
                                                                                            gtgtCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTc                                                     <  2:1759852/149‑1 (MQ=255)
                                                                                            gtgtCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTc                                                     >  2:3441311/1‑149 (MQ=255)
                                                                                             tgtCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCt                                                    >  2:2275360/1‑149 (MQ=255)
                                                                                              gtCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCta                                                   >  1:3936345/1‑149 (MQ=255)
                                                                                                  tttGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATaaa                                               >  1:878383/1‑149 (MQ=255)
                                                                                                  tttGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATaaa                                               >  1:878241/1‑149 (MQ=255)
                                                                                                   ttGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAg                                              >  1:1614211/1‑149 (MQ=255)
                                                                                                   ttGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAg                                              >  1:770086/1‑149 (MQ=255)
                                                                                                   ttGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAg                                              >  1:496064/1‑149 (MQ=255)
                                                                                                   ttGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAg                                              <  1:3717539/149‑1 (MQ=255)
                                                                                                      ttttGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGa                                           >  2:3656487/1‑149 (MQ=255)
                                                                                                           aCCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAg                                      <  2:1839399/149‑1 (MQ=255)
                                                                                                                  aTAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCAtc                               <  2:2896485/149‑1 (MQ=255)
                                                                                                                   tAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGAGGAGTTAAATGAAGAGTGCTAATGCCAGGTACTTTTCTATAAAGTGAACGAGCTACAtct                              >  2:3800307/1‑149 (MQ=255)
                                                                                                                     gAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATctct                            >  2:1526345/1‑149 (MQ=255)
                                                                                                                     gAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATctct                            <  2:718670/149‑1 (MQ=255)
                                                                                                                     gAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATctct                            >  1:1097438/1‑149 (MQ=255)
                                                                                                                     gAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATctct                            <  2:2541269/149‑1 (MQ=255)
                                                                                                                     gAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATctct                            >  1:3027375/1‑149 (MQ=255)
                                                                                                                        ttattGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGAt                         >  2:950790/1‑149 (MQ=255)
                                                                                                                         tattGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATa                        <  1:3423538/149‑1 (MQ=255)
                                                                                                                           ttGACGATGAATTTGCAGACACTTATAACGAATTTATACTAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGt                      >  2:3504117/1‑149 (MQ=255)
                                                                                                                           ttGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGt                      <  1:368996/149‑1 (MQ=255)
                                                                                                                              aCGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGa                   >  1:420804/1‑149 (MQ=255)
                                                                                                                              aCGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGa                   <  1:1794866/149‑1 (MQ=255)
                                                                                                                               cGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGaa                  >  1:3331591/1‑149 (MQ=255)
                                                                                                                               cGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGaa                  >  1:1280266/1‑149 (MQ=255)
                                                                                                                                 aTGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAAtt                >  2:1710041/1‑149 (MQ=255)
                                                                                                                                 aTGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAAtt                >  2:888474/1‑149 (MQ=255)
                                                                                                                                 aTGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAAtt                >  2:2946156/1‑149 (MQ=255)
                                                                                                                                 aTGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAAtt                <  2:3738793/149‑1 (MQ=255)
                                                                                                                                  tgtaTTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAAttt               <  2:3343760/146‑1 (MQ=255)
                                                                                                                                   gAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTa              >  1:3158718/1‑149 (MQ=255)
                                                                                                                                   gAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTa              <  1:1487028/149‑1 (MQ=255)
                                                                                                                                    aaTTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTaa             >  1:3240138/1‑149 (MQ=255)
                                                                                                                                      tttGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACt           <  1:3876741/149‑1 (MQ=255)
                                                                                                                                        tGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACttt         <  2:1159448/149‑1 (MQ=255)
                                                                                                                                         gCAGACACTTATAACGAATTTATAATAACGCGTATTTTTACCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACtttt        <  2:1038892/149‑1 (MQ=255)
                                                                                                                                           aGACACTTATAACGAATTTATAATTACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCa      <  2:3073692/149‑1 (MQ=255)
                                                                                                                                           aGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCa      <  1:2570232/149‑1 (MQ=255)
                                                                                                                                              cacTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAATAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGc   <  2:443598/149‑1 (MQ=255)
                                                                                                                                               acTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGCt  >  1:1642252/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGCT  >  CP009273/4567132‑4567423

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