Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,054,222 T→A 8.3% V139D (GTT→GAT)  yihL → putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,054,2220TA 8.3% 1043.9 / 67.1 313V139D (GTT→GAT) yihLputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (130/157);  new base A (11/15);  total (141/172)
Fisher's exact test for biased strand distribution p-value = 8.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

Only 100 of 318 total aligned reads displayed.

AATTGAGCAAGGGCGAGAGCCTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAACCCGTCAGAATGCCGG  >  CP009273/4054073‑4054371
                                                                                                                                                      |                                                                                                                                                     
aaTTGAGCAAGGGCGAGAGCCTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGt                                                                                                                                                       <  1:3880088/150‑1 (MQ=255)
 aTTGAGCAAGGGCGAGAGCCTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGtt                                                                                                                                                      <  2:4727380/150‑1 (MQ=255)
    gAGCAAGGGCGAGAGCCTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTc                                                                                                                                                   <  2:526463/150‑1 (MQ=255)
    gAGCAAGGGCGAGAGCCTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTc                                                                                                                                                   <  2:592595/150‑1 (MQ=255)
    gAGCAAGGGCGAGAGCCTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGATTTc                                                                                                                                                   <  1:4793539/150‑1 (MQ=255)
          gggCGAGAGCCTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCAc                                                                                                                                             <  1:4623073/150‑1 (MQ=255)
           ggCGAGAGGCTTCCTGGGGGTTTACCGAGTAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCAct                                                                                                                                            <  2:2702698/150‑2 (MQ=255)
             cGAGAGCCTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGaa                                                                                                                                          >  1:105594/1‑150 (MQ=255)
              gagaGCCTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGaaa                                                                                                                                         >  2:2432476/1‑150 (MQ=255)
              gagaGCCTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGaaa                                                                                                                                         >  1:2637335/1‑150 (MQ=255)
               agagCCTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAAc                                                                                                                                        <  2:4845606/150‑1 (MQ=255)
               agagCCTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAAc                                                                                                                                        >  2:3251690/1‑150 (MQ=255)
                gagCCTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACa                                                                                                                                       <  2:105576/150‑1 (MQ=255)
                   ccTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACatat                                                                                                                                    <  1:3508129/150‑1 (MQ=255)
                   ccTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGATTTCTATCACGAAACatat                                                                                                                                    <  2:4402558/150‑1 (MQ=255)
                      tCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATAtt                                                                                                                                 >  2:411185/1‑150 (MQ=255)
                       ccTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTa                                                                                                                                <  1:1123305/150‑1 (MQ=255)
                         tGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAAt                                                                                                                              <  2:5031876/150‑1 (MQ=255)
                           ggggTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATcc                                                                                                                            >  1:3398315/1‑150 (MQ=255)
                               tttACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGaa                                                                                                                        <  2:3749397/150‑1 (MQ=255)
                                ttACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAg                                                                                                                       <  2:1337816/150‑1 (MQ=255)
                                ttACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAg                                                                                                                       >  1:1718465/1‑150 (MQ=255)
                                 tACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGt                                                                                                                      >  2:2897080/1‑150 (MQ=255)
                                  aCCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGtt                                                                                                                     <  1:1343373/150‑1 (MQ=255)
                                   ccGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTg                                                                                                                    >  1:2120332/1‑150 (MQ=255)
                                    cGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGATTTCTATCACGAAACATATATTAATCCTGAAGTTGc                                                                                                                   >  1:736908/1‑150 (MQ=255)
                                     gagaAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCt                                                                                                                  >  2:4980231/1‑150 (MQ=255)
                                     gagaAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGATTTCTATCACGAAACATATATTAATCCTGAAGTTGCt                                                                                                                  >  2:2880580/1‑150 (MQ=255)
                                      agaAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTc                                                                                                                 <  2:3873742/150‑1 (MQ=255)
                                       gaAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTcc                                                                                                                <  1:4836595/150‑1 (MQ=255)
                                         aaaaCCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCgg                                                                                                              <  1:2313538/150‑1 (MQ=255)
                                         aaaaCCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCgg                                                                                                              <  2:1166663/150‑1 (MQ=255)
                                             ccGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAAATGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCtt                                                                                                          <  1:3208334/150‑1 (MQ=255)
                                              cGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCttt                                                                                                         <  2:1571665/150‑1 (MQ=255)
                                               gTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTa                                                                                                        >  2:1632129/1‑150 (MQ=255)
                                                tACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTAt                                                                                                       <  2:1016208/150‑1 (MQ=255)
                                                tACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTAt                                                                                                       <  1:4714958/150‑1 (MQ=255)
                                                 aCCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTAtt                                                                                                      >  1:4804980/1‑150 (MQ=255)
                                                   cAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGa                                                                                                    <  1:2713038/150‑1 (MQ=255)
                                                   cAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGa                                                                                                    <  1:1214734/150‑1 (MQ=255)
                                                      cGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGATTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGaaca                                                                                                 >  1:1078448/1‑150 (MQ=255)
                                                       gATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGaacaa                                                                                                >  2:640353/1‑150 (MQ=255)
                                                       gATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGaacaa                                                                                                <  1:3542794/150‑1 (MQ=255)
                                                       gATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGaacaa                                                                                                >  2:4262514/1‑150 (MQ=255)
                                                        atatTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGaacaac                                                                                               <  2:1478165/150‑1 (MQ=255)
                                                         tatTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACt                                                                                              >  1:3949563/1‑150 (MQ=255)
                                                          atTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACtt                                                                                             <  1:164946/150‑1 (MQ=255)
                                                             cccGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGaa                                                                                          >  1:654980/1‑150 (MQ=255)
                                                               cGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGaaaa                                                                                        >  2:4362383/1‑150 (MQ=255)
                                                                gagtCGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCTTCACGGGCTGGGGGGCGCTGGAAGGACGTAAAGTTTTCTATCACGATACATTTATTATTCCGGAAGTTGCTCCGGGCTTTTTTGAACAACTTGAAAAc                                                                                       <  2:1709819/150‑1 (MQ=255)
                                                                gagaCGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAAc                                                                                       <  2:3970213/150‑1 (MQ=255)
                                                                gagaCGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAAc                                                                                       >  2:2776236/1‑150 (MQ=255)
                                                                 agaCGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAAcc                                                                                      <  1:1891679/150‑1 (MQ=255)
                                                                  gaCGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCa                                                                                     >  1:992755/1‑150 (MQ=255)
                                                                   aCGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCAc                                                                                    <  1:5071984/150‑1 (MQ=255)
                                                                    cGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACt                                                                                   <  1:4432016/150‑1 (MQ=255)
                                                                    cGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACt                                                                                   >  1:2918888/1‑150 (MQ=255)
                                                                      cTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCa                                                                                 >  1:1447840/1‑150 (MQ=255)
                                                                        cgcgCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCAtt                                                                               <  1:816235/150‑1 (MQ=255)
                                                                        cgcgCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCAtt                                                                               <  2:2374704/150‑1 (MQ=255)
                                                                        cgcgCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCAtt                                                                               >  2:818603/1‑150 (MQ=255)
                                                                        cgcgCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCAtt                                                                               <  2:835722/150‑1 (MQ=255)
                                                                         gcgcCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCAttt                                                                              <  1:2292015/150‑1 (MQ=255)
                                                                           gcCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTc                                                                            <  1:2783173/150‑1 (MQ=255)
                                                                                  aTTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTc                                                                     <  2:3645188/150‑1 (MQ=255)
                                                                                   ttGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCt                                                                    <  2:3397586/150‑1 (MQ=255)
                                                                                    tGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTg                                                                   <  2:2456210/150‑1 (MQ=255)
                                                                                       aGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGa                                                                >  1:2038083/1‑150 (MQ=255)
                                                                                         tGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGaaa                                                              <  1:1384652/150‑1 (MQ=255)
                                                                                          gACGCCATCAACTCAACTCATCGGAATCACGGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGaaaa                                                             <  1:3607753/150‑1 (MQ=255)
                                                                                           aCGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAg                                                            >  2:683066/1‑150 (MQ=255)
                                                                                               cATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCt                                                        >  1:1035512/1‑150 (MQ=255)
                                                                                               cATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCt                                                        >  1:1972552/1‑150 (MQ=255)
                                                                                                aTCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTa                                                       >  2:3312203/1‑150 (MQ=255)
                                                                                                aTCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTa                                                       >  2:3312247/1‑150 (MQ=255)
                                                                                                   aaCTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCa                                                    <  1:907811/150‑1 (MQ=255)
                                                                                                         aCTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAaga                                              <  2:208892/150‑1 (MQ=255)
                                                                                                           tctATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAgaga                                            >  1:1096373/1‑150 (MQ=255)
                                                                                                            ctATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGAc                                           <  2:1961339/150‑1 (MQ=255)
                                                                                                             tATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACg                                          <  1:4251724/150‑1 (MQ=255)
                                                                                                              aTCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACgg                                         >  1:308944/1‑150 (MQ=255)
                                                                                                                cGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAAATTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGgta                                       <  2:3009506/150‑1 (MQ=255)
                                                                                                                  gCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGgtag                                      >  1:878503/1‑150 (MQ=255)
                                                                                                                    aTCACCGGCTGGGGGGCGCTGGAAGGACATAAAGATTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGgtagta                                    <  2:2306653/150‑1 (MQ=255)
                                                                                                                     tCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTaa                                   >  1:679747/1‑150 (MQ=255)
                                                                                                                      cACCGGCTGGGGGGCGATGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTaaa                                  >  1:2356286/1‑150 (MQ=255)
                                                                                                                             tGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTg                           <  1:4150695/150‑1 (MQ=255)
                                                                                                                              ggggggCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGa                          >  1:1008515/1‑150 (MQ=255)
                                                                                                                               gggggCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGAt                         >  2:4585246/1‑150 (MQ=255)
                                                                                                                                  ggCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGAtttt                      <  1:4915041/150‑1 (MQ=255)
                                                                                                                                    cgcTGGTAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCa                    <  1:4202473/150‑1 (MQ=255)
                                                                                                                                      cTGGAAGGACATAAAGATTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCaaa                  <  1:3214856/150‑1 (MQ=255)
                                                                                                                                       tGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAAc                 <  1:188387/150‑1 (MQ=255)
                                                                                                                                         gAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAAccc               >  2:3072968/1‑150 (MQ=255)
                                                                                                                                          aaGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAACCCg              <  1:4127192/150‑1 (MQ=255)
                                                                                                                                          aaGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAACCCg              >  1:2867470/1‑150 (MQ=255)
                                                                                                                                              aCATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAACCCGTCAg          >  2:2632722/1‑150 (MQ=255)
                                                                                                                                              aCATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACAACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTCGTAAAAAAATTGATTTTCAAACCCGTCAg          >  2:3532246/1‑150 (MQ=255)
                                                                                                                                                  aaaGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAACCCGTCAGAATg      >  2:4929312/1‑150 (MQ=255)
                                                                                                                                                     gTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAACCCGTCAGAATGCCg   >  2:4236826/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
AATTGAGCAAGGGCGAGAGCCTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAACCCGTCAGAATGCCGG  >  CP009273/4054073‑4054371

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: