Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,025,678 A→T 8.8% 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,6780AT 8.8% 905.9 / 53.7 288Q318L (CAA→CTA) uacTuric acid permease
Reads supporting (aligned to +/- strand):  ref base A (133/129);  new base T (15/11);  total (148/140)
Fisher's exact test for biased strand distribution p-value = 5.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 288 total aligned reads displayed.

AGTCAATGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTA  >  CP009273/3025534‑3025826
                                                                                                                                                |                                                                                                                                                    
 gTCAATGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGt                                                                                                                                                >  1:1388858/1‑150 (MQ=255)
  tCAATGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGtt                                                                                                                                               <  1:3098883/150‑1 (MQ=255)
  tCAATGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGtt                                                                                                                                               >  2:2694635/1‑150 (MQ=255)
   cAATGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTg                                                                                                                                              <  1:2969642/150‑1 (MQ=255)
   cAATGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTg                                                                                                                                              >  2:4214500/1‑150 (MQ=255)
           aTGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGtt                                                                                                                                      >  1:4200277/1‑150 (MQ=255)
           aTGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGAGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAAAGTTGGCCTGGtt                                                                                                                                      >  1:4200983/1‑150 (MQ=255)
            tGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTa                                                                                                                                     <  1:2892229/150‑1 (MQ=255)
              ttCCTGGCACTGTGTGATATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTTTCGATGGCGTATGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGc                                                                                                                                   <  1:2055715/150‑1 (MQ=255)
               tCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCACAACGTTGGCCTGGTTAGCg                                                                                                                                  >  2:1361010/1‑150 (MQ=255)
               tCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCg                                                                                                                                  <  2:3114646/150‑1 (MQ=255)
                 cTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTg                                                                                                                                >  1:760278/1‑150 (MQ=255)
                 cTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCAATTTCTCAAAACGTTGGCCAGGTTAGCGTg                                                                                                                                >  1:2638659/1‑150 (MQ=255)
                   ggCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGAc                                                                                                                              >  1:3730555/1‑150 (MQ=255)
                    gCACTGGGTGAAATAGTCGGTCGTAAACTCTCGTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGAcg                                                                                                                             <  1:2172176/150‑1 (MQ=255)
                          ggTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGtt                                                                                                                       >  1:2620873/1‑150 (MQ=255)
                           gTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTc                                                                                                                      <  1:876644/150‑1 (MQ=255)
                           gTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTc                                                                                                                      <  1:3527599/150‑1 (MQ=255)
                           gTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTc                                                                                                                      <  1:1967871/150‑1 (MQ=255)
                            tGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCa                                                                                                                     <  1:564907/150‑1 (MQ=255)
                            tGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCa                                                                                                                     <  2:1086743/150‑1 (MQ=255)
                             gAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCAt                                                                                                                    <  2:280320/150‑1 (MQ=255)
                              aaaTAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATa                                                                                                                   >  1:707363/1‑150 (MQ=255)
                              aaaTAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATa                                                                                                                   <  2:365768/150‑1 (MQ=255)
                                 tAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGcc                                                                                                                <  1:1257055/150‑1 (MQ=255)
                                 tAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGcc                                                                                                                <  1:351529/150‑1 (MQ=255)
                                    tcggtcgTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGcc                                                                                                             <  2:2082389/150‑2 (MQ=255)
                                    tcggtcgTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCt                                                                                                             <  2:2443999/150‑1 (MQ=255)
                                     cggtcgTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTg                                                                                                            <  2:4209591/150‑1 (MQ=255)
                                     cggtcgTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTg                                                                                                            <  1:303936/150‑1 (MQ=255)
                                       gtcgTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTggg                                                                                                          >  2:1236041/1‑150 (MQ=255)
                                        tcgTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGgt                                                                                                         >  2:1633975/1‑150 (MQ=255)
                                          gTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGgtgt                                                                                                       <  1:4117903/150‑1 (MQ=255)
                                          gTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGgtgt                                                                                                       <  2:2864170/150‑1 (MQ=255)
                                          gTAAACACTATTCGCACGAAATAATTCGCGGGCTGCGTGTCGATGGCGAAGGGACAAAGATAGGAGGCACGTTTAAAAGCTTCCACCACACGACATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCAGCTGGgtgt                                                                                                       >  2:1322059/1‑150 (MQ=255)
                                           tAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGgtgtg                                                                                                      >  2:4103122/1‑150 (MQ=255)
                                             aaCTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTa                                                                                                    <  2:4025966/150‑1 (MQ=255)
                                              aCTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTAt                                                                                                   <  2:2149519/150‑1 (MQ=255)
                                                  tctTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATttct                                                                                               >  2:2613772/1‑150 (MQ=255)
                                                    ttCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATttcttc                                                                                             <  1:1088489/150‑1 (MQ=255)
                                                    ttCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATttcttc                                                                                             >  2:3036423/1‑150 (MQ=255)
                                                    ttCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATttcttc                                                                                             <  2:1103440/150‑1 (MQ=255)
                                                    ttCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATttcttc                                                                                             >  2:2307950/1‑150 (MQ=255)
                                                    ttCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATttcttc                                                                                             <  1:1715084/150‑1 (MQ=255)
                                                     tCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCg                                                                                            <  2:2257042/150‑1 (MQ=255)
                                                      cGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCgg                                                                                           <  1:4293405/150‑1 (MQ=255)
                                                      cGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCgg                                                                                           >  2:3790035/1‑150 (MQ=255)
                                                        cACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGa                                                                                         <  1:1819068/150‑1 (MQ=255)
                                                        cACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGa                                                                                         <  1:3477107/150‑1 (MQ=255)
                                                           gATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAAtt                                                                                      <  2:4076042/150‑1 (MQ=255)
                                                            atatTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATta                                                                                     >  1:847234/1‑150 (MQ=255)
                                                             tattatTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTAGGGAATtat                                                                                    >  2:4162116/1‑150 (MQ=255)
                                                                 attCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTa                                                                                >  1:3650991/1‑150 (MQ=255)
                                                                 attCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTa                                                                                >  1:451340/1‑150 (MQ=255)
                                                                    cgcgGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATc                                                                             >  1:462136/1‑150 (MQ=255)
                                                                        ggCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGt                                                                         <  2:3653942/150‑1 (MQ=255)
                                                                          cTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGtta                                                                       >  1:2749178/1‑149 (MQ=255)
                                                                          cTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGttt                                                                       >  1:3808308/1‑150 (MQ=255)
                                                                             cGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGGCGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGc                                                                    >  2:2980188/1‑150 (MQ=255)
                                                                              gtgtCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCa                                                                   >  1:4248825/1‑150 (MQ=255)
                                                                                gtCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATg                                                                 >  2:176527/1‑150 (MQ=255)
                                                                                gtCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATg                                                                 <  2:4510123/150‑1 (MQ=255)
                                                                                 tCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATgg                                                                >  1:163168/1‑150 (MQ=255)
                                                                                  cGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGt                                                               <  2:4529534/150‑1 (MQ=255)
                                                                                        cGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCaaa                                                         <  1:2993554/150‑1 (MQ=255)
                                                                                         gTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCaaaa                                                        <  2:1388858/150‑1 (MQ=255)
                                                                                         gTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCaaaa                                                        >  2:4246808/1‑150 (MQ=255)
                                                                                            gggACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATg                                                     >  1:2280301/1‑150 (MQ=255)
                                                                                              gACAATGATAGGCGGCACGATTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGc                                                   <  2:3257081/150‑1 (MQ=255)
                                                                                                cAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATCGTGCCAAAAATGGCgg                                                 >  2:1448112/1‑150 (MQ=255)
                                                                                                 aaTGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGt                                                <  1:2757511/150‑1 (MQ=255)
                                                                                                   tGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGc                                              >  1:3449607/1‑150 (MQ=255)
                                                                                                     agaGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTg                                            <  2:2803548/148‑1 (MQ=255)
                                                                                                     aTAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTg                                            >  2:2850971/1‑150 (MQ=255)
                                                                                                      tAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTgg                                           >  1:808616/1‑150 (MQ=255)
                                                                                                        ggcggcACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTa                                         <  1:2087377/150‑1 (MQ=255)
                                                                                                        ggcggcACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTa                                         <  1:323584/150‑1 (MQ=255)
                                                                                                        ggcggcACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTa                                         <  2:3143478/150‑1 (MQ=255)
                                                                                                            gcACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCt                                     <  1:4382030/150‑1 (MQ=255)
                                                                                                             cACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCAGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGACAAAAATGGCGGTGCTGGTAGCCTc                                    >  2:623971/1‑150 (MQ=255)
                                                                                                                gTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCAt                                 >  1:1217772/1‑150 (MQ=255)
                                                                                                                     aCAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGc                            >  2:3789596/1‑150 (MQ=255)
                                                                                                                     aCAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGc                            >  1:396315/1‑150 (MQ=255)
                                                                                                                        gCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAAt                         <  1:1540518/150‑1 (MQ=255)
                                                                                                                        gCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAAt                         <  1:2757059/150‑1 (MQ=255)
                                                                                                                        gCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAAt                         <  1:3601189/150‑1 (MQ=255)
                                                                                                                        gCTGCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAAt                         <  2:3353108/150‑1 (MQ=255)
                                                                                                                         cTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAAtt                        >  2:898723/1‑150 (MQ=255)
                                                                                                                          ttCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAAttt                       <  1:3189384/150‑1 (MQ=255)
                                                                                                                            cccccACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTtgt                     >  1:1830330/1‑150 (MQ=255)
                                                                                                                                     gTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGggcgg            >  1:639467/1‑150 (MQ=255)
                                                                                                                                     gTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGggcgg            >  1:2488026/1‑150 (MQ=255)
                                                                                                                                     gTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGggcgg            >  2:4424756/1‑150 (MQ=255)
                                                                                                                                     gTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGggcgg            <  2:1218091/150‑1 (MQ=255)
                                                                                                                                      tCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGggcggc           >  1:3843814/1‑150 (MQ=255)
                                                                                                                                          tttCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTg       >  1:1885641/1‑150 (MQ=255)
                                                                                                                                          tttCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGATGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTg       <  1:2807737/150‑1 (MQ=255)
                                                                                                                                             ctcAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTc    >  1:2855335/1‑150 (MQ=255)
                                                                                                                                               cAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTa  >  1:945574/1‑150 (MQ=255)
                                                                                                                                               cAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTa  <  2:3614495/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
AGTCAATGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTA  >  CP009273/3025534‑3025826

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