Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,025,678 A→T 20.2% 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,6780AT20.2% 672.9 / 158.9 270Q318L (CAA→CTA) uacTuric acid permease
Reads supporting (aligned to +/- strand):  ref base A (115/100);  new base T (29/26);  total (144/126)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 281 total aligned reads displayed.

TATCGAGTCAATGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAG  >  CP009273/3025529‑3025827
                                                                                                                                                     |                                                                                                                                                     
tATCGAGTCAATGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCa                                                                                                                                                       >  1:1221724/1‑150 (MQ=255)
tATCGAGTCAATGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCa                                                                                                                                                       <  1:3993269/150‑1 (MQ=255)
tATCGAGTCAATGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCa                                                                                                                                                       <  1:1804414/150‑1 (MQ=255)
   cGAGTCAATGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCaaaa                                                                                                                                                    >  2:2546915/1‑150 (MQ=255)
    gAGTCAATGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAAc                                                                                                                                                   >  1:1738826/1‑150 (MQ=255)
      gTCAATGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGt                                                                                                                                                 >  2:1269573/1‑150 (MQ=255)
      gTCAATGGGGATGTTCATGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCACCACACGTCATATTCTCAAAACGt                                                                                                                                                 >  1:1483170/1‑150 (MQ=255)
       tCAATGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGtt                                                                                                                                                <  2:2924813/150‑1 (MQ=255)
        cAATGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTg                                                                                                                                               <  1:3011199/150‑1 (MQ=255)
         aaTGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTgg                                                                                                                                              >  2:1679318/1‑150 (MQ=255)
         aaTGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTgg                                                                                                                                              <  2:1854751/150‑1 (MQ=255)
           tGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGcc                                                                                                                                            <  1:345679/150‑1 (MQ=255)
               gATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGt                                                                                                                                        >  2:4207850/1‑150 (MQ=255)
               gATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGt                                                                                                                                        >  2:3068157/1‑150 (MQ=255)
                aTGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGtt                                                                                                                                       >  1:1149277/1‑150 (MQ=255)
                    tCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCg                                                                                                                                   <  2:422620/150‑1 (MQ=255)
                     ccTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTAGTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGt                                                                                                                                  <  1:498828/150‑1 (MQ=255)
                       tGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGa                                                                                                                                >  1:1448493/1‑150 (MQ=255)
                       tGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGa                                                                                                                                <  2:1661330/150‑1 (MQ=255)
                          cACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGAcgc                                                                                                                             >  2:808152/1‑150 (MQ=255)
                               ggTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGtt                                                                                                                        >  1:3779924/1‑150 (MQ=255)
                                gTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTc                                                                                                                       <  1:2097914/150‑1 (MQ=255)
                                 tGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCa                                                                                                                      >  2:4080055/1‑150 (MQ=255)
                                    aaTAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAg                                                                                                                   >  1:385706/1‑150 (MQ=255)
                                    aaTAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCGGGTTAGCGTGACGCGCGTTCATAg                                                                                                                   >  1:2142952/1‑150 (MQ=255)
                                     aTAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACCTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGc                                                                                                                  >  1:839777/1‑150 (MQ=255)
                                     aTAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGGTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGc                                                                                                                  >  2:904354/1‑150 (MQ=255)
                                        gtcggtcgTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGc                                                                                                               >  1:542505/1‑150 (MQ=255)
                                          cggtcgTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTg                                                                                                             <  1:2897893/150‑1 (MQ=255)
                                          cggtcgTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTg                                                                                                             <  2:3711557/150‑1 (MQ=255)
                                           ggtcgTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTgg                                                                                                            <  2:1068695/150‑1 (MQ=255)
                                              cgTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGgtg                                                                                                         >  1:2452729/1‑150 (MQ=255)
                                              cgTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGgtg                                                                                                         <  1:79138/150‑1 (MQ=255)
                                                tAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGgtgtg                                                                                                       <  1:3239132/150‑1 (MQ=255)
                                                tAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGgtgtg                                                                                                       <  2:3099488/150‑1 (MQ=255)
                                                tAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGgtgtg                                                                                                       >  1:2753682/1‑150 (MQ=255)
                                                    ctctctTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTAtt                                                                                                   <  2:154208/150‑1 (MQ=255)
                                                     tctctTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTAttt                                                                                                  >  1:1574168/1‑150 (MQ=255)
                                                      ctctTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATttc                                                                                                 >  2:2584606/1‑150 (MQ=255)
                                                       tctTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATttct                                                                                                <  1:1324691/150‑1 (MQ=255)
                                                       tctTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATttct                                                                                                <  2:1652305/150‑1 (MQ=255)
                                                          tCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCg                                                                                             <  2:2019644/150‑1 (MQ=255)
                                                           cGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCgg                                                                                            >  2:378011/1‑150 (MQ=255)
                                                           cGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCgg                                                                                            >  2:376819/1‑150 (MQ=255)
                                                           cGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCgg                                                                                            <  2:1179136/150‑1 (MQ=255)
                                                            gCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCggg                                                                                           >  2:3879801/1‑150 (MQ=255)
                                                                gATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAAtt                                                                                       >  2:2299702/1‑150 (MQ=255)
                                                                 atatTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATta                                                                                      <  1:1456749/150‑1 (MQ=255)
                                                                 atatTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATta                                                                                      >  2:1387922/1‑150 (MQ=255)
                                                                 atatTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATta                                                                                      >  2:666457/1‑150 (MQ=255)
                                                                   attattCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATtata                                                                                    <  1:1848043/150‑1 (MQ=255)
                                                                    ttattCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCACAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGACTTatat                                                                                   >  1:2693201/1‑150 (MQ=255)
                                                                     tattCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATAtt                                                                                  >  1:4177804/1‑150 (MQ=255)
                                                                           cgGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCt                                                                            <  2:1763502/150‑1 (MQ=255)
                                                                               cTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGttt                                                                        >  1:630816/1‑150 (MQ=255)
                                                                               cTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGttt                                                                        <  1:2819295/150‑1 (MQ=255)
                                                                               cTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGttt                                                                        >  1:2615294/1‑150 (MQ=255)
                                                                                tGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCTCGTTTAAGAGCGTCCCCCACACGTCATTTTCTCAAATCGTGGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTTATCCTGTTTg                                                                       <  2:4349709/150‑1 (MQ=255)
                                                                                tGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTg                                                                       >  2:1294510/1‑150 (MQ=255)
                                                                                 gCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTgg                                                                      <  2:595037/150‑1 (MQ=255)
                                                                                  cGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCGTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGc                                                                     <  2:4068336/150‑1 (MQ=255)
                                                                                   gtgtCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCa                                                                    <  2:2557891/150‑1 (MQ=255)
                                                                                        gATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTg                                                               >  1:902444/1‑150 (MQ=255)
                                                                                                aGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAAt                                                       <  1:3133907/150‑1 (MQ=255)
                                                                                                   gACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGc                                                    >  1:4034596/1‑150 (MQ=255)
                                                                                                     cAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCgg                                                  <  2:3153615/150‑1 (MQ=255)
                                                                                                      aaTGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGt                                                 >  1:3305217/1‑150 (MQ=255)
                                                                                                      aaTGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGt                                                 <  1:2289465/150‑1 (MQ=255)
                                                                                                               cggcACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGc                                        >  1:4309976/1‑150 (MQ=255)
                                                                                                                ggcACGTTTAACAGCTTCCCCCACACGTCATTGTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGcc                                       >  1:3955381/1‑150 (MQ=255)
                                                                                                                   aCGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCGTGGTGCCAAAAATGGCGGTGCTGGTAGCCTcc                                    >  2:1241018/1‑150 (MQ=255)
                                                                                                                   aCGTTTAACAGCTTCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTcc                                    >  1:4104492/1‑150 (MQ=255)
                                                                                                                    cGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCa                                   <  2:445324/150‑1 (MQ=255)
                                                                                                                     gTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTCTTTCTTCGGGACTTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCAt                                  >  2:3487131/1‑150 (MQ=255)
                                                                                                                      tttAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCAtt                                 >  1:1529853/1‑150 (MQ=255)
                                                                                                                      tttAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCAtt                                 >  2:2901475/1‑150 (MQ=255)
                                                                                                                        tAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTcc                               >  1:3533765/1‑150 (MQ=255)
                                                                                                                          aCAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGc                             >  2:4008294/1‑150 (MQ=255)
                                                                                                                           cAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCa                            <  1:1039960/150‑1 (MQ=255)
                                                                                                                             gCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAAt                          <  1:3617640/150‑1 (MQ=255)
                                                                                                                              cTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAAtt                         >  1:66116/1‑150 (MQ=255)
                                                                                                                                tCCCCCACACGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCC‑AAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTtgt                      >  2:114297/1‑150 (MQ=255)
                                                                                                                                  ccccACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCAGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTtgtg                     <  1:1269956/150‑1 (MQ=255)
                                                                                                                                      acacGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTgg                 <  2:3626163/150‑1 (MQ=255)
                                                                                                                                      acacGTCATTTTCTCTAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTgg                 <  2:2828332/150‑1 (MQ=255)
                                                                                                                                       cacGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTggg                >  2:2590096/1‑150 (MQ=255)
                                                                                                                                       cacGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTggg                <  1:971223/150‑1 (MQ=255)
                                                                                                                                       cacGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCACAATTTGTGCTggg                >  2:1047935/1‑150 (MQ=255)
                                                                                                                                       cacGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTATATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTggg                <  2:2520001/150‑1 (MQ=255)
                                                                                                                                            cATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGgcggcg           <  1:92434/150‑1 (MQ=255)
                                                                                                                                            cATTTTCTCAAAACGTTGGCCTGGTAAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGgcggcg           >  1:216228/1‑150 (MQ=255)
                                                                                                                                            cATTTTCTCAAAACGTTCGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTCGCCTCCATTCCGCAGTTTGTGCTGGgcggcg           >  2:368895/1‑150 (MQ=255)
                                                                                                                                             aTTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGgcgc          >  2:2659560/1‑150 (MQ=255)
                                                                                                                                             aTTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGgcgc          >  2:487346/1‑150 (MQ=255)
                                                                                                                                                 tctcAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGt      <  1:1963713/150‑1 (MQ=255)
                                                                                                                                                 tctcAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGt      >  1:3814061/1‑150 (MQ=255)
                                                                                                                                                  ctcAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTc     <  2:1483170/150‑1 (MQ=255)
                                                                                                                                                  ctcAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTc     >  2:3556722/1‑150 (MQ=255)
                                                                                                                                                   tcAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCt    <  2:3445361/150‑1 (MQ=255)
                                                                                                                                                     aaaaCGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAg  <  2:4310466/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TATCGAGTCAATGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAG  >  CP009273/3025529‑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: ‑