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,2780GT96.7% 664.4 / 11.8 215I166I (ATC→ATA) mcrC5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrC
Reads supporting (aligned to +/- strand):  ref base G (4/3);  new base T (115/93);  total (119/96)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 217 total aligned reads displayed.

ACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGCTT  >  CP009273/4567133‑4567424
                                                                                                                                                 |                                                                                                                                                  
aCGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTAt                                                                                                                                                  >  1:1899684/1‑148 (MQ=255)
 cGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTtata                                                                                                                                                 >  1:1400901/1‑148 (MQ=255)
  gTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATaa                                                                                                                                                <  2:2408891/148‑1 (MQ=255)
    tAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACg                                                                                                                                              >  1:3222064/1‑148 (MQ=255)
       ttCACGACGGCAAAATTCATAAGGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAAt                                                                                                                                           <  2:2035236/148‑1 (MQ=255)
       ttCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAAt                                                                                                                                           >  1:2951516/1‑148 (MQ=255)
        tCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAAtt                                                                                                                                          >  1:1685005/1‑148 (MQ=255)
           cgacgGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTtat                                                                                                                                       <  2:783683/148‑1 (MQ=255)
            gacgGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTtata                                                                                                                                      <  2:1185379/148‑1 (MQ=255)
             acgGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTataa                                                                                                                                     >  2:827955/1‑148 (MQ=255)
              cgGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTataat                                                                                                                                    >  1:1311865/1‑148 (MQ=255)
                 cAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAAc                                                                                                                                 >  1:2195584/1‑148 (MQ=255)
                   aaaTTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAAcgc                                                                                                                               >  2:938720/1‑148 (MQ=255)
                     aTTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGt                                                                                                                             >  1:829931/1‑148 (MQ=255)
                     aTTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGt                                                                                                                             <  1:2581996/148‑1 (MQ=255)
                      ttCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTa                                                                                                                            >  2:918122/1‑148 (MQ=255)
                      ttCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTa                                                                                                                            >  2:1394675/1‑148 (MQ=255)
                      ttCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTa                                                                                                                            >  2:1570565/1‑148 (MQ=255)
                      ttCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTa                                                                                                                            >  2:1183939/1‑148 (MQ=255)
                      ttCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTa                                                                                                                            >  1:2606459/1‑148 (MQ=255)
                          tAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTAtttt                                                                                                                        <  1:42922/148‑1 (MQ=255)
                           aaaGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTAttttt                                                                                                                       <  2:1246404/148‑1 (MQ=255)
                            aaGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTAtttttt                                                                                                                      <  2:2495564/148‑1 (MQ=255)
                            aaGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTAtttttt                                                                                                                      >  2:531735/1‑148 (MQ=255)
                             aGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTc                                                                                                                     <  1:755399/148‑1 (MQ=255)
                              gAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTcc                                                                                                                    >  1:2304626/1‑148 (MQ=255)
                               aaaCTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCGTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCg                                                                                                                   >  2:1989325/1‑148 (MQ=255)
                                    tttGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCAtt                                                                                                              <  2:2116565/148‑1 (MQ=255)
                                    tttGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCAtt                                                                                                              >  2:316146/1‑148 (MQ=255)
                                    tttGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCAtt                                                                                                              <  1:1003470/148‑1 (MQ=255)
                                    tttGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCAtt                                                                                                              >  2:2619471/1‑148 (MQ=255)
                                    tttGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAAATTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCAtt                                                                                                              >  1:2990277/1‑148 (MQ=255)
                                     ttGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTc                                                                                                             <  1:2262892/148‑1 (MQ=255)
                                     ttGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCAGCCATTc                                                                                                             >  2:2726672/1‑148 (MQ=255)
                                          aaaGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATa                                                                                                        >  1:3050192/1‑148 (MQ=255)
                                           aaGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAg                                                                                                       <  1:1570565/148‑1 (MQ=255)
                                            aGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGc                                                                                                      >  1:3261506/1‑148 (MQ=255)
                                               aaTGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATATTTTCCGCCATTCAGATAGCTGa                                                                                                   >  2:1523563/1‑148 (MQ=255)
                                                       ctctTTTTCGTTTCTTTCAAAATCATCGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAAtgctgc                                                                                           <  2:2343729/148‑1 (MQ=255)
                                                         ctTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCgg                                                                                         >  2:2805456/1‑148 (MQ=255)
                                                         ctTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCgg                                                                                         >  1:2515649/1‑148 (MQ=255)
                                                           ttttCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAg                                                                                       >  2:1679438/1‑148 (MQ=255)
                                                              tCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTa                                                                                    >  2:2611067/1‑148 (MQ=255)
                                                               cGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTaa                                                                                   >  2:461429/1‑148 (MQ=255)
                                                               cGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTaa                                                                                   <  2:1963292/148‑1 (MQ=255)
                                                                gTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTaaa                                                                                  <  2:2951516/148‑1 (MQ=255)
                                                                 tttctttcAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGGTAGCTGAAATGCTGCGGAGTTAAAt                                                                                 >  1:365945/1‑148 (MQ=255)
                                                                  ttctttcAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATg                                                                                >  2:3287115/1‑148 (MQ=255)
                                                                      ttcAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAaga                                                                            >  2:1698630/1‑148 (MQ=255)
                                                                      ttcAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAaga                                                                            <  1:2549506/148‑1 (MQ=255)
                                                                       tcAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAagag                                                                           <  1:889233/148‑1 (MQ=255)
                                                                        cAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGt                                                                          >  1:558295/1‑148 (MQ=255)
                                                                          aaaTCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGc                                                                        >  1:2393875/1‑148 (MQ=255)
                                                                            aTCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTa                                                                      >  1:2902331/1‑148 (MQ=255)
                                                                             tCATAGAAACGGTAGTGTCCTTTGTTTTGACCTTGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTaa                                                                     <  1:3121153/148‑1 (MQ=255)
                                                                             tCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTaa                                                                     >  2:588084/1‑148 (MQ=255)
                                                                               aTAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATg                                                                   <  2:942611/148‑1 (MQ=255)
                                                                                   aaaCGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCg                                                               <  1:3282119/148‑1 (MQ=255)
                                                                                    aaCGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCgg                                                              <  1:276058/148‑1 (MQ=255)
                                                                                      cGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTa                                                            >  1:1548599/1‑148 (MQ=255)
                                                                                      cGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTa                                                            >  1:2661955/1‑148 (MQ=255)
                                                                                      cGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTa                                                            <  1:1184268/148‑1 (MQ=255)
                                                                                         tAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAAtt                                                         >  1:2905373/1‑148 (MQ=255)
                                                                                            tgtCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAAAGTGCTAATGCCCGGTAATTTTc                                                      >  2:2443699/1‑148 (MQ=255)
                                                                                                 tttGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATaa                                                 >  1:2560885/1‑148 (MQ=255)
                                                                                                      tttGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCAGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGa                                            <  2:1899684/148‑1 (MQ=255)
                                                                                                         gACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCTCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACg                                         <  2:3050290/148‑1 (MQ=255)
                                                                                                          aCCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGa                                        <  2:3384470/148‑1 (MQ=255)
                                                                                                               gAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTc                                   >  2:2529392/1‑148 (MQ=255)
                                                                                                               gAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTc                                   <  2:389130/148‑1 (MQ=255)
                                                                                                                aaTAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCa                                  <  2:1248494/148‑1 (MQ=255)
                                                                                                                aaTAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCa                                  >  1:3326721/1‑148 (MQ=255)
                                                                                                                   aGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCAtct                               <  1:813326/148‑1 (MQ=255)
                                                                                                                   aGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCAtct                               <  2:258225/148‑1 (MQ=255)
                                                                                                                   aGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAAAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCAtct                               <  1:1565041/148‑1 (MQ=255)
                                                                                                                        tattGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGAt                          <  1:3283762/148‑1 (MQ=255)
                                                                                                                         attGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATa                         <  1:1146590/148‑1 (MQ=255)
                                                                                                                          ttGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAAAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAg                        >  2:1466108/1‑148 (MQ=255)
                                                                                                                           tGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGt                       <  2:1426667/148‑1 (MQ=255)
                                                                                                                           tGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGt                       >  2:1267443/1‑148 (MQ=255)
                                                                                                                           tGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGt                       >  1:992397/1‑148 (MQ=255)
                                                                                                                           tGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGt                       <  1:937771/148‑1 (MQ=255)
                                                                                                                             aCGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTg                     <  2:1589620/148‑1 (MQ=255)
                                                                                                                              ctaTGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGa                    <  1:2993118/146‑1 (MQ=255)
                                                                                                                              cGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGa                    <  1:151422/148‑1 (MQ=255)
                                                                                                                              cGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGa                    <  2:3241841/148‑1 (MQ=255)
                                                                                                                               gATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGaa                   >  2:2885517/1‑148 (MQ=255)
                                                                                                                               gATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGaa                   >  1:794520/1‑148 (MQ=255)
                                                                                                                                  gAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGATTGATAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAAttt                >  1:363398/1‑148 (MQ=255)
                                                                                                                                  gAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAAttt                >  1:2859954/1‑148 (MQ=255)
                                                                                                                                   aaTTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTa               <  1:1857807/148‑1 (MQ=255)
                                                                                                                                    aTTTGCATACACTTATAACGAATTTATAATTACGCGTATTTTTTCCGCAATTCAGATATCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTTATAGTTGAATTTaa              <  1:2444757/148‑1 (MQ=255)
                                                                                                                                     tttGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAAc             <  2:2404492/148‑1 (MQ=255)
                                                                                                                                     tttGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAAAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAAc             <  1:1323158/148‑1 (MQ=255)
                                                                                                                                       tGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACtt           <  2:575519/148‑1 (MQ=255)
                                                                                                                                           gACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCa       <  1:1955417/148‑1 (MQ=255)
                                                                                                                                           gACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCa       <  2:1633235/148‑1 (MQ=255)
                                                                                                                                            acacTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCAt      <  1:1679438/148‑1 (MQ=255)
                                                                                                                                            acacTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCAt      >  2:1805432/1‑148 (MQ=255)
                                                                                                                                                tGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGCtt  >  2:3478032/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
ACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGCTT  >  CP009273/4567133‑4567424

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