Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,025,678 A→T 27.9% Q318L (CAA→CTA)  uacT → uric acid permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,025,6780AT27.9% 502.7 / 253.2 280Q318L (CAA→CTA) uacTuric acid permease
Reads supporting (aligned to +/- strand):  ref base A (91/110);  new base T (39/39);  total (130/150)
Fisher's exact test for biased strand distribution p-value = 5.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

Only 100 of 283 total aligned reads displayed.

TCGAGTCAATGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAG  >  CP009273/3025531‑3025827
                                                                                                                                                   |                                                                                                                                                     
tCGAGTCAATGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCaaa                                                                                                                                                     <  1:3081567/150‑1 (MQ=255)
  gAGTCAATGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAAc                                                                                                                                                   <  2:231072/150‑1 (MQ=255)
  gAGTCAATGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAAc                                                                                                                                                   >  1:1691165/1‑150 (MQ=255)
      cAATGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTg                                                                                                                                               <  2:2377907/150‑1 (MQ=255)
      cAATGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTg                                                                                                                                               <  2:3487311/150‑1 (MQ=255)
      cAATGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTg                                                                                                                                               <  2:8766/150‑1 (MQ=255)
       aaTGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTgg                                                                                                                                              >  2:3985450/1‑150 (MQ=255)
        aTGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGc                                                                                                                                             >  1:1033953/1‑150 (MQ=255)
          ggggATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCt                                                                                                                                           >  2:4203639/1‑150 (MQ=255)
           gggATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTg                                                                                                                                          >  1:1835848/1‑150 (MQ=255)
           gggATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTg                                                                                                                                          <  2:3551797/150‑1 (MQ=255)
            ggATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTgg                                                                                                                                         <  2:3471382/150‑1 (MQ=255)
               tGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTa                                                                                                                                      >  2:4282028/1‑150 (MQ=255)
                         aCTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACCGCTTCCCCCACACGTCATTGTCTCAAAACGTTGGCCTGGTTAGCGTGAcgcg                                                                                                                            >  2:1671368/1‑150 (MQ=255)
                          cTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGAcgcgc                                                                                                                           >  2:2146996/1‑150 (MQ=255)
                            gggTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGt                                                                                                                         <  2:4109653/150‑1 (MQ=255)
                            gggTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGt                                                                                                                         <  2:1831889/150‑1 (MQ=255)
                             ggTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGtt                                                                                                                        >  1:4118137/1‑150 (MQ=255)
                                gAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCAt                                                                                                                     >  2:750133/1‑150 (MQ=255)
                                gAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCAt                                                                                                                     >  2:1019822/1‑150 (MQ=255)
                                 aaaTAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATa                                                                                                                    >  2:4402529/1‑150 (MQ=255)
                                    tAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGcc                                                                                                                 >  1:2351755/1‑150 (MQ=255)
                                    tAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGcc                                                                                                                 >  1:952471/1‑150 (MQ=255)
                                    tAGTCGGTCGGAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGcc                                                                                                                 <  2:2818519/150‑1 (MQ=255)
                                         ggtcgTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTgg                                                                                                            <  2:2191674/150‑1 (MQ=255)
                                            cgTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGgtg                                                                                                         >  1:4264926/1‑150 (MQ=255)
                                            cgTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGgtg                                                                                                         <  1:570941/150‑1 (MQ=255)
                                             gTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGgtgt                                                                                                        >  1:1055187/1‑150 (MQ=255)
                                              tAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGgtgtg                                                                                                       <  1:3852801/150‑1 (MQ=255)
                                               aaaCTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGgtgtgt                                                                                                      <  2:3779718/150‑1 (MQ=255)
                                                   tctctTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACAAGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTAttt                                                                                                  <  2:1554226/150‑1 (MQ=255)
                                                    ctctTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATttc                                                                                                 <  2:1199133/150‑1 (MQ=255)
                                                    ctctTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATttc                                                                                                 <  2:1429534/150‑1 (MQ=255)
                                                     tctTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATttct                                                                                                >  1:3685426/1‑150 (MQ=255)
                                                     tctTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATttct                                                                                                >  1:4079688/1‑150 (MQ=255)
                                                      ctTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATttctt                                                                                               >  1:724584/1‑150 (MQ=255)
                                                       ttCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATttcttc                                                                                              <  1:2047531/150‑1 (MQ=255)
                                                       ttCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATttcttc                                                                                              <  1:3195522/150‑1 (MQ=255)
                                                        tCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCg                                                                                             <  2:2532386/150‑1 (MQ=255)
                                                         cGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCgg                                                                                            >  1:4196266/1‑150 (MQ=255)
                                                         cGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCgg                                                                                            >  1:3262170/1‑150 (MQ=255)
                                                          gCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCggg                                                                                           <  1:1483325/150‑1 (MQ=255)
                                                          gCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCggg                                                                                           <  1:3290586/150‑1 (MQ=255)
                                                           cACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGa                                                                                          <  2:2767981/150‑1 (MQ=255)
                                                            aCGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGaa                                                                                         >  2:2470048/1‑150 (MQ=255)
                                                              gATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAAtt                                                                                       <  1:3143179/150‑1 (MQ=255)
                                                               atatTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATta                                                                                      >  1:3016130/1‑150 (MQ=255)
                                                                    attCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTa                                                                                 >  1:3926077/1‑150 (MQ=255)
                                                                        gcgGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATcc                                                                             <  1:2581464/150‑1 (MQ=255)
                                                                         cgGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCt                                                                            <  2:4486227/150‑1 (MQ=255)
                                                                         cgGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCt                                                                            >  1:1579592/1‑150 (MQ=255)
                                                                          gggCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTg                                                                           >  2:1310218/1‑150 (MQ=255)
                                                                             cTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCTCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGttt                                                                        <  2:370598/150‑1 (MQ=255)
                                                                             cTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGttt                                                                        >  1:1373644/1‑150 (MQ=255)
                                                                             cTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGttt                                                                        >  1:3658574/1‑150 (MQ=255)
                                                                             cTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGttt                                                                        <  1:3775534/150‑1 (MQ=255)
                                                                             cTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGttt                                                                        <  2:3083333/150‑1 (MQ=255)
                                                                              tGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTg                                                                       <  1:1914391/150‑1 (MQ=255)
                                                                               gCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTgg                                                                      <  1:1576494/150‑1 (MQ=255)
                                                                               gCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTgg                                                                      <  1:4362696/150‑1 (MQ=255)
                                                                                  tgtCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCAt                                                                   >  2:4522051/1‑150 (MQ=255)
                                                                                  tgtCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCAt                                                                   >  2:539805/1‑150 (MQ=255)
                                                                                  tgtCGAAGGCGTAGGGACAATGATAGGCGGCACGATTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTATCGTGACGCGCGTACATAGCCGCATGGTGTGAATTTCTTCGGGAATTATAATAATCCTGTTTGGCAt                                                                   >  2:4521767/1‑150 (MQ=255)
                                                                                       aTGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGc                                                              <  2:1033953/150‑1 (MQ=255)
                                                                                         ggCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCa                                                            <  1:1779654/150‑1 (MQ=255)
                                                                                          gCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCaa                                                           <  2:4477998/150‑1 (MQ=255)
                                                                                          gCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCaa                                                           <  2:881549/150‑1 (MQ=255)
                                                                                             tAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCaaaaa                                                        <  1:3225150/150‑1 (MQ=255)
                                                                                              aGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAAt                                                       <  1:64354/150‑1 (MQ=255)
                                                                                                ggACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATgg                                                     >  2:3277647/1‑150 (MQ=255)
                                                                                                 gACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGc                                                    <  1:4427564/150‑1 (MQ=255)
                                                                                                    aaTGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGt                                                 >  1:1130901/1‑150 (MQ=255)
                                                                                                     aTGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTg                                                >  2:2670784/1‑150 (MQ=255)
                                                                                                     aTGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTg                                                <  1:4185278/150‑1 (MQ=255)
                                                                                                     aTGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTg                                                <  1:2572326/150‑1 (MQ=255)
                                                                                                      tGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGc                                               <  1:2293738/150‑1 (MQ=255)
                                                                                                          aGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGt                                           >  2:3194992/1‑150 (MQ=255)
                                                                                                             cggcACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGc                                        >  1:4061889/1‑150 (MQ=255)
                                                                                                                cACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTc                                     >  1:4488753/1‑150 (MQ=255)
                                                                                                                 aCGTTTAACAGCTTCCCCCCCACGTCATTTTCTCAACACGTTGGCCTGGTGAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTcc                                    <  2:1657155/150‑1 (MQ=255)
                                                                                                                 aCGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTcc                                    >  2:4310732/1‑150 (MQ=255)
                                                                                                                  cGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCa                                   <  2:1937885/150‑1 (MQ=255)
                                                                                                                   gTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCAt                                  <  1:2470048/150‑1 (MQ=255)
                                                                                                                              tCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTtg                       <  1:2164834/150‑1 (MQ=255)
                                                                                                                               cccccACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTtgt                      <  1:2536363/150‑1 (MQ=255)
                                                                                                                               cccccACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTtgt                      <  2:3137491/150‑1 (MQ=255)
                                                                                                                                ccccACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGACTTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTtgtg                     <  2:170346/150‑1 (MQ=255)
                                                                                                                                ccccACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTtgtg                     <  2:2960492/150‑1 (MQ=255)
                                                                                                                                  ccACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCt                   >  1:169781/1‑150 (MQ=255)
                                                                                                                                     cacGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTAAATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTggg                >  1:1973630/1‑150 (MQ=255)
                                                                                                                                      acGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGggc               >  1:2376171/1‑150 (MQ=255)
                                                                                                                                       cGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGTATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGggcg              >  2:3821941/1‑150 (MQ=255)
                                                                                                                                        gTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCGTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGggcgg             <  2:3673364/150‑1 (MQ=255)
                                                                                                                                           aTTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGgcgc          >  1:2693957/1‑150 (MQ=255)
                                                                                                                                           aTTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGgcgc          >  1:53300/1‑150 (MQ=255)
                                                                                                                                            ttttCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCt         <  2:4336956/150‑1 (MQ=255)
                                                                                                                                                 tcAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCt    <  1:3909892/150‑1 (MQ=255)
                                                                                                                                                 tcAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCt    <  2:274140/150‑1 (MQ=255)
                                                                                                                                                  cAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTa   >  1:2911339/1‑150 (MQ=255)
                                                                                                                                                  cAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTa   >  1:3240893/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TCGAGTCAATGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAG  >  CP009273/3025531‑3025827

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