Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,296,060 C→T 39.6% 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,0600CT39.6% 36.3 / 76.7 76intergenic (+266/+376)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base C (25/19);  new base T (18/11);  total (46/30)
Fisher's exact test for biased strand distribution p-value = 8.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.06e-02

GGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCT  >  NC_000913/4295915‑4296208
                                                                                                                                                 |                                                                                                                                                    
ggATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAAGCCt                                                                                                                                                  >  1:14352/1‑150 (MQ=17)
   aatggCGTTTATGCCGCTTCCGTCTTCAAGGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGTATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGAt                                                                                                                                               <  2:406646/147‑1 (MQ=11)
     tggCTTTTATTCCGCTTCAGTCTTCAACTCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACGCCAGACAGAGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGc                                                                                                                                             <  1:832458/149‑1 (MQ=11)
     tggCGTTTATGCCGCATCAGTTTTCAATGACTGATGCGACGCTTTACGCGTATTATCAGGCCTACGCCAGACAGCGCATTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGc                                                                                                                                             <  2:390067/149‑1 (MQ=11)
        ggTTTATTCCGCATACGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGAc                                                                                                                                          <  2:602138/149‑1 (MQ=17)
        cGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGAc                                                                                                                                          >  1:314371/1‑150 (MQ=17)
           ttATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCt                                                                                                                                       <  2:729624/150‑1 (MQ=32)
            tATGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGAGTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCTTCCGACATCACTGCCTGATGCGACGCtt                                                                                                                                      >  2:679655/1‑150 (MQ=11)
            tATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTAAGAAAGACAGCGCAATAGCATGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAAGCCTGATGCGACGCtt                                                                                                                                      >  1:516243/1‑150 (MQ=11)
             ttGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTg                                                                                                                                     <  2:385622/149‑1 (MQ=17)
             accccGCACCCGACCTCAACGCCTGAAGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACATCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTg                                                                                                                                     <  1:764352/147‑1 (MQ=11)
             aTGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTg                                                                                                                                     >  2:899279/1‑150 (MQ=17)
               gCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGcc                                                                                                                                   >  2:486059/1‑150 (MQ=21)
               gCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTgac                                                                                                                                   >  1:1038324/1‑148 (MQ=17)
                ccGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCcg                                                                                                                                  >  1:386518/1‑150 (MQ=25)
                   cATCCGCCATCAACGCCTGTTGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGCCCGCGCAATTGCCTGATTTTGCGTGATTTTGTAGGTAGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGCCGCTTGCCGCGt                                                                                                                               >  2:926925/1‑150 (MQ=11)
                   cATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGt                                                                                                                               >  1:19137/1‑150 (MQ=25)
                   cATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCATGCcgcgg                                                                                                                               >  1:226806/1‑149 (MQ=25)
                   cATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGt                                                                                                                               >  1:296810/1‑150 (MQ=37)
                        gtcATAAAAGCCTGATGCGACGCTTAACGCGTCTTATCAGGCATAAGCCAGACAGCGCAATTGCATGATTTAGCGTGATTTTGTAGGTCGGATAAGGCTTTTACGCCGCATAAGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                                                          <  2:823236/148‑1 (MQ=11)
                        gACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGAGACGCTTGCCGCGTCTTAt                                                                                                                          >  2:982022/1‑150 (MQ=25)
                          cATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCa                                                                                                                        >  2:43255/1‑150 (MQ=25)
                            tCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                                                                                                      <  1:1067498/150‑1 (MQ=37)
                             caagggatGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCATTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCTACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                     <  1:905884/143‑1 (MQ=11)
                              aaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAgggc                                                                                                                    >  1:21539/1‑148 (MQ=37)
                               aCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCt                                                                                                                   >  1:1255542/1‑150 (MQ=25)
                                cGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTa                                                                                                                  >  2:1198937/1‑150 (MQ=25)
                                  ccTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGCCGCTTGCCGCGTCTTATCAGGCCTatc                                                                                                                >  2:676690/1‑148 (MQ=16)
                                   cTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACga                                                                                                               >  1:305209/1‑149 (MQ=37)
                                    tGATGCGACGCTTAACGCGTCTTATCAGGGCTACGCCAGACAGAGCTATTGTCTGATTTAGCTTGATTTTGTAGGGCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGcc                                                                                                              <  2:549546/150‑1 (MQ=11)
                                     gATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGc                                                                                                                                                    <  1:1203448/111‑1 (MQ=17)
                                     gATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGc                                                                                                                                                    <  1:1198408/111‑1 (MQ=17)
                                     gATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGc                                                                                                                                                    >  2:1203448/1‑111 (MQ=32)
                                     gATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGc                                                                                                                                                    >  2:1198408/1‑111 (MQ=32)
                                      ccgCGACGCGTTACGCGCCTTATCAGGCCTACGCCACACAGCGCAATCGCCTGATTTAGCGTGATTTTGTATTGCGGTTATGGCTTTTACGCCGCTTAAGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAg                                                                                                            <  2:176851/148‑1 (MQ=11)
                                         cGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGAc                                                                                                                                          >  2:928670/1‑117 (MQ=17)
                                         cGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGAc                                                                                                                                          <  1:928670/117‑1 (MQ=17)
                                          gACGCTTTCCGGGTCTTATGAGGCCTACGCCAGACAGCGCAATATCCTGATTTAGCGTGATTTTGTAGGTCGTATAAGGCGTTTACGCCGGATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAg                                                                                                        <  2:785131/150‑1 (MQ=11)
                                          gACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAg                                                                                                        <  1:274794/150‑1 (MQ=37)
                                                tAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTa                                                                                                                           <  1:196115/125‑1 (MQ=25)
                                                tAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTa                                                                                                                           <  1:196116/125‑1 (MQ=25)
                                                tAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTa                                                                                                                           >  2:196116/1‑125 (MQ=25)
                                                tAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAAGCCTGATGCGACGCTTGCCGCGTCTTa                                                                                                                           >  2:196115/1‑125 (MQ=38)
                                                  aCGCGTCTTATCAGGCCTACGCCAGACAGCGCAAAAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGATTAAGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCTCGTAATATAAGGCCTAAACAACACTGCGCAATAg                                                                                                >  1:958311/1‑150 (MQ=11)
                                                   cgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTc                                                                                                                              <  2:217408/119‑1 (MQ=37)
                                                   cgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTc                                                                                                                              >  1:217408/1‑119 (MQ=37)
                                                     cgTCTTATGAGGCATACGCGAGACAGCGAAATTGCGTTATTTAGGGTGATTTTGTATGTAGTTTATGGCTTTTACCCCGCCTCAGTTATTATCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCt                                                                                             <  2:231627/150‑1 (MQ=2)
                                                      gTCTTATAAGGCCTACGCCAGACATCGCAATAGACTGATTTAGCGTGATTTTGTTTGTCGTATATGGCGTTTACTCCGCATAATTGTTCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTg                                                                                            <  2:773782/150‑1 (MQ=11)
                                                          tATCAGGCCTACGCCATACAGCGAAATTGCCTTATTTAGCGTGATTTTGTAGGTAGGATAAGGCTTTTACTCCGCATCCGTCATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGAttt                                                                                        <  2:1109371/150‑1 (MQ=11)
                                                            tCAGGCCTACTCCATACAGTGAAATATCCTTATTTATCGTGATTTTGTATGTGTGATATGGGTTTTACTCCGCTTAAGTCATGAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAg                                                                                      <  2:299403/150‑1 (MQ=11)
                                                                       ccAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCATACCGCTCAATAGCATTATATATCGTGAATTTGt                                                                           >  1:124255/1‑150 (MQ=11)
                                                                          gACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCAAGACAGCGCAATAGCATGATTTAGCGTGATTTTGTAgg                                                                        >  1:375041/1‑150 (MQ=255)
                                                                            cAGCGCAATAGCCTGATTTAGCGTGATTTTGTCGGTCGGATAAGGCGTTTACGCCGCATCCGACATCACCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCCGACAGCGCAAAAGACTGATTTAGCTTGAGTTTGTATGTc                                                                      >  1:328843/1‑150 (MQ=11)
                                                                               cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAg                                                                                                                       >  1:517164/1‑98 (MQ=16)
                                                                               cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAg                                                                                                                       <  2:517164/98‑1 (MQ=255)
                                                                                gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACCGAAGACAGAGCAATAGACTTAATTAGCGTCATTATGTAGATCGGaa                                                                  >  2:1058982/1‑149 (MQ=11)
                                                                                        cttaTTTAGCGTGATTTTGTAGGTCGGATACGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGtt                                                          <  1:405043/147‑1 (MQ=255)
                                                                                               aGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCg                                                   >  2:962565/1‑150 (MQ=37)
                                                                                                 cGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCATTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCa                                                 >  1:906383/1‑150 (MQ=37)
                                                                                                     aTTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCg                                             >  2:948712/1‑150 (MQ=21)
                                                                                                         tGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCt                                                                                                                             <  1:160605/66‑1 (MQ=255)
                                                                                                         tGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCt                                                                                                                             >  2:160605/1‑66 (MQ=255)
                                                                                                         tGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACACCGCATTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATACGGCGTTTATGCCGCATCCGACAt                                         >  1:331616/1‑150 (MQ=255)
                                                                                                         tGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATACGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACAt                                         >  2:1112658/1‑150 (MQ=37)
                                                                                                         tGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATACGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACAt                                         >  2:1117841/1‑150 (MQ=37)
                                                                                                             ggTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGAGCAATAGCATGATTTAGCGTGAGTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCaaa                                     >  2:271665/1‑149 (MQ=255)
                                                                                                                cGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGcc                                  >  1:408143/1‑150 (MQ=37)
                                                                                                                cGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTAAGCTTGATTATGTAGTTCGGATAAGTCGTTTATGCAGCATACGACATCAAAGcc                                  >  2:349788/1‑150 (MQ=11)
                                                                                                                  gATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCAAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAAGCCGCATCCGACATCAACGCCTg                                >  2:170132/1‑150 (MQ=255)
                                                                                                                     atggCGTTTACGCCGCTTCAGACTTCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATg                             <  2:1149927/148‑1 (MQ=255)
                                                                                                                               cgccgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTac                   >  1:1141151/1‑149 (MQ=25)
                                                                                                                                      tCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCAGCATCCGACATCAACGCCTGATGCGAAGCTTAAAGCGTCt            >  1:166366/1‑150 (MQ=255)
                                                                                                                                          aCATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACTCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTATGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATc        <  2:21539/150‑1 (MQ=21)
                                                                                                                                           aaTCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCa       <  2:1156223/149‑1 (MQ=21)
                                                                                                                                               aaCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAAGACAGACAGCGCAATAGCCAGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGAAGCTTAACGCGTCTTATCAGGcc   >  2:537108/1‑150 (MQ=255)
                                                                                                                                                aCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCTATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCt  >  1:29347/1‑150 (MQ=14)
                                                                                                                                                 |                                                                                                                                                    
GGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCT  >  NC_000913/4295915‑4296208

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: