Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,340,373 A→C 8.9% intergenic (+84/+44) 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,3730AC 8.9% 230.0 / 14.5 79intergenic (+84/+44)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 A (37/35);  new base C (2/5);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.33e-01

AATACTTTTAAACTCGGCCCCGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTGCCACGCCACGCTGCCGCGCAGGGACACTCGCTGACTGATA  >  NC_000913/2340227‑2340519
                                                                                                                                                  |                                                                                                                                                  
aaTACTTTTAAACTCGGCCCCGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGc                                                                                                                                                 >  1:838523/1‑150 (MQ=255)
 aTACTTTTAAACTCGGCCCCGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCg                                                                                                                                                <  1:1126418/150‑1 (MQ=255)
   aCTTTTAAACTCGGCCCCGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGtt                                                                                                                                              <  2:543915/150‑1 (MQ=255)
   aCTTTTAAACTCGGCCCCGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGtt                                                                                                                                              <  1:1337162/150‑1 (MQ=255)
     ttttAAACTCGGCCCCGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTa                                                                                                                                            >  1:100912/1‑150 (MQ=255)
     ttttAAACTCGGCCCCGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTa                                                                                                                                            >  1:596645/1‑150 (MQ=255)
        tAAACTCGGCCCCGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTAcgc                                                                                                                                         <  2:757272/150‑1 (MQ=255)
         aaaCTCGGCCCCGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTAcgcc                                                                                                                                        <  1:713208/150‑1 (MQ=255)
           aCTCGGCCCCGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTAcgccgc                                                                                                                                      <  1:362508/150‑1 (MQ=255)
               ggCCCCGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                  <  1:941310/150‑1 (MQ=38)
               ggCCCCGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                  <  1:936620/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                                                                                                    <  1:629753/150‑1 (MQ=255)
                                             cacaCGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTAcc                                                                                                    <  2:100912/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:386010/150‑1 (MQ=255)
                                                                                                       ccGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTg                                          >  2:696456/1‑150 (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       <  1:363708/150‑1 (MQ=255)
                                                                                                                                          gCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTGCCACGCCACGCTGCCGCGCAGGGACACTCGctgac       >  2:307440/1‑150 (MQ=255)
                                                                                                                                            cGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTGCCACGCCACGCTGCCGCGCAGGGACACTCGctgactg     <  1:403882/150‑1 (MQ=255)
                                                                                                                                            cGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTGCCACGCCACGCTGCCGCGCAGGGACACTCGctgactg     <  1:684691/150‑1 (MQ=255)
                                                                                                                                               aTAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTGCCACGCCACGCTGCCGCGCAGGGACACTCGCTGACTGata  >  2:652479/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AATACTTTTAAACTCGGCCCCGGCGTGGATGTGAACTATATGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAATGCTCCCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCACCGTCATCGACAAAAACCCTGCCGTCTGGGCAAAATCATCACTCCCTTTCTGCCACGCCACGCTGCCGCGCAGGGACACTCGCTGACTGATA  >  NC_000913/2340227‑2340519

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: