Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,296,060 C→T 40.3% intergenic (+266/+376) gltP → / ← yjcO glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,296,0600CT40.3% 54.2 / 97.0 91intergenic (+266/+376)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base C (27/26);  new base T (26/10);  total (55/36)
Fisher's exact test for biased strand distribution p-value = 5.10e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.60e-01

GTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCC  >  NC_000913/4295912‑4296207
                                                                                                                                                    |                                                                                                                                                   
gTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACg                                                                                                                                                    >  2:1519108/1‑150 (MQ=32)
  gggATAAGGCGTTTTTGCCGGATCAGAAATCAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGcc                                                                                                                                                  <  2:606228/149‑1 (MQ=11)
   ggATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCt                                                                                                                                                 <  2:1754014/150‑1 (MQ=32)
    gATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCAGACATCAATGCCTg                                                                                                                                                >  2:310485/1‑150 (MQ=17)
       aacgagTTTATGCCGAAAACGACATCAACGCCTGATCCGACGCTTCACGCGTCTAATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGGCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATg                                                                                                                                             <  1:1792551/145‑1 (MQ=11)
        aGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGc                                                                                                                                            <  2:2004317/150‑1 (MQ=17)
        aGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGc                                                                                                                                            >  1:178065/1‑150 (MQ=17)
        aGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGc                                                                                                                                            <  1:859825/150‑1 (MQ=32)
           cGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGAc                                                                                                                                         >  1:793382/1‑150 (MQ=17)
           cGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGACAGACAGCGCAATAGCCAGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCAGCATCAGACATCAATGCCAGAtaagac                                                                                                                                         >  1:114749/1‑145 (MQ=11)
              ttATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTATGTAGGTCGGATAAGGAGTTTACGCCCCATCCGAAATCAACGCCAGAGGCCAAGCt                                                                                                                                      >  1:1112783/1‑150 (MQ=11)
              ttATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTCAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGTTGCGACGCt                                                                                                                                      >  2:703381/1‑150 (MQ=12)
              ttATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAAAGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAGTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCTActct                                                                                                                                      >  1:1273838/1‑150 (MQ=11)
                aTGCCGCATACGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTg                                                                                                                                    <  1:1718296/150‑1 (MQ=17)
                 tGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGAAAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGc                                                                                                                                   >  1:1525148/1‑150 (MQ=18)
                  gCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATATCCTGATTTAGCGTGATTTTGTAGGTCGCATAAGGCGTTTACGACGCATCCGACATCAATGCCTGATGCGACGCTTgac                                                                                                                                  >  2:635000/1‑148 (MQ=255)
                  gCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGcc                                                                                                                                  >  1:1946226/1‑150 (MQ=21)
                  gCCGCATCCGACATAACAGCCTGATGGGAGTCTGTAGGCGTCTTTTGTGGGCTACGCCTGTCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGct                                                                                                                                  <  1:116768/150‑2 (MQ=11)
                   ccGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCcg                                                                                                                                 >  1:1723667/1‑150 (MQ=25)
                    cGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCcgc                                                                                                                                >  1:402573/1‑150 (MQ=35)
                        tCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCt                                                                                                                            >  1:451283/1‑150 (MQ=25)
                        tCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCGGATATAGGGTGATTTTGTAGGTCGGATAAGGCGTTTAAGCCGCAAACGACATGAACGCCTGATGCGACGCGTGCCGCGTCt                                                                                                                            >  1:1290343/1‑150 (MQ=11)
                         ccGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCCGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAAGCCTGATGCGACGCTTGCCGCGTCtt                                                                                                                           >  1:1928440/1‑150 (MQ=255)
                         ccGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGAGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGACGCGTCtt                                                                                                                           >  2:1598487/1‑150 (MQ=12)
                          cGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTa                                                                                                                          >  2:1660302/1‑150 (MQ=25)
                          cGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAAGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTa                                                                                                                          >  2:1742412/1‑150 (MQ=17)
                           gACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                                                         >  2:797541/1‑150 (MQ=25)
                            tcATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACCGCTCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATc                                                                                                                        <  1:1220508/149‑1 (MQ=255)
                              aTCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAg                                                                                                                      >  1:485261/1‑150 (MQ=25)
                               tCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGACAGACAGCGCAATAGCCTGATTTAGCGTGAGTTTGTAGGTCGGATAAGGCGTTTACGCAGAATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                                                                                                     >  1:1032708/1‑150 (MQ=255)
                                cAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                    >  1:1600354/1‑150 (MQ=25)
                                cAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                    >  2:823829/1‑150 (MQ=25)
                                cAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGGCGGATCAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                    <  1:702142/150‑1 (MQ=255)
                                 aaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCcgc                                                                                                                                <  1:886401/137‑1 (MQ=25)
                                 aaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCcgc                                                                                                                                >  2:886401/1‑137 (MQ=25)
                                 aaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGACTACGCCAGCCAGCGCAATAGCCTGATTTAGAGTGATTTTGTAGGTCGGATAAGTCGTTTAAGCAGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCGGGcg                                                                                                                   >  2:1294603/1‑149 (MQ=11)
                                      cTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGAATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGACTatct                                                                                                              >  1:1882640/1‑147 (MQ=255)
                                       tGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatctt                                                                                                             >  1:44031/1‑146 (MQ=17)
                                       tGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatctt                                                                                                             >  1:1234120/1‑146 (MQ=17)
                                       tGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGcc                                                                                                             <  2:534973/150‑1 (MQ=37)
                                         aTGCGACGCTTGACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAg                                                                                                           <  2:710857/150‑1 (MQ=37)
                                         aTGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCTCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAg                                                                                                           <  1:1434950/150‑1 (MQ=255)
                                            cGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                   >  1:512510/1‑139 (MQ=25)
                                            cGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                   <  2:512510/139‑1 (MQ=25)
                                            cGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAgccc                                                                                                        >  1:1145467/1‑147 (MQ=37)
                                                    aaCGCGTCTTTTCCGGCCTACTCTAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATa                                                                                                <  1:1286485/150‑1 (MQ=255)
                                                       gcgTCTTATCAGGCCTACGCCAGCCTGATCAATAGCCTGATTTAGCGTGATTTTGTTGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGcc                                                                                             <  1:1925013/150‑1 (MQ=11)
                                                            ttATCAGGCCTACGCCAGACAGAGCTATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCTTTTACTCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGATAGCGCAATAGCCTGAtt                                                                                        <  2:1997147/150‑1 (MQ=11)
                                                                cAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                   <  2:691589/119‑1 (MQ=37)
                                                                cAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGc                                                                                    <  1:675794/150‑1 (MQ=37)
                                                                cAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAAGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAACAAAAACAGCGCAATATCCTGaaaaagc                                                                                    >  2:1571838/1‑144 (MQ=11)
                                                                cAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAAGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                   >  1:691589/1‑119 (MQ=16)
                                                                  ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATc                                                                                                                        <  2:492143/112‑1 (MQ=37)
                                                                  ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATc                                                                                                                        >  1:492143/1‑112 (MQ=39)
                                                                     cTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGcc                                                                                                                                                  <  2:2008566/83‑1 (MQ=18)
                                                                     cTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGcc                                                                                                                                                  >  1:2008566/1‑83 (MQ=18)
                                                                         gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCATGATGCGACGCTTGCCGCGTCTTATCAGGCCTAAGAAAGACCGCGGATTAGAGTGAATTCTCGTGagaatg                                                                           >  1:125033/1‑145 (MQ=11)
                                                                             gACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAgg                                                                       <  2:788762/150‑1 (MQ=37)
                                                                                   gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGc                                                                                                              >  2:1407311/1‑105 (MQ=40)
                                                                                   gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGc                                                                                                              <  1:1407311/105‑1 (MQ=40)
                                                                                       tatcttGATTTTGCTTGATTATGTAGGGCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAgg                                                             <  2:1112783/145‑1 (MQ=11)
                                                                                         gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTc                                                                                                                             <  2:1555292/84‑1 (MQ=255)
                                                                                         gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTc                                                                                                                             >  1:1555292/1‑84 (MQ=255)
                                                                                                 tAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                   <  1:678812/86‑1 (MQ=255)
                                                                                                 tAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                   >  2:678812/1‑86 (MQ=255)
                                                                                                   gCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTTCAGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCAGc                                                 >  2:1817576/1‑150 (MQ=255)
                                                                                                    ggTGATTTTGTTTTTCTGTTACGGCGTTTACTCCGGTTACGTCTTCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCa                                                <  2:1271004/149‑1 (MQ=11)
                                                                                                      tGATTTTGTTGGGCGGATATGGGGTTTACTCCGCTTCAGTCATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATc                                              <  2:1519388/150‑1 (MQ=11)
                                                                                                       gATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTg                                                                           <  1:1163911/120‑1 (MQ=17)
                                                                                                             gTAGGGCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCTTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATc                                       <  1:1519108/150‑1 (MQ=255)
                                                                                                              tAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAAGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAAGCCGCATCCGACATCa                                      >  2:518256/1‑150 (MQ=255)
                                                                                                                 gTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACg                                   >  1:1493842/1‑150 (MQ=37)
                                                                                                                   cGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAggt                                                                                                                    <  1:133486/67‑2 (MQ=255)
                                                                                                                   cGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                    >  2:133486/1‑67 (MQ=255)
                                                                                                                   cGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTAGGATAAGGCGTTTATGCCGCATCCGACATCAAAGcc                                 >  1:2004317/1‑150 (MQ=255)
                                                                                                                   cGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCGTGCCGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGGCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGcc                                 >  1:891106/1‑150 (MQ=255)
                                                                                                                     gATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGCTTATGCCGCATCCGACATCAACGCCTg                               >  1:1927510/1‑150 (MQ=37)
                                                                                                                       tAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAt                             >  1:848439/1‑150 (MQ=37)
                                                                                                                            cGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                   <  2:277280/59‑1 (MQ=255)
                                                                                                                            cGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                   >  1:277280/1‑59 (MQ=255)
                                                                                                                               ttACGCCGCATCCTACTTCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCt                     <  2:1127234/150‑1 (MQ=255)
                                                                                                                                     cgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACCGCGCAATAGCCTGATTTAGCGTGAGTTTCTAGGTCGGATAAGGCGTTTAAGCCGCATCCCACATCAACGCCTGATGCGACGCTTACcgc               >  1:116777/1‑150 (MQ=11)
                                                                                                                                      gcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAAcgcg              >  1:1757048/1‑150 (MQ=34)
                                                                                                                                        ctCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTc            <  1:1294603/149‑1 (MQ=21)
                                                                                                                                          ccGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCTTGAATATGTAGGTCGGATAAGGCGTTTAAGCCGCATCCGACATCAAAGCCTGATGCTACGCTTAACGCGTCtt          >  2:691634/1‑150 (MQ=11)
                                                                                                                                            gACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTaa        >  2:715584/1‑149 (MQ=21)
                                                                                                                                              cATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCa      >  1:630122/1‑150 (MQ=34)
                                                                                                                                               aTCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGcc                                                                                             <  1:663740/62‑1 (MQ=40)
                                                                                                                                                tCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGttt                                                        >  2:1301627/1‑98 (MQ=35)
                                                                                                                                                 cAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGc   >  2:808521/1‑150 (MQ=32)
                                                                                                                                                  aaCGCCTGATGCGTCGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGcc  <  1:1420828/150‑1 (MQ=17)
                                                                                                                                                    |                                                                                                                                                   
GTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCC  >  NC_000913/4295912‑4296207

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