Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,271,128 A→C 100% T145P (ACC→CCC)  yhaV → toxin of the SohB(PrlF)‑YhaV toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,271,1280AC97.1% 941.0 / 4.3 235T145P (ACC→CCC) yhaVtoxin of the SohB(PrlF)‑YhaV toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base A (6/1);  new base C (129/99);  total (135/100)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 236 total aligned reads displayed.

TTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAACCCCAGCCTTTCGCAGTCCTTCCAGGCTAT  >  CP009273/3270984‑3271270
                                                                                                                                                |                                                                                                                                              
tttCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAcccc                                                                                                                                             <  2:3225335/148‑1 (MQ=255)
 ttCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCt                                                                                                                                            >  1:2453835/1‑148 (MQ=255)
   cGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCa                                                                                                                                          >  1:3794386/1‑148 (MQ=255)
     ttATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCAcc                                                                                                                                        >  1:396171/1‑148 (MQ=255)
         aGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCgag                                                                                                                                    >  2:2315229/1‑148 (MQ=255)
           tGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGaa                                                                                                                                  <  1:350447/148‑1 (MQ=255)
            gAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGaaa                                                                                                                                 <  1:1492037/148‑1 (MQ=255)
            gAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGaaa                                                                                                                                 <  1:1552825/148‑1 (MQ=255)
            gAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGaaa                                                                                                                                 <  2:3641281/148‑1 (MQ=255)
            gAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCAGCCGACTGGGAACCCCTCACCCGAGaaa                                                                                                                                 <  2:3509461/148‑1 (MQ=255)
               aaaGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGACCCCCTCACCCGAGAAACag                                                                                                                              >  2:2593841/1‑148 (MQ=255)
               aaaGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACag                                                                                                                              >  1:1321193/1‑148 (MQ=255)
               aaaGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACag                                                                                                                              >  2:2647644/1‑148 (MQ=255)
                 agagAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACagaa                                                                                                                            >  1:1389502/1‑148 (MQ=255)
                  gagaAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACagaag                                                                                                                           >  2:1127549/1‑148 (MQ=255)
                   agaAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACagaaga                                                                                                                          <  2:1208105/148‑1 (MQ=255)
                       aGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAAcc                                                                                                                      <  1:3333947/148‑1 (MQ=255)
                       aGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAAcc                                                                                                                      <  2:21780/148‑1 (MQ=255)
                        gTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAAccc                                                                                                                     >  2:2012382/1‑148 (MQ=255)
                         tcatcaTTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCa                                                                                                                    >  1:2330500/1‑148 (MQ=255)
                         tcatcaTTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCa                                                                                                                    <  1:2520869/148‑1 (MQ=255)
                         tcatcaTTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCa                                                                                                                    >  1:3708769/1‑148 (MQ=255)
                          catcatTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCAt                                                                                                                   <  1:2051990/148‑1 (MQ=255)
                           atcatTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCAtt                                                                                                                  >  1:1742659/1‑148 (MQ=255)
                            tcatTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTg                                                                                                                 <  2:202006/148‑1 (MQ=255)
                            tcatTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTg                                                                                                                 >  2:3020336/1‑148 (MQ=255)
                             catTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGa                                                                                                                >  2:3441963/1‑148 (MQ=255)
                               ttCTTGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATTCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATg                                                                                                              <  1:622226/148‑1 (MQ=255)
                               ttCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATg                                                                                                              <  1:449999/148‑1 (MQ=255)
                               ttCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATg                                                                                                              <  1:170038/148‑1 (MQ=255)
                                   gggATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGgtgt                                                                                                          >  2:1840154/1‑148 (MQ=255)
                                   gggATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGgtgt                                                                                                          <  2:327996/148‑1 (MQ=255)
                                     gatggatgAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCACTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTc                                                                                                        <  1:2832211/148‑1 (MQ=255)
                                      atggatgAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTca                                                                                                       >  1:712290/1‑148 (MQ=255)
                                         gatgAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACaa                                                                                                    >  1:2822087/1‑148 (MQ=255)
                                             aaCGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGcc                                                                                                >  1:1333695/1‑148 (MQ=255)
                                               cGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCccg                                                                                              >  2:2365262/1‑148 (MQ=255)
                                                   gAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGa                                                                                          >  1:2393882/1‑148 (MQ=255)
                                                   gAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGa                                                                                          <  1:194723/148‑1 (MQ=255)
                                                      aaCACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACg                                                                                       >  1:2658846/1‑148 (MQ=255)
                                                        cacTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCa                                                                                     >  1:1609015/1‑148 (MQ=255)
                                                         acTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCAt                                                                                    <  2:2544032/148‑1 (MQ=255)
                                                          ctctGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCAtt                                                                                   >  2:13539/1‑148 (MQ=255)
                                                          ctctGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCAtt                                                                                   >  2:13471/1‑148 (MQ=255)
                                                            ctGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCg                                                                                 >  2:3153064/1‑148 (MQ=255)
                                                            ctGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCg                                                                                 <  1:2593881/148‑1 (MQ=255)
                                                            ctGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCg                                                                                 <  1:2593717/148‑1 (MQ=255)
                                                             tGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTTGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCgt                                                                                <  2:1025006/148‑1 (MQ=255)
                                                              gcgcACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCgtg                                                                               <  2:458306/148‑1 (MQ=255)
                                                                 cACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTg                                                                            <  1:3135900/148‑1 (MQ=255)
                                                                  aCCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGa                                                                           <  2:798747/148‑1 (MQ=255)
                                                                  aCCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGa                                                                           <  1:2105918/148‑1 (MQ=255)
                                                                  aCCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGa                                                                           <  1:2160915/148‑1 (MQ=255)
                                                                  aCCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGa                                                                           <  2:798846/148‑1 (MQ=255)
                                                                   ccTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGAc                                                                          <  2:55228/148‑1 (MQ=255)
                                                                     tACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCa                                                                        >  1:412871/1‑148 (MQ=255)
                                                                        ggTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCg                                                                     <  2:3805037/148‑1 (MQ=255)
                                                                        ggTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCg                                                                     <  1:53861/148‑1 (MQ=255)
                                                                           aaaaaaaCAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGc                                                                  >  1:989584/1‑148 (MQ=255)
                                                                                  cAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCttt                                                           <  1:1018771/148‑1 (MQ=255)
                                                                                  cAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCttt                                                           <  2:3765769/148‑1 (MQ=255)
                                                                                    gATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTtct                                                         >  1:2326542/1‑148 (MQ=255)
                                                                                           ataCCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTccc                                                  >  2:2070090/1‑148 (MQ=255)
                                                                                           ataCCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTccc                                                  <  1:2631471/148‑1 (MQ=255)
                                                                                           ataCCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTccc                                                  >  2:3407835/1‑148 (MQ=255)
                                                                                               cGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGAc                                              >  2:1492420/1‑148 (MQ=255)
                                                                                                 tATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCa                                            >  1:1775172/1‑148 (MQ=255)
                                                                                                 tATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCa                                            >  1:3664619/1‑148 (MQ=255)
                                                                                                   ttCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGa                                          <  2:1660071/148‑1 (MQ=255)
                                                                                                    tCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGaa                                         <  1:594727/148‑1 (MQ=255)
                                                                                                    tCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGaa                                         >  1:1871895/1‑148 (MQ=255)
                                                                                                        cAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCAc                                     >  2:1482256/1‑148 (MQ=255)
                                                                                                        cAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCAc                                     >  2:2914024/1‑148 (MQ=255)
                                                                                                              gTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAAc                               >  2:907133/1‑148 (MQ=255)
                                                                                                              gTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAAc                               >  1:2775594/1‑148 (MQ=255)
                                                                                                               ttAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAAcc                              <  2:3387738/148‑1 (MQ=255)
                                                                                                                tAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAAccc                             >  1:2016515/1‑148 (MQ=255)
                                                                                                                  aaaaGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAACCCCa                           >  2:1841919/1‑148 (MQ=255)
                                                                                                                     agagGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAACCCCAGcc                        >  2:1105404/1‑148 (MQ=255)
                                                                                                                     agagGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAACCCCAGcc                        >  1:238239/1‑148 (MQ=255)
                                                                                                                      gagGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAACCCCAGCCt                       >  2:3417215/1‑148 (MQ=255)
                                                                                                                      gagGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAACCCCAGCCt                       >  1:2718751/1‑148 (MQ=255)
                                                                                                                       agGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAACCCCAGCCtt                      >  2:845532/1‑148 (MQ=255)
                                                                                                                        ggACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAACCCCAGCCttt                     <  1:554826/148‑1 (MQ=255)
                                                                                                                        ggACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAACCCCAGCCttt                     <  1:1999657/148‑1 (MQ=255)
                                                                                                                            aTCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAACCCCAGCCTTTCGCa                 >  1:2473704/1‑148 (MQ=255)
                                                                                                                            aTCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAACCCCAGCCTTTCGCa                 >  2:884045/1‑148 (MQ=255)
                                                                                                                             tcctccTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAACCCCAGCCTTTCGCAg                >  1:2197186/1‑148 (MQ=255)
                                                                                                                              cctcctGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACACTGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAACCCCAGCCTTTCGCAGt               <  2:2677121/148‑1 (MQ=255)
                                                                                                                              cctcctGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAACCCCAGCCTTTCGCAGt               <  2:371022/148‑1 (MQ=255)
                                                                                                                               ctcctGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAACCCCAGCCTTTCGCAGTc              <  1:1461982/148‑1 (MQ=255)
                                                                                                                               ctcctGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCAAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAACCCCAGCCTTTCGCAGTc              <  2:2825166/148‑1 (MQ=255)
                                                                                                                               ctcctGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCAAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAACCCCAGCCTTTCGCAGTc              <  2:1329732/148‑1 (MQ=255)
                                                                                                                                 cctGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAACCCCAGCCTTTCGCAGTCCt            >  1:169563/1‑148 (MQ=255)
                                                                                                                                 cctGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAACCCCAGCCTTTCGCAGTCCt            >  2:1011917/1‑148 (MQ=255)
                                                                                                                                  ctGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAACCCCAGCCTTTCGCAGTCCtt           <  1:3173750/148‑1 (MQ=255)
                                                                                                                                   tGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAACCCCAGCCTTTCGCAGTCCTTc          <  1:2344322/148‑1 (MQ=255)
                                                                                                                                       gACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAACCCCAGCCTTTCGCAGTCCTTCCAgg      >  1:2710396/1‑148 (MQ=255)
                                                                                                                                        aCTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAACCCCAGCCTTTCGCAGTCCTTCCAGGc     >  2:2723464/1‑148 (MQ=255)
                                                                                                                                           gggAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAACCCCAGCCTTTCGCAGTCCTTCCAGGCTAt  <  1:3270530/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                              
TTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAACCCCAGCCTTTCGCAGTCCTTCCAGGCTAT  >  CP009273/3270984‑3271270

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: