Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,815,043 C→A 16.9% R86S (CGT→AGT)  gmk → guanylate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,815,0430CA16.9% 977.4 / 118.9 312R86S (CGT→AGT) gmkguanylate kinase
Reads supporting (aligned to +/- strand):  ref base C (152/107);  new base A (33/20);  total (185/127)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 319 total aligned reads displayed.

GTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGACCGCCGTCTA  >  CP009273/3814896‑3815192
                                                                                                                                                   |                                                                                                                                                     
gTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGt                                                                                                                                                     >  1:301600/1‑150 (MQ=255)
gTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGAGt                                                                                                                                                     >  2:3104867/1‑150 (MQ=255)
 tttCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGAGTg                                                                                                                                                    <  2:3586910/150‑1 (MQ=255)
      cacaccacacGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGcc                                                                                                                                               >  1:145847/1‑150 (MQ=255)
          ccacacGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTg                                                                                                                                           <  1:2425441/150‑1 (MQ=255)
            acacGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGAGTGAGGCCATTGAg                                                                                                                                         >  2:4536934/1‑150 (MQ=255)
               cGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCaa                                                                                                                                      >  1:1645371/1‑150 (MQ=255)
                  cAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTa                                                                                                                                   >  2:2815786/1‑150 (MQ=255)
                    aCCGCGTCCTGGTGATTTCCTCGGTGAACTTTTTTTCTTTGTTAATCATGATGAATTTAAAGAAATGTTTAGCTGAGATGCGTTCCTCGTACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGATCAAGTACt                                                                                                                                 <  1:1408778/150‑1 (MQ=255)
                       gcgTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGc                                                                                                                              <  2:1447387/150‑1 (MQ=255)
                         gTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGa                                                                                                                            >  2:153324/1‑150 (MQ=255)
                           ccTGGTGAAGTCCACGGTGAACATTATTTATTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGAcc                                                                                                                          <  2:4308493/150‑1 (MQ=255)
                            cTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCg                                                                                                                         <  2:2683102/150‑1 (MQ=255)
                              ggTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGgt                                                                                                                       <  1:1879707/150‑1 (MQ=255)
                               gTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGgtg                                                                                                                      <  1:1733994/150‑1 (MQ=255)
                                  aaGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCg                                                                                                                   >  2:1857817/1‑150 (MQ=255)
                                  aaGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCg                                                                                                                   >  1:2657789/1‑150 (MQ=255)
                                     tCCACGGTGAACATTATTTCTTTGTTAATCCTGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATg                                                                                                                >  1:2828535/1‑150 (MQ=255)
                                      ccACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGt                                                                                                               >  2:771309/1‑150 (MQ=255)
                                         cGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGtttt                                                                                                            >  1:513977/1‑150 (MQ=255)
                                          ggTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGttttt                                                                                                           >  1:2276234/1‑150 (MQ=255)
                                           gTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTtc                                                                                                          >  1:3338651/1‑150 (MQ=255)
                                           gTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGAGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTtc                                                                                                          <  2:3856890/150‑1 (MQ=255)
                                            tGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTtct                                                                                                         >  2:3004673/1‑150 (MQ=255)
                                              aaCATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCg                                                                                                       <  2:1566891/150‑1 (MQ=255)
                                                cATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGat                                                                                                     >  1:3741666/1‑150 (MQ=255)
                                                 attattTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGata                                                                                                    >  1:854446/1‑150 (MQ=255)
                                                  ttattTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGatat                                                                                                   >  1:3863369/1‑150 (MQ=255)
                                                     tttCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGa                                                                                                >  1:3416532/1‑150 (MQ=255)
                                                     tttCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGAGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGa                                                                                                >  2:736707/1‑150 (MQ=255)
                                                     tttCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACGATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGa                                                                                                <  1:4311163/150‑1 (MQ=255)
                                                       tCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACt                                                                                              >  1:2543168/1‑150 (MQ=255)
                                                       tCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACt                                                                                              <  1:4636259/150‑1 (MQ=255)
                                                       tCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACt                                                                                              <  2:4063262/150‑1 (MQ=255)
                                                        cTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTg                                                                                             >  2:3931480/1‑150 (MQ=255)
                                                             ttAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTNAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGNCTGGCAgg                                                                                        >  1:2890596/1‑150 (MQ=255)
                                                             ttAATCATGATGAATTTAAAGAAATGATAAGCAGAGATGCGTTACTCGAACACGCACAAGTTTTTCGTAATTACTATGGCACATCGAGTGAGGCCATTGAGCAAGTACTGGCTACCGGTGTCGATGTTTTTCTCGATATAGACTGGCAgg                                                                                        >  2:808797/1‑150 (MQ=255)
                                                               aaTCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGgc                                                                                      <  2:2919763/150‑1 (MQ=255)
                                                                 tCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGAGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGgcgc                                                                                    >  1:4349798/1‑150 (MQ=255)
                                                                  cATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGAGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGgcgcg                                                                                   >  1:421565/1‑150 (MQ=255)
                                                                    tgatgaATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCgca                                                                                 <  2:567526/150‑1 (MQ=255)
                                                                    tgatgaATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCgca                                                                                 >  2:4342975/1‑150 (MQ=255)
                                                                     gatgaATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACGATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCgcag                                                                                <  1:1421492/150‑1 (MQ=255)
                                                                      atgaATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCgcagc                                                                               >  2:1666451/1‑150 (MQ=255)
                                                                        gaATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCaa                                                                             >  1:1711716/1‑150 (MQ=255)
                                                                         aaTTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCaaa                                                                            >  1:1247875/1‑150 (MQ=255)
                                                                           tttAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAAtt                                                                          <  1:3920311/150‑1 (MQ=255)
                                                                           tttAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAAtt                                                                          <  1:3828337/150‑1 (MQ=255)
                                                                           tttAAAGAAATGATTAGCAGAGATGCGGTCCTCGAACACGCAGAAGTTTTTGGTAATTACGATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAAtt                                                                          <  2:2328218/150‑1 (MQ=255)
                                                                            ttAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTc                                                                         >  2:714806/1‑150 (MQ=255)
                                                                             tAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCg                                                                        <  1:4027569/150‑1 (MQ=255)
                                                                             tAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCg                                                                        <  1:964856/150‑1 (MQ=255)
                                                                             tAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCg                                                                        <  1:3514153/150‑1 (MQ=255)
                                                                                aGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCa                                                                     >  2:3321568/1‑150 (MQ=255)
                                                                                aGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGAGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCa                                                                     <  1:303952/150‑1 (MQ=255)
                                                                                 gAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCag                                                                    <  1:280830/150‑1 (MQ=255)
                                                                                     tGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCagaaga                                                                >  1:2192722/1‑150 (MQ=255)
                                                                                      gATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGAt                                                               <  1:4522220/150‑1 (MQ=255)
                                                                                       aTTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTTGTAATTACTATGGCACATCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATAGACTGGCAGGGCGCGCCGCAAATTCGCCAGAAGATg                                                              >  1:2719686/1‑150 (MQ=255)
                                                                                       aTTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATg                                                              >  2:1094888/1‑150 (MQ=255)
                                                                                       aTTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACGTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATg                                                              <  2:4498839/150‑1 (MQ=255)
                                                                                          aGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCg                                                           <  2:3152949/150‑1 (MQ=255)
                                                                                                gaTGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCAcgcg                                                     >  2:440375/1‑150 (MQ=255)
                                                                                                   gCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCgg                                                  <  2:2060879/150‑1 (MQ=255)
                                                                                                    cGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGAGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCCCGCGCGGa                                                 >  1:1178337/1‑150 (MQ=255)
                                                                                                          tCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGAGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCt                                           >  1:2896125/1‑150 (MQ=255)
                                                                                                            gAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCttt                                         >  2:1758566/1‑150 (MQ=255)
                                                                                                            gAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGAGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCttt                                         <  1:2582556/150‑1 (MQ=255)
                                                                                                             aaCACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTa                                        >  1:996011/1‑150 (MQ=255)
                                                                                                             aaCACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTa                                        <  2:3492843/150‑1 (MQ=255)
                                                                                                               cacGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTAtt                                      <  2:1414875/150‑1 (MQ=255)
                                                                                                                 cGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTAtttt                                    <  1:3993546/150‑1 (MQ=255)
                                                                                                                 cGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTAtttt                                    <  1:3973295/150‑1 (MQ=255)
                                                                                                                 cGCAGAAGTTTTTGGTAATTACTATGGCACTTCGAGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTAtttt                                    <  1:3149058/150‑1 (MQ=255)
                                                                                                                  gCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTa                                   >  2:2447750/1‑150 (MQ=255)
                                                                                                                  gCAGAAGTTTTTGGTAATTACTATGGCACTTCGAGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTa                                   <  1:1939496/150‑1 (MQ=255)
                                                                                                                     gAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCCCGCGCGGAGTATCTTTATTTTAccg                                >  1:4401975/1‑150 (MQ=255)
                                                                                                                        gTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTAccgccg                             <  1:3785759/150‑1 (MQ=255)
                                                                                                                        gTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTAccgccg                             >  2:2228232/1‑150 (MQ=255)
                                                                                                                         tttttGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGt                            >  1:622986/1‑150 (MQ=255)
                                                                                                                         tttttGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGt                            <  2:1859327/150‑1 (MQ=255)
                                                                                                                             tGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCaa                        <  1:4130186/150‑1 (MQ=255)
                                                                                                                             tGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTA‑TTTACCGCCGTCCaaa                       >  2:1657482/1‑150 (MQ=255)
                                                                                                                                tAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAAt                     >  1:2055206/1‑150 (MQ=255)
                                                                                                                                  aTTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTg                   <  2:3658544/150‑1 (MQ=255)
                                                                                                                                  aTTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTg                   >  1:956922/1‑150 (MQ=255)
                                                                                                                                      cTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACt               >  1:589289/1‑150 (MQ=255)
                                                                                                                                      cTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACt               >  2:1739848/1‑150 (MQ=255)
                                                                                                                                        aTGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTgg             >  2:3040/1‑150 (MQ=255)
                                                                                                                                        aTGGCACTTCGAGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTgg             <  1:839581/150‑1 (MQ=255)
                                                                                                                                          ggCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGAc           <  2:854446/150‑1 (MQ=255)
                                                                                                                                           gCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGAcc          <  2:1373007/150‑1 (MQ=255)
                                                                                                                                            cACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGAccg         >  1:1681689/1‑150 (MQ=255)
                                                                                                                                            cACTTCGAGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGAccg         >  2:2127286/1‑150 (MQ=255)
                                                                                                                                            cACTTCGAGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGAccg         >  2:1969914/1‑150 (MQ=255)
                                                                                                                                             aCTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGAccgc        >  1:1019688/1‑150 (MQ=255)
                                                                                                                                                tCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGACCGCCGt     >  2:4165055/1‑150 (MQ=255)
                                                                                                                                                tCGAGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGACCGCCGt     >  1:749398/1‑150 (MQ=255)
                                                                                                                                                   cgTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGACCGCCGTCTa  >  1:829984/1‑150 (MQ=255)
                                                                                                                                                   cgTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGACCGCCGTCTa  <  1:4600081/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGACCGCCGTCTA  >  CP009273/3814896‑3815192

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