Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,340,386 G→T 10.5% intergenic (+97/+31) ubiG → / ← yfaL bifunctional 3‑demethylubiquinone‑8 3‑O‑methyltransferase and 2‑octaprenyl‑6‑hydroxyphenol methylase/putative autotransporter adhesin YfaL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,340,3860GT10.5% 222.0 / 10.9 72intergenic (+97/+31)ubiG/yfaLbifunctional 3‑demethylubiquinone‑8 3‑O‑methyltransferase and 2‑octaprenyl‑6‑hydroxyphenol methylase/putative autotransporter adhesin YfaL
Reads supporting (aligned to +/- strand):  ref base G (35/29);  new base T (2/6);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 1.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

ACTCGGCCCCGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTGCCACGCCACGCTGCCGCGCAGGGACACTCGCTGACTGATATTGCCCGTGACT  >  NC_000913/2340238‑2340531
                                                                                                                                                    |                                                                                                                                                 
aCTCGGCCCCGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTAcgccgc                                                                                                                                                  <  1:362508/150‑1 (MQ=255)
    ggCCCCGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                              <  1:936620/150‑1 (MQ=38)
    ggCCCCGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                              <  1:941310/150‑1 (MQ=38)
     gCCCCGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCg                                                                                                                                             >  2:125315/1‑150 (MQ=255)
      ccccGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                                                                                            <  1:1064032/150‑1 (MQ=38)
      ccccGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                                                                                            >  1:754488/1‑150 (MQ=38)
      ccccGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                                                                                                                                            <  1:406347/150‑1 (MQ=255)
       cccGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                                                                                           <  1:1157470/150‑1 (MQ=255)
       cccGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                                                                                           >  1:580776/1‑150 (MQ=255)
        ccGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                          <  1:1415064/150‑2 (MQ=38)
        ccGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                          <  2:754488/150‑2 (MQ=38)
                  tgtgAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGa                                                                                                                                >  2:534337/1‑150 (MQ=255)
                    tgAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCTTAGGCCGGATAAGGCGTTTACGCCTCATCCGGCAACCGTTCCGACt                                                                                                                              <  2:1153244/150‑1 (MQ=255)
                    tgAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACt                                                                                                                              <  1:1233644/150‑1 (MQ=255)
                      aaCTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAg                                                                                                                            <  1:1072062/150‑1 (MQ=255)
                         tataTGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACa                                                                                                                         >  2:481785/1‑150 (MQ=255)
                          atatGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTTGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAg                                                                                                                        >  2:413323/1‑150 (MQ=255)
                           tatGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGt                                                                                                                       >  1:386010/1‑150 (MQ=255)
                            atGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTg                                                                                                                      <  2:454327/150‑1 (MQ=255)
                             tgctgcACACGCAGAATAAGTGAGGTTCATGTTTGGACGCGCCAATGCCTGATGCGAAGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGa                                                                                                                     >  2:950917/1‑150 (MQ=255)
                              gctgtACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGAt                                                                                                                    >  2:317319/1‑150 (MQ=255)
                              gctgcACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGAt                                                                                                                    >  2:808352/1‑150 (MQ=255)
                              gctgcACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCATGCCGACTAGACAGTGAt                                                                                                                    >  2:363706/1‑150 (MQ=255)
                                  cacaCGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTAcc                                                                                                                <  2:100912/150‑1 (MQ=255)
                                  cacaCGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTAcc                                                                                                                <  1:629753/150‑1 (MQ=255)
                                      cGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCAttt                                                                                                            >  1:871534/1‑150 (MQ=255)
                                                 tGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCg                                                                                                 <  2:871534/150‑1 (MQ=255)
                                                    ggTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACa                                                                                              >  2:523327/1‑150 (MQ=255)
                                                    ggTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACa                                                                                              >  2:238391/1‑150 (MQ=255)
                                                       tGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACaaaa                                                                                           >  2:854721/1‑150 (MQ=255)
                                                       agATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACaaaa                                                                                           >  2:403882/2‑150 (MQ=255)
                                                              tGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGc                                                                                    >  1:1147390/1‑150 (MQ=255)
                                                                    cgcCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTg                                                                              <  1:1162315/150‑1 (MQ=255)
                                                                    cgcCAATGCCTGATGCGACGCTTGCATCGTCTTATCAGGCCTACAAATGCTCCCCGTAGTCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTg                                                                              <  2:1017955/150‑1 (MQ=255)
                                                                     gcCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTgg                                                                             <  2:729227/150‑1 (MQ=255)
                                                                      ccAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTggg                                                                            <  2:596645/150‑1 (MQ=255)
                                                                      ccAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTggg                                                                            <  1:1040845/150‑1 (MQ=255)
                                                                      ccAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTggg                                                                            <  2:127595/150‑1 (MQ=255)
                                                                       cAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGc                                                                           >  1:940643/1‑150 (MQ=255)
                                                                          tGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCaaa                                                                        <  1:950917/150‑1 (MQ=255)
                                                                            ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAat                                                                      >  1:297359/1‑150 (MQ=255)
                                                                                aTGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAatcatc                                                                  <  1:125315/150‑1 (MQ=255)
                                                                                  gCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGAc                                                                                                                          >  2:761253/1‑92 (MQ=255)
                                                                                  gCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGAc                                                                                                                          <  1:761253/92‑1 (MQ=255)
                                                                                    gACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTc                                                              >  1:1176126/1‑150 (MQ=255)
                                                                                     aCGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTcc                                                             >  2:347855/1‑150 (MQ=255)
                                                                                        cTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACAGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCtt                                                          >  2:1296500/1‑150 (MQ=255)
                                                                                           gCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGACAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAatcatc                                                                  >  2:365698/1‑139 (MQ=255)
                                                                                           gCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGACAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAatcatc                                                                  <  1:365698/139‑1 (MQ=255)
                                                                                           gCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCAGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGACATCGACAAAAACCCTTACGTCTTGGCAAAATCATCACTAACTTTCt                                                       >  2:875530/1‑150 (MQ=38)
                                                                                            ccGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTg                                                      >  2:696456/1‑150 (MQ=255)
                                                                                            ccGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTg                                                      <  2:386010/150‑1 (MQ=255)
                                                                                             cgcgTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTgc                                                     <  1:808352/150‑1 (MQ=255)
                                                                                              gcgTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTgcc                                                    >  2:289223/1‑150 (MQ=255)
                                                                                              gcgTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCATGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTgcc                                                    <  1:363706/150‑1 (MQ=255)
                                                                                                  cTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTgccacgc                                                >  2:707362/1‑150 (MQ=255)
                                                                                                   ttATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTATACAGTGATTAACATTTCACCGACATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTgccacgcc                                               >  2:591019/1‑150 (MQ=255)
                                                                                                   ttATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTgccacgcc                                               <  1:1300869/150‑1 (MQ=255)
                                                                                                    tATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTgccacgcca                                              >  1:1308305/1‑150 (MQ=255)
                                                                                                        aGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTGCCACGCCACGCt                                          <  1:317319/150‑1 (MQ=255)
                                                                                                            cTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTGCCACGCCACGCTGCcg                                      <  2:1147390/150‑1 (MQ=255)
                                                                                                               cAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTGCCACGCCACGCTGCCgcgc                                   >  2:153298/1‑150 (MQ=255)
                                                                                                                   tGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTGCCACGCCACGCTGCCGCGCAggg                               <  2:1308305/150‑1 (MQ=255)
                                                                                                                        cccGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTGCCACGCCACGCTGCCGCGCAGGGACACt                          >  1:609589/1‑150 (MQ=255)
                                                                                                                          cGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTGCCACGCCACGCTGCCGCGCAGGGACACTCg                        >  2:523551/1‑150 (MQ=255)
                                                                                                                           ggttgCCGGATGCGGCGTAAACGCCTTATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTGCCACGCCACGCTGCCGCGCAGGGACACTCGc                       >  1:387525/5‑150 (MQ=255)
                                                                                                                               gCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTGCCACGCCACGCTGCCGCGCAGGGACACTCGctgac                   >  2:307440/1‑150 (MQ=255)
                                                                                                                               gCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTGCCACGCCACGCTGCCGCGCAGGGACACTCGctgac                   <  1:363708/150‑1 (MQ=255)
                                                                                                                                 cGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTGCCACGCCACGCTGCCGCGCAGGGACACTCGctgactg                 <  1:684691/150‑1 (MQ=255)
                                                                                                                                 cGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTGCCACGCCACGCTGCCGCGCAGGGACACTCGctgactg                 <  1:403882/150‑1 (MQ=255)
                                                                                                                                    aTAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTGCCACGCCACGCTGCCGCGCAGGGACACTCGCTGACTGata              >  2:652479/1‑150 (MQ=255)
                                                                                                                                            tttACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTGCCACGCCACGCTGCCGCGCAGGGACACTCGCTGACTGATATTGCCCGt      <  1:981398/150‑1 (MQ=255)
                                                                                                                                              tACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTGCCACGCCACGCTGCCGCGCAGGGACACTCGCTGACTGATATTGCCCGTGa    >  2:459393/1‑150 (MQ=255)
                                                                                                                                                cgccgcATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTGCCACGCCACGCTGCCGCGCAGGGACACTCGCTGACTGATATTGCCCGTGACt  >  2:79518/1‑150 (MQ=255)
                                                                                                                                                cgccgcATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTGCCACGCCACGCTGCCGCGCAGGGACACTCGCTGACTGATATTGCCCGTGACt  <  2:838523/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
ACTCGGCCCCGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTGCCACGCCACGCTGCCGCGCAGGGACACTCGCTGACTGATATTGCCCGTGACT  >  NC_000913/2340238‑2340531

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: