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.5% 834.3 / 15.4 270I166I (ATC→ATA) mcrC5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrC
Reads supporting (aligned to +/- strand):  ref base G (8/2);  new base T (157/103);  total (165/105)
Fisher's exact test for biased strand distribution p-value = 3.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 271 total aligned reads displayed.

CGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGC  >  CP009273/4567134‑4567422
                                                                                                                                                |                                                                                                                                                
 gTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCGTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATaa                                                                                                                                              >  1:2109787/1‑148 (MQ=255)
 gTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATaa                                                                                                                                              <  2:301078/148‑1 (MQ=255)
     aTTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGaa                                                                                                                                          >  1:936643/1‑148 (MQ=255)
        cACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAAttt                                                                                                                                       <  2:649035/148‑1 (MQ=255)
         acgacgGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTta                                                                                                                                      <  2:1278381/148‑1 (MQ=255)
               gCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTAtaataa                                                                                                                                <  1:1136124/148‑1 (MQ=255)
                 aaaaTTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAAcg                                                                                                                              <  2:1201741/148‑1 (MQ=255)
                 aaaaTTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAAcg                                                                                                                              >  2:1906607/1‑148 (MQ=255)
                   aaTTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAAcgcg                                                                                                                            >  1:2180473/1‑148 (MQ=255)
                   aaTTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATGTATAATAAcgcg                                                                                                                            >  1:1820200/1‑148 (MQ=255)
                    aTTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGt                                                                                                                           <  1:3440397/148‑1 (MQ=255)
                     ttCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTa                                                                                                                          >  2:585963/1‑148 (MQ=255)
                      tCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATACTAACGCGTAt                                                                                                                         >  1:2331455/1‑148 (MQ=255)
                      tCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTAt                                                                                                                         <  2:1749340/148‑1 (MQ=255)
                         tAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTAtttt                                                                                                                      >  2:2613836/1‑148 (MQ=255)
                          aaaGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTAttttt                                                                                                                     >  2:2888228/1‑148 (MQ=255)
                          aaaGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTAttttt                                                                                                                     >  2:1114066/1‑148 (MQ=255)
                           aaGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTAtttttt                                                                                                                    >  2:2657247/1‑148 (MQ=255)
                            aGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGAAGACACTTATAACGAATTTATAATAACGCGTATTTTTTc                                                                                                                   >  2:2791430/1‑148 (MQ=255)
                              aaaCTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCg                                                                                                                 >  2:2350270/1‑148 (MQ=255)
                              aaaCTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCg                                                                                                                 >  2:2935246/1‑148 (MQ=255)
                               aaCTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGc                                                                                                                >  1:3232704/1‑148 (MQ=255)
                                 cTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCa                                                                                                              >  2:506400/1‑148 (MQ=255)
                                     tGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCa                                                                                                          >  2:3362089/1‑148 (MQ=255)
                                      gATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAg                                                                                                         >  1:2525159/1‑148 (MQ=255)
                                       aTAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGa                                                                                                        <  1:1181189/148‑1 (MQ=255)
                                        tAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGAt                                                                                                       >  1:2314956/1‑148 (MQ=255)
                                          aaGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAg                                                                                                     <  2:956192/148‑1 (MQ=255)
                                            gTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCt                                                                                                   >  1:591830/1‑148 (MQ=255)
                                              aaTGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGa                                                                                                 >  2:2451957/1‑148 (MQ=255)
                                               aTGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGaa                                                                                                >  2:3452751/1‑148 (MQ=255)
                                               aTGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGaa                                                                                                >  1:431038/1‑148 (MQ=255)
                                                  aCATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAAtg                                                                                             <  2:1018876/148‑1 (MQ=255)
                                                  aCATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAAtg                                                                                             <  1:127792/148‑1 (MQ=255)
                                                   cATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAAtgc                                                                                            >  1:1246404/1‑148 (MQ=255)
                                                   cATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAAtgc                                                                                            >  1:3017508/1‑148 (MQ=255)
                                                      ctctTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAAtgctgc                                                                                         >  2:2066105/1‑148 (MQ=255)
                                                        ctTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCgg                                                                                       >  1:2809382/1‑148 (MQ=255)
                                                         tttttCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGa                                                                                      >  2:1053387/1‑148 (MQ=255)
                                                            ttCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGtt                                                                                   <  1:2996180/148‑1 (MQ=255)
                                                              cGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATGTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTaa                                                                                 <  2:3257886/148‑1 (MQ=255)
                                                               gTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTaaa                                                                                >  2:753339/1‑148 (MQ=255)
                                                               gTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTaaa                                                                                >  2:2010703/1‑148 (MQ=255)
                                                                tttctttcAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAAt                                                                               >  2:545271/1‑148 (MQ=255)
                                                                tttctttcAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAAt                                                                               >  1:1033683/1‑148 (MQ=255)
                                                                 ttctttcAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATg                                                                              >  2:3443985/1‑148 (MQ=255)
                                                                      tcAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAagag                                                                         >  1:1396151/1‑148 (MQ=255)
                                                                        aaaaTCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTg                                                                       <  1:2619067/148‑1 (MQ=255)
                                                                         aaaTCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGc                                                                      >  2:2369409/1‑148 (MQ=255)
                                                                          aaTCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATCACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCt                                                                     >  2:1820771/1‑148 (MQ=255)
                                                                            tCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTaa                                                                   <  1:2164334/148‑1 (MQ=255)
                                                                            tCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTaa                                                                   >  2:2731396/1‑148 (MQ=255)
                                                                              aTAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATg                                                                 <  2:1999909/148‑1 (MQ=255)
                                                                               tAGAAACGGTAGTGTCCTTTGTTTTGACCTGGGATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGc                                                                >  1:3397852/1‑148 (MQ=255)
                                                                               tAGAAACGGTAGTGTCCTTTGTTTTGACCTGGGATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGc                                                                >  1:1638532/1‑148 (MQ=255)
                                                                                aGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGcc                                                               >  2:2548640/1‑148 (MQ=255)
                                                                                aGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGcc                                                               >  2:2457732/1‑148 (MQ=255)
                                                                                  aaaCGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCg                                                             >  1:1738886/1‑148 (MQ=255)
                                                                                  aaaCGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCg                                                             >  2:450450/1‑148 (MQ=255)
                                                                                    aCGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGt                                                           >  2:815906/1‑148 (MQ=255)
                                                                                      ggTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTaa                                                         >  1:3549210/1‑148 (MQ=255)
                                                                                      ggTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTaa                                                         <  1:2449088/148‑1 (MQ=255)
                                                                                       gTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAAt                                                        >  1:2245581/1‑148 (MQ=255)
                                                                                         aGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAAttt                                                      >  2:3095635/1‑148 (MQ=255)
                                                                                          gtgtCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAAtttt                                                     >  1:839060/1‑148 (MQ=255)
                                                                                              ccTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCtat                                                 >  1:3162139/1‑148 (MQ=255)
                                                                                                tttGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATaa                                               >  1:2864039/1‑148 (MQ=255)
                                                                                                tttGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATaa                                               >  1:3179992/1‑148 (MQ=255)
                                                                                                 ttGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATaaa                                              >  1:1390533/1‑148 (MQ=255)
                                                                                                  tGTTTTGACCTGGAATAGTATTATTGACGATGAATTTGCGGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAg                                             <  2:2395222/148‑1 (MQ=255)
                                                                                                   gTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGt                                            <  2:2845411/148‑1 (MQ=255)
                                                                                                     tttGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGa                                          <  1:3277231/148‑1 (MQ=255)
                                                                                                       tGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAAc                                        >  2:985081/1‑148 (MQ=255)
                                                                                                         aCCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGa                                      <  1:1533698/148‑1 (MQ=255)
                                                                                                          ccTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAg                                     <  2:433115/148‑1 (MQ=255)
                                                                                                           cTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATCACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGc                                    >  1:99518/1‑148 (MQ=255)
                                                                                                               aaTAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAAAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCa                                >  1:717693/1‑148 (MQ=255)
                                                                                                                  aGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCAtct                             <  2:837119/148‑1 (MQ=255)
                                                                                                                  aGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCAtct                             >  1:2604799/1‑148 (MQ=255)
                                                                                                                  aGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCAtct                             >  1:2753026/1‑148 (MQ=255)
                                                                                                                    aaTTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATctct                           >  2:1260752/1‑148 (MQ=255)
                                                                                                                    aaTTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATctct                           >  2:3458122/1‑148 (MQ=255)
                                                                                                                      ttattGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGa                         <  2:2755720/148‑1 (MQ=255)
                                                                                                                      ttattGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGa                         <  2:1711763/148‑1 (MQ=255)
                                                                                                                       tattGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGAt                        >  1:649034/1‑148 (MQ=255)
                                                                                                                        attGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATa                       >  1:2414251/1‑148 (MQ=255)
                                                                                                                        attGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATa                       >  1:628538/1‑148 (MQ=255)
                                                                                                                         ttGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAg                      >  2:2096501/1‑148 (MQ=255)
                                                                                                                                 gAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAAttt              >  1:1102664/1‑148 (MQ=255)
                                                                                                                                 gAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAAttt              >  1:1374018/1‑148 (MQ=255)
                                                                                                                                 gAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAAttt              >  1:2573317/1‑148 (MQ=255)
                                                                                                                                 gAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAAttt              <  2:2573317/148‑1 (MQ=255)
                                                                                                                                  aaTTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTa             >  1:3136973/1‑148 (MQ=255)
                                                                                                                                    tttGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAAc           <  1:517264/148‑1 (MQ=255)
                                                                                                                                    tttGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAAc           >  1:969925/1‑148 (MQ=255)
                                                                                                                                      tGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACtt         <  1:1967463/148‑1 (MQ=255)
                                                                                                                                         aGACTCTTATAAAGAATTTATAATATCGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTc      <  1:1926796/148‑1 (MQ=255)
                                                                                                                                          gACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCa     >  2:3259413/1‑148 (MQ=255)
                                                                                                                                           acacTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCAt    <  2:3183695/148‑1 (MQ=255)
                                                                                                                                           acacTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCAt    >  1:1339412/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGC  >  CP009273/4567134‑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: ‑