Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,025,678 A→T 41.7% 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,6780AT41.7% 146.0 / inf 230Q318L (CAA→CTA) uacTuric acid permease
Reads supporting (aligned to +/- strand):  ref base A (59/75);  new base T (40/55);  total (100/130)
Fisher's exact test for biased strand distribution p-value = 7.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.37e-01

Only 100 of 235 total aligned reads displayed.

TATCGAGTCAATGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTA  >  CP009273/3025529‑3025826
                                                                                                                                                     |                                                                                                                                                    
   cGAGTCAATGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCtaaa                                                                                                                                                   <  1:3132383/150‑1 (MQ=255)
    gAGTCAATGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAAc                                                                                                                                                  <  1:103269/150‑1 (MQ=255)
    gAGTCAATGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAAc                                                                                                                                                  >  1:3328428/1‑150 (MQ=255)
    gAGTCAATGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAAc                                                                                                                                                  >  1:233/1‑150 (MQ=255)
         aaTGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATATTCTCAAAACGTTgg                                                                                                                                             >  2:2647307/1‑150 (MQ=255)
          aTGGGGATGTTCCTGGCTCTGGGTGAAATAGTCGGTCGTTATCTCTCTTCGCTCGATATTAGTCGCGGGCTGCGTGTGGATGGCGTAGGGACAATGTTAGGCGTCAAGTTTTACAGCTTCCCCCACAAGTCATTTTCTCTAAACGTTGGc                                                                                                                                            <  1:1291199/150‑1 (MQ=255)
            ggggATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCt                                                                                                                                          <  2:3617273/150‑1 (MQ=255)
            ggggATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCt                                                                                                                                          <  2:3617234/150‑1 (MQ=255)
             gggATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTg                                                                                                                                         <  2:678422/150‑1 (MQ=255)
              ggATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTgg                                                                                                                                        <  2:3486549/150‑1 (MQ=255)
                 tGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTa                                                                                                                                     <  2:3614706/150‑1 (MQ=255)
                 tGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTa                                                                                                                                     <  2:4005216/150‑1 (MQ=255)
                      cTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTg                                                                                                                                <  2:2698300/150‑1 (MQ=255)
                       tGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGa                                                                                                                               >  1:3242782/1‑150 (MQ=255)
                       tGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGAACGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGa                                                                                                                               >  1:2287841/1‑150 (MQ=255)
                        ggCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGAc                                                                                                                              >  1:2693327/1‑150 (MQ=255)
                             tGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGAcgcgcg                                                                                                                         >  1:3565523/1‑150 (MQ=255)
                                gTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTc                                                                                                                      >  2:250908/1‑150 (MQ=255)
                                gTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTc                                                                                                                      <  1:375293/150‑1 (MQ=255)
                                     aTAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGc                                                                                                                 <  1:490632/150‑1 (MQ=255)
                                      tAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGcc                                                                                                                <  1:604339/150‑1 (MQ=255)
                                       aGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCg                                                                                                               <  1:2920614/150‑1 (MQ=255)
                                       aGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCg                                                                                                               >  1:2046994/1‑150 (MQ=255)
                                        gtcggtcgTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGc                                                                                                              >  1:126262/1‑150 (MQ=255)
                                         tcggtcgTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCt                                                                                                             >  2:4145652/1‑150 (MQ=255)
                                           ggtcgTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTgg                                                                                                           <  2:2020661/150‑1 (MQ=255)
                                                 aaaCTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGgggtgt                                                                                                     >  2:1520532/1‑150 (MQ=255)
                                                  aaCTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTa                                                                                                    >  1:92646/1‑150 (MQ=255)
                                                  aaCTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTa                                                                                                    >  1:717427/1‑150 (MQ=255)
                                                   aCTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTAt                                                                                                   >  2:445256/1‑150 (MQ=255)
                                                     tctctTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTAttt                                                                                                 <  2:317102/150‑1 (MQ=255)
                                                      ctctTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATttc                                                                                                <  2:1809902/150‑1 (MQ=255)
                                                       tctTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATttct                                                                                               <  2:2155324/150‑1 (MQ=255)
                                                       tctTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATttct                                                                                               <  2:2582670/150‑1 (MQ=255)
                                                        ctTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATttctt                                                                                              >  2:1062460/1‑150 (MQ=255)
                                                           cGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCgg                                                                                           >  2:3927800/1‑150 (MQ=255)
                                                           cGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCgg                                                                                           >  1:1532733/1‑150 (MQ=255)
                                                            gCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCggg                                                                                          >  1:627538/1‑150 (MQ=255)
                                                            gCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCggg                                                                                          <  1:460183/150‑1 (MQ=255)
                                                            gCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCggg                                                                                          <  1:1115503/150‑1 (MQ=255)
                                                              aCGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCGTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGaa                                                                                        <  1:2502263/150‑1 (MQ=255)
                                                               cGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAAt                                                                                       <  1:869356/150‑1 (MQ=255)
                                                                gATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAAtt                                                                                      >  1:3897336/1‑150 (MQ=255)
                                                                 atatTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATta                                                                                     <  1:1701250/150‑1 (MQ=255)
                                                                 atatTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATta                                                                                     <  1:3799503/150‑1 (MQ=255)
                                                                   attattCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATtata                                                                                   <  1:610708/150‑1 (MQ=255)
                                                                   attattCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATtata                                                                                   <  2:2816987/150‑1 (MQ=255)
                                                                    ttagtCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTatat                                                                                  <  2:1694916/150‑1 (MQ=255)
                                                                      attCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTa                                                                                <  2:3513851/150‑1 (MQ=255)
                                                                         cgcgGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATc                                                                             <  2:2387738/150‑1 (MQ=255)
                                                                           cgGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCt                                                                           >  1:2901285/1‑150 (MQ=255)
                                                                           cgGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCt                                                                           <  2:2453514/150‑1 (MQ=255)
                                                                            gggCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTg                                                                          >  1:2975047/1‑150 (MQ=255)
                                                                            gggCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATGTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTg                                                                          >  2:862193/1‑150 (MQ=255)
                                                                              gCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGtt                                                                        <  2:3052171/150‑1 (MQ=255)
                                                                              gCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGtt                                                                        <  2:708890/150‑1 (MQ=255)
                                                                                tGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTg                                                                      >  1:3691432/1‑150 (MQ=255)
                                                                                   gtgtCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATGTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCa                                                                   >  2:34922/1‑150 (MQ=255)
                                                                                      tCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATgg                                                                >  1:1971495/1‑150 (MQ=255)
                                                                                        gATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTg                                                              <  1:933993/150‑1 (MQ=255)
                                                                                        gATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTg                                                              >  1:1416610/1‑150 (MQ=255)
                                                                                           ggCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCa                                                           >  2:778734/1‑150 (MQ=255)
                                                                                           ggCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCa                                                           <  1:2525069/150‑1 (MQ=255)
                                                                                           ggCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCa                                                           >  2:1100581/1‑150 (MQ=255)
                                                                                             cGTAGGGACAATGATAGGCGGCACGTTTAACCGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCaaa                                                         <  1:898523/150‑1 (MQ=255)
                                                                                             cGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCaaa                                                         <  2:3241595/150‑1 (MQ=255)
                                                                                              gTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCaaaa                                                        >  1:2463044/1‑150 (MQ=255)
                                                                                                aGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAAt                                                      <  1:80633/150‑1 (MQ=255)
                                                                                                    aCAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCg                                                  <  2:2039054/150‑1 (MQ=255)
                                                                                                       aTGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTg                                               <  2:3426627/150‑1 (MQ=255)
                                                                                                        tGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGc                                              <  2:3001903/150‑1 (MQ=255)
                                                                                                         gATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCt                                             <  1:1579168/150‑1 (MQ=255)
                                                                                                          aTAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTCTATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTg                                            >  2:3138890/1‑150 (MQ=255)
                                                                                                          aTAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTg                                            >  1:613976/1‑150 (MQ=255)
                                                                                                          aTAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTg                                            >  2:1048986/1‑150 (MQ=255)
                                                                                                           tAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTgg                                           >  2:758624/1‑150 (MQ=255)
                                                                                                           tAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTgg                                           <  2:2271240/150‑1 (MQ=255)
                                                                                                            aGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGt                                          <  1:568302/150‑1 (MQ=255)
                                                                                                            aGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGt                                          <  2:3358642/150‑1 (MQ=255)
                                                                                                               cggcACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGc                                       >  1:3211222/1‑150 (MQ=255)
                                                                                                                ggcACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGcc                                      <  1:1606964/150‑1 (MQ=255)
                                                                                                                 gcACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCt                                     <  2:1620327/150‑1 (MQ=255)
                                                                                                                 gcACGTTTAACAGCGTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCt                                     <  1:272383/150‑1 (MQ=255)
                                                                                                                   aCGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTcc                                   <  2:3326206/150‑1 (MQ=255)
                                                                                                                        tAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTcc                              <  1:3991897/150‑1 (MQ=255)
                                                                                                                          aCAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGc                            >  2:1270525/1‑150 (MQ=255)
                                                                                                                                      acacGTCATTTTCTCGAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTgg                >  1:4085521/1‑150 (MQ=255)
                                                                                                                                      acacGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGACTTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTgg                >  2:56806/1‑150 (MQ=255)
                                                                                                                                       cacGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTggg               <  1:285235/150‑1 (MQ=255)
                                                                                                                                            cATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGgcggcg          >  1:2938812/1‑150 (MQ=255)
                                                                                                                                             aTTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGgcgc         >  1:3504707/1‑150 (MQ=255)
                                                                                                                                              ttttCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGTGCt        >  1:775530/1‑150 (MQ=255)
                                                                                                                                               tttCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTg       >  2:3902038/1‑150 (MQ=255)
                                                                                                                                               tttCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTg       >  1:1851359/1‑150 (MQ=255)
                                                                                                                                                ttCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTgg      <  1:440699/150‑1 (MQ=255)
                                                                                                                                                 tctcTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGt     <  2:1921993/150‑1 (MQ=255)
                                                                                                                                                 tctcTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGt     <  2:2836576/150‑1 (MQ=255)
                                                                                                                                                  ctcAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTc    >  1:2951614/1‑150 (MQ=255)
                                                                                                                                                   tcAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCt   <  1:1634311/150‑1 (MQ=255)
                                                                                                                                                   tcAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCt   <  1:700279/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TATCGAGTCAATGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTA  >  CP009273/3025529‑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: