Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,268,967 A→T 11.7% S187S (TCA→TCT yjcD → inner membrane putative guanine permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,268,9670AT11.7% 794.7 / 81.3 264S187S (TCA→TCTyjcDinner membrane putative guanine permease
Reads supporting (aligned to +/- strand):  ref base A (118/115);  new base T (19/12);  total (137/127)
Fisher's exact test for biased strand distribution p-value = 3.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01

Only 100 of 265 total aligned reads displayed.

CTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCG  >  CP009273/4268819‑4269115
                                                                                                                                                    |                                                                                                                                                    
 tCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACt                                                                                                                                                    >  2:331897/1‑150 (MQ=255)
 tCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACt                                                                                                                                                    >  2:1387739/1‑150 (MQ=255)
  cACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCTCTg                                                                                                                                                   >  2:2656307/1‑150 (MQ=255)
    cGGTGTGGCGCACGGCACGGGGATTGTTATCGGTCTGTTCCTTCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCAGTTTCGCTGGGTGATTTCGCTACCTTCCCGGTGATTATGTCACTGGt                                                                                                                                                 <  2:4334174/150‑1 (MQ=255)
     ggTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTa                                                                                                                                                >  1:3169756/1‑150 (MQ=255)
      gtgtGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAg                                                                                                                                               <  1:3288523/150‑1 (MQ=255)
      gtgtGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAg                                                                                                                                               <  1:1938733/150‑1 (MQ=255)
      gtgtGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAg                                                                                                                                               <  2:1402599/150‑1 (MQ=255)
                cggcacggGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGt                                                                                                                                     >  2:3559136/1‑150 (MQ=255)
                 ggcacggGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGTTg                                                                                                                                    >  2:3985858/1‑150 (MQ=255)
                 ggcacggGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTg                                                                                                                                    <  1:2836025/150‑1 (MQ=255)
                      ggggATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGatcat                                                                                                                               >  2:2244449/1‑150 (MQ=255)
                      ggggATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGatcat                                                                                                                               >  2:2835404/1‑150 (MQ=255)
                      ggggATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGatcat                                                                                                                               >  1:2128308/1‑150 (MQ=255)
                      ggggATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGatcat                                                                                                                               <  2:1586276/150‑1 (MQ=255)
                      ggggATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGatcat                                                                                                                               >  2:4537531/1‑150 (MQ=255)
                      ggggATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATATCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGatcat                                                                                                                               >  2:2481128/1‑150 (MQ=255)
                         gATTGGTATCGGTCTGGTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGAGTAAAAACCCGCTGGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCgg                                                                                                                            <  1:4324442/150‑1 (MQ=255)
                          aTTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGc                                                                                                                           <  1:206691/150‑1 (MQ=255)
                          aTTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGc                                                                                                                           <  1:206743/150‑1 (MQ=255)
                          aTTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGc                                                                                                                           <  1:1100946/150‑1 (MQ=255)
                            tGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCt                                                                                                                         >  2:1769193/1‑150 (MQ=255)
                               tATCGGTCTGTTCCTGCTGCTCATTGCCGCTATCGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGa                                                                                                                      <  2:3895772/150‑1 (MQ=255)
                               tATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGa                                                                                                                      <  2:951336/150‑1 (MQ=255)
                                        gTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGaa                                                                                                             >  1:1704522/1‑150 (MQ=255)
                                         ttCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGaaa                                                                                                            >  1:3194428/1‑150 (MQ=255)
                                         ttCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGaaa                                                                                                            >  1:1793657/1‑150 (MQ=255)
                                          tCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAg                                                                                                           <  2:3512396/150‑1 (MQ=255)
                                                 gctCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCTCTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCtgg                                                                                                    <  2:3813670/150‑1 (MQ=255)
                                                 gctCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCtgg                                                                                                    <  2:1212587/150‑1 (MQ=255)
                                                   tCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCtggtg                                                                                                  <  2:2183851/150‑1 (MQ=255)
                                                   tCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCtggtg                                                                                                  <  2:728429/150‑1 (MQ=255)
                                                       tGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCAt                                                                                              <  1:2972304/150‑1 (MQ=255)
                                                        gCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCAtt                                                                                             <  2:1683193/150‑1 (MQ=255)
                                                         ccGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTc                                                                                            >  2:3259735/1‑150 (MQ=255)
                                                         ccGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTc                                                                                            >  2:4364314/1‑150 (MQ=255)
                                                          cGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTct                                                                                           <  1:3097191/150‑1 (MQ=255)
                                                           gCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTctg                                                                                          >  1:3565043/1‑150 (MQ=255)
                                                              aaCGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTctgctg                                                                                       <  1:1123571/150‑1 (MQ=255)
                                                                   tgtCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCTCTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCAt                                                                                  <  2:3363343/150‑1 (MQ=255)
                                                                   tgtCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCTCTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCAt                                                                                  <  2:1089672/150‑1 (MQ=255)
                                                                    gtCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCAtt                                                                                 >  1:3511840/1‑150 (MQ=255)
                                                                    gtCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCAtt                                                                                 >  1:1240530/1‑150 (MQ=255)
                                                                       ggTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATc                                                                              >  1:3139422/1‑150 (MQ=255)
                                                                       ggTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATc                                                                              <  2:77217/150‑1 (MQ=255)
                                                                         tCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCTCTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCgg                                                                            >  2:217189/1‑150 (MQ=255)
                                                                          cTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGt                                                                           <  2:2520145/150‑1 (MQ=255)
                                                                             gTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTAtc                                                                        <  1:1446334/150‑1 (MQ=255)
                                                                              tGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCTCTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTAtct                                                                       >  2:503795/1‑150 (MQ=255)
                                                                              tGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTAtct                                                                       >  2:1294181/1‑150 (MQ=255)
                                                                              tGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTAtct                                                                       <  1:1043972/150‑1 (MQ=255)
                                                                                aTTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCTCTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCa                                                                     >  2:4076919/1‑150 (MQ=255)
                                                                                    aaaaCCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTg                                                                 >  1:95464/1‑150 (MQ=255)
                                                                                    aaaaCCCGCTTGATGGTCTGCCCGCTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTg                                                                 <  1:3954285/150‑1 (MQ=255)
                                                                                     aaaCCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGt                                                                >  1:1169574/1‑150 (MQ=255)
                                                                                        cccGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCgg                                                             <  1:4537481/150‑1 (MQ=255)
                                                                                         ccGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGt                                                            >  2:1825366/1‑150 (MQ=255)
                                                                                           gCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGttt                                                          <  1:4144866/150‑1 (MQ=255)
                                                                                            cTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCTCTGGTAGGTCTGGCGGTGATCATCGGCTTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTg                                                         >  1:975092/1‑150 (MQ=255)
                                                                                            cTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTg                                                         <  1:1527165/150‑1 (MQ=255)
                                                                                             ttGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGa                                                        <  2:4292694/150‑1 (MQ=255)
                                                                                             ttGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGa                                                        <  2:2044570/150‑1 (MQ=255)
                                                                                               gATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCACTTGTCGGTTTGATc                                                      <  2:3521983/150‑1 (MQ=255)
                                                                                               gATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATc                                                      <  1:412230/150‑1 (MQ=255)
                                                                                               gATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATc                                                      <  2:2978179/150‑1 (MQ=255)
                                                                                                   gTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCg                                                  <  1:4289327/150‑1 (MQ=255)
                                                                                                   gTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCg                                                  >  1:589408/1‑150 (MQ=255)
                                                                                                   gTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCg                                                  <  1:1347671/150‑1 (MQ=255)
                                                                                                   gTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCg                                                  >  1:1000166/1‑150 (MQ=255)
                                                                                                     cTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGAt                                                <  1:3559118/150‑1 (MQ=255)
                                                                                                      tGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATc                                               >  1:4079313/1‑150 (MQ=255)
                                                                                                      tGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATc                                               <  1:2044274/150‑1 (MQ=255)
                                                                                                          cGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTaa                                           >  2:1429655/1‑150 (MQ=255)
                                                                                                          cGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTaa                                           >  2:3694190/1‑150 (MQ=255)
                                                                                                           gTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCTCTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAAc                                          >  2:4380652/1‑150 (MQ=255)
                                                                                                           gTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCTCTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAAc                                          >  2:4380470/1‑150 (MQ=255)
                                                                                                              gcgcTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTc                                       >  1:736804/1‑150 (MQ=255)
                                                                                                                  tGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCAtt                                   >  2:2675035/1‑150 (MQ=255)
                                                                                                                     gTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTtct                                <  2:473836/150‑1 (MQ=255)
                                                                                                                       gATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTcc                              >  1:2813377/1‑150 (MQ=255)
                                                                                                                       gATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTcc                              >  1:529140/1‑150 (MQ=255)
                                                                                                                       gATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCAATCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTcc                              >  1:3695325/1‑150 (MQ=255)
                                                                                                                         tttCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCgg                            >  1:3357865/1‑150 (MQ=255)
                                                                                                                               gACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGtttt                      >  2:4438120/1‑150 (MQ=255)
                                                                                                                               gACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGtttt                      >  2:1809331/1‑150 (MQ=255)
                                                                                                                               gACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGtttt                      >  1:4006571/1‑150 (MQ=255)
                                                                                                                                aCCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTc                     <  1:2271446/150‑1 (MQ=255)
                                                                                                                                 ccTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCg                    <  1:2481128/150‑1 (MQ=255)
                                                                                                                                     cccGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCAt                <  2:392816/150‑1 (MQ=255)
                                                                                                                                        ggTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGcc             >  2:1158748/1‑150 (MQ=255)
                                                                                                                                         gTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCt            >  2:795352/1‑150 (MQ=255)
                                                                                                                                          tGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCtt           >  1:3390836/1‑150 (MQ=255)
                                                                                                                                          tGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGATTTCGCCATGCCtt           >  2:4040794/1‑150 (MQ=255)
                                                                                                                                           gATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTc          <  1:3794291/150‑1 (MQ=255)
                                                                                                                                               aTGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTg      >  2:4504846/1‑150 (MQ=255)
                                                                                                                                               aTGTCACGGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTg      <  2:3336961/150‑1 (MQ=255)
                                                                                                                                                 gTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAg    <  1:197969/150‑1 (MQ=255)
                                                                                                                                                 gTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAg    <  2:3954392/150‑1 (MQ=255)
                                                                                                                                                 gTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAg    <  2:652435/150‑1 (MQ=255)
                                                                                                                                                 gTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAg    >  2:282598/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCG  >  CP009273/4268819‑4269115

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: