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% 837.0 / 16.9 272I166I (ATC→ATAmcrC5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrC
Reads supporting (aligned to +/- strand):  ref base G (3/6);  new base T (123/140);  total (126/146)
Fisher's exact test for biased strand distribution p-value = 5.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 274 total aligned reads displayed.

ACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGC  >  CP009273/4567133‑4567422
                                                                                                                                                 |                                                                                                                                                
aCGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTAt                                                                                                                                                >  2:2136248/1‑148 (MQ=255)
 cGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTtata                                                                                                                                               <  2:2094581/148‑1 (MQ=255)
     aaTTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGa                                                                                                                                           <  1:2375889/148‑1 (MQ=255)
      aTTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGaa                                                                                                                                          <  1:334800/148‑1 (MQ=255)
       ttCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAAt                                                                                                                                         <  1:2631381/148‑1 (MQ=255)
       ttCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAAt                                                                                                                                         <  2:142979/148‑1 (MQ=255)
         cACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAAttt                                                                                                                                       <  2:2603534/148‑1 (MQ=255)
           cgacgGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTtat                                                                                                                                     >  1:2620129/1‑148 (MQ=255)
              cgGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTataat                                                                                                                                  <  2:2549764/148‑1 (MQ=255)
                gCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTAtaataa                                                                                                                                <  2:1961267/148‑1 (MQ=255)
                    aaTTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAAcgcg                                                                                                                            <  2:801221/148‑1 (MQ=255)
                     aTTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGt                                                                                                                           >  1:1268240/1‑148 (MQ=255)
                      ttCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTa                                                                                                                          <  2:911768/148‑1 (MQ=255)
                       tCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTAt                                                                                                                         >  2:682412/1‑148 (MQ=255)
                       tCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTAt                                                                                                                         <  1:116983/148‑1 (MQ=255)
                        cATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTAtt                                                                                                                        <  1:279991/148‑1 (MQ=255)
                           aaaGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACTCGTAttttt                                                                                                                     >  1:2549609/1‑148 (MQ=255)
                            aaGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTAtttttt                                                                                                                    <  1:2153664/148‑1 (MQ=255)
                            aaGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTAtttttt                                                                                                                    <  1:2058997/148‑1 (MQ=255)
                             aGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTc                                                                                                                   <  2:1678539/148‑1 (MQ=255)
                             aGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTc                                                                                                                   <  1:2990681/148‑1 (MQ=255)
                              gAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTcc                                                                                                                  >  2:2683814/1‑148 (MQ=255)
                                 aCTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGcc                                                                                                               >  1:2236139/1‑148 (MQ=255)
                                    tttGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCAtt                                                                                                            >  2:2721124/1‑148 (MQ=255)
                                    tttGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCAtt                                                                                                            <  1:2298237/148‑1 (MQ=255)
                                     ttGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTc                                                                                                           <  2:313713/148‑1 (MQ=255)
                                      tGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCa                                                                                                          >  2:567378/1‑148 (MQ=255)
                                       gATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAg                                                                                                         <  1:1501166/148‑1 (MQ=255)
                                        aTAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGa                                                                                                        <  1:2680714/148‑1 (MQ=255)
                                         tAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGAt                                                                                                       <  2:479333/148‑1 (MQ=255)
                                         tAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGAt                                                                                                       <  2:1140044/148‑1 (MQ=255)
                                          aaaGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATa                                                                                                      >  2:2498896/1‑148 (MQ=255)
                                          aaaGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATa                                                                                                      >  1:981088/1‑148 (MQ=255)
                                          aaaGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATa                                                                                                      >  2:1792363/1‑148 (MQ=255)
                                            aGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGc                                                                                                    >  1:302167/1‑148 (MQ=255)
                                                aTGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCAATTCAGATAGCTGaa                                                                                                >  2:3366071/1‑148 (MQ=255)
                                                   aCATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAAtg                                                                                             <  2:2885777/148‑1 (MQ=255)
                                                   aCATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAAtg                                                                                             >  1:3069566/1‑148 (MQ=255)
                                                    cATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAAtgc                                                                                            <  1:2011634/148‑1 (MQ=255)
                                                     aTCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAAtgct                                                                                           <  2:808701/148‑1 (MQ=255)
                                                            tttCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGt                                                                                    <  1:2330296/148‑1 (MQ=255)
                                                            tttCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGt                                                                                    <  1:1792065/148‑1 (MQ=255)
                                                              tCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTa                                                                                  >  2:650506/1‑148 (MQ=255)
                                                                 tttctttcAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAAt                                                                               >  2:229492/1‑148 (MQ=255)
                                                                 tttctttcAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAAt                                                                               >  1:248441/1‑148 (MQ=255)
                                                                 tttctttcAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAAt                                                                               >  2:1053932/1‑148 (MQ=255)
                                                                 tttctttcAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAAt                                                                               >  1:1077925/1‑148 (MQ=255)
                                                                  ttctttcAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATg                                                                              >  1:242601/1‑148 (MQ=255)
                                                                   tctttcAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGa                                                                             <  2:611451/148‑1 (MQ=255)
                                                                      ttcAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAaga                                                                          <  1:2183861/148‑1 (MQ=255)
                                                                        cAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGt                                                                        <  1:1018003/148‑1 (MQ=255)
                                                                         aaaaTCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTg                                                                       >  2:2626443/1‑148 (MQ=255)
                                                                         aaaaTCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTg                                                                       <  1:2659368/148‑1 (MQ=255)
                                                                          aaaTCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGc                                                                      >  1:891041/1‑148 (MQ=255)
                                                                          aaaTCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGc                                                                      <  2:768867/148‑1 (MQ=255)
                                                                          aaaTCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGc                                                                      <  1:2545609/148‑1 (MQ=255)
                                                                          aaaTCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGc                                                                      <  1:3100114/148‑1 (MQ=255)
                                                                            aTCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTa                                                                    >  1:2828624/1‑148 (MQ=255)
                                                                             tCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTaa                                                                   >  2:2411869/1‑148 (MQ=255)
                                                                                tAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGc                                                                <  1:3195823/148‑1 (MQ=255)
                                                                                  gAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGccc                                                              >  1:3225634/1‑148 (MQ=255)
                                                                                    aaCGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCgg                                                            <  1:2995738/148‑1 (MQ=255)
                                                                                      cGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTa                                                          <  2:2355307/148‑1 (MQ=255)
                                                                                         tAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAAtt                                                       >  1:399352/1‑148 (MQ=255)
                                                                                           gtgtCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAAtttt                                                     <  2:3011648/148‑1 (MQ=255)
                                                                                           gtgtCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAAtttt                                                     <  1:567378/148‑1 (MQ=255)
                                                                                            tgtCCTTTGTTTTGACCTGGAATAGAATTATTGACGATAAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTc                                                    >  2:2841469/1‑148 (MQ=255)
                                                                                             gtCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCt                                                   <  1:437342/148‑1 (MQ=255)
                                                                                              tCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCta                                                  >  2:2133866/1‑148 (MQ=255)
                                                                                                 tttGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATaa                                               <  1:650506/148‑1 (MQ=255)
                                                                                                 tttGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATaa                                               <  1:700170/148‑1 (MQ=255)
                                                                                                 tttGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATaa                                               <  2:2534740/148‑1 (MQ=255)
                                                                                                      tttGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGa                                          >  2:964024/1‑148 (MQ=255)
                                                                                                       ttGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGaa                                         >  1:3426499/1‑148 (MQ=255)
                                                                                                       ttGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGaa                                         <  2:2438636/148‑1 (MQ=255)
                                                                                                          aCCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGa                                      <  2:248446/148‑1 (MQ=255)
                                                                                                           ccTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAg                                     <  2:3263596/148‑1 (MQ=255)
                                                                                                            cTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGc                                    <  1:14618/148‑1 (MQ=255)
                                                                                                            cTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGc                                    <  1:2102543/148‑1 (MQ=255)
                                                                                                              ggAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCtt                                  <  1:3563608/148‑1 (MQ=255)
                                                                                                              ggAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATNCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCtt                                  >  1:1807027/1‑148 (MQ=255)
                                                                                                                   aGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCAtct                             <  1:3536005/148‑1 (MQ=255)
                                                                                                                      attattGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTg                          <  2:3436318/148‑1 (MQ=255)
                                                                                                                        tattGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGAt                        >  2:1144444/1‑148 (MQ=255)
                                                                                                                          ttGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAg                      <  2:131880/148‑1 (MQ=255)
                                                                                                                          ttGAAGATGAATTTGCAGACACTTATAAAGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAg                      <  2:3564417/148‑1 (MQ=255)
                                                                                                                           tGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGt                     <  2:2632577/148‑1 (MQ=255)
                                                                                                                             aCGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTg                   >  2:2186305/1‑148 (MQ=255)
                                                                                                                              cGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGa                  <  1:1363633/148‑1 (MQ=255)
                                                                                                                               gATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGaa                 <  1:2509179/148‑1 (MQ=255)
                                                                                                                                   aaTTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTa             <  1:258938/148‑1 (MQ=255)
                                                                                                                                   aaTTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTa             <  1:2503837/148‑1 (MQ=255)
                                                                                                                                          aGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTc      <  2:3544496/148‑1 (MQ=255)
                                                                                                                                          aGACACTTATAACGAAATTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTc      <  2:2929474/148‑1 (MQ=255)
                                                                                                                                           gACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAAATTTTCa     >  2:2395987/1‑148 (MQ=255)
                                                                                                                                             cacTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATg   <  2:2507550/148‑1 (MQ=255)
                                                                                                                                             cacTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTACTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATg   >  1:779578/1‑148 (MQ=255)
                                                                                                                                              acTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGc  >  2:562837/1‑148 (MQ=255)
                                                                                                                                              acTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGc  <  1:1938751/148‑1 (MQ=255)
                                                                                                                                              acTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGc  <  1:2977323/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
ACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGC  >  CP009273/4567133‑4567422

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: