Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,296,060 C→T 38.2% 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,0600CT38.2% 55.2 / 78.5 76intergenic (+266/+376)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base C (24/23);  new base T (18/11);  total (42/34)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

GTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCT  >  NC_000913/4295912‑4296208
                                                                                                                                                    |                                                                                                                                                    
gTCGGATCAGGCGTTTATGCCGCTTCAGTCATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAc                                                                                                                                                      <  2:90493/149‑1 (MQ=255)
    gATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCt                                                                                                                                                  <  1:238649/149‑1 (MQ=17)
      tAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGa                                                                                                                                                <  6:348679/149‑1 (MQ=32)
         ggCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCATGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGc                                                                                                                                             >  6:244723/1‑149 (MQ=17)
          gCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGc                                                                                                                                             >  3:152614/1‑148 (MQ=17)
           cGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGa                                                                                                                                           >  3:38638/1‑149 (MQ=17)
            gTTTATGCCGCAGCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGAc                                                                                                                                          <  4:183166/149‑1 (MQ=17)
                 tGCCGCATCCGAAATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTg                                                                                                                                     <  8:257810/149‑1 (MQ=17)
                   ccGCATCCGTCATCAAGGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGc                                                                                                                                    <  5:266428/148‑1 (MQ=18)
                   ccGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGcc                                                                                                                                   >  5:233327/1‑149 (MQ=21)
                   ccGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGcc                                                                                                                                   >  2:79317/1‑149 (MQ=21)
                     gCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCcgc                                                                                                                                 <  1:79317/149‑1 (MQ=25)
                     gCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCcgc                                                                                                                                 >  4:113140/1‑149 (MQ=25)
                        tCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTc                                                                                                                              >  8:240696/1‑149 (MQ=37)
                          cGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCtt                                                                                                                            <  7:221965/149‑1 (MQ=25)
                            aCATCAACGACTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCCTGATTTTGTAGGTCGGATAAGGCGTAAACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCCTCTTAt                                                                                                                          >  3:227709/1‑149 (MQ=11)
                             cttCAACGCCTGATGCGTCGATTAATGTGTATTATCAGGCCTACGCCTTACAGCGTAATAGCATGATTTAGCTTGATTTTGTTTGTCGGCTATGGGGTTTACGCCGGTTCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                                                          <  2:227752/146‑1 (MQ=11)
                              aTCAACGTGTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCTCAATAGCCTGATTTAGCGTGATTTTGTTGTTCGGATAAGGCTTTTACGCCGCTTCAGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATc                                                                                                                         <  8:375658/148‑1 (MQ=11)
                              aTCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                                                          >  3:194338/1‑147 (MQ=25)
                              aTCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                                                          >  7:124242/1‑147 (MQ=25)
                              aTCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCa                                                                                                                        <  7:240696/149‑1 (MQ=37)
                               tCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACGGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAg                                                                                                                       <  8:166557/149‑1 (MQ=255)
                                 aaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACTACAGACAGCGCAATAGACTGATTTCGCGTGATTTTGGAGTTAGGATAAGGAGTTTACGCGGCATCCGACAGCAGCGCCGTATGCGAAGCTTGGAGAGTCTTTTCAGGc                                                                                                                     >  7:272075/1‑149 (MQ=11)
                                  aCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                    >  8:263993/1‑149 (MQ=25)
                                   cGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                    >  1:349688/1‑148 (MQ=25)
                                     ccTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                                                                                                                 <  5:86702/149‑2 (MQ=21)
                                      ccgacgCGACGCGTAACGCGTCTTCTCAGGCCTACGCCAGACCGCGCAATCGCCTGATTTAGCGGGATTTTGTCGGGCGGATAAGGCGGTTACGCCGCATCCGACATCAATGCCTGATGCGCCGCTTGCCGCGTCTTATCAGGCCTatc                                                                                                                <  5:341615/144‑3 (MQ=11)
                                        gATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatct                                                                                                               >  1:247318/1‑145 (MQ=17)
                                           gCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGAGCTATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGcgcgc                                                                                                           >  3:6697/1‑145 (MQ=255)
                                            cGACGCTTCCCGCGTCTATTCCGGCATACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTCGGGCGGATATGGCTTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGAc                                                                                                          <  7:338463/149‑1 (MQ=11)
                                            cGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGAc                                                                                                          >  5:170622/1‑149 (MQ=37)
                                              aCGCTTAACGCGTCTTATCAGTCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGCTCGTATACGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAg                                                                                                        <  2:15706/149‑1 (MQ=255)
                                                gCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCTGATAAGGCGTTTAGGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAgcg                                                                                                      <  7:95307/149‑1 (MQ=21)
                                                      cgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAt                                                                                                                                               <  1:156193/102‑1 (MQ=32)
                                                      cgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCGGAt                                                                                                                                               >  2:156193/1‑102 (MQ=17)
                                                       gcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGc                                                                                               <  6:29274/149‑1 (MQ=37)
                                                          tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTc                                                                                                                              >  4:292226/1‑115 (MQ=39)
                                                          tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTc                                                                                                                              <  3:292226/115‑1 (MQ=39)
                                                            ttATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                     >  2:304461/1‑122 (MQ=25)
                                                            ttATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                     <  1:304461/122‑1 (MQ=25)
                                                                cAGGCCTACGCCAGACACCTCACTAGCCTGATTTAGCTTGATTTTGTACTGCTGCTATGGCGTTTACGCCGCCTAAGATATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTa                                                                                       <  8:108555/148‑1 (MQ=11)
                                                                  ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGAc                                                                                                                                          >  6:40428/1‑95 (MQ=18)
                                                                  ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGAc                                                                                                                                          <  5:40428/95‑1 (MQ=18)
                                                                  ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGa                                                                                                                                                <  1:141761/89‑1 (MQ=33)
                                                                  ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCcg                                                                                                                                  >  7:32746/1‑103 (MQ=38)
                                                                  ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCcg                                                                                                                                  <  8:32746/103‑1 (MQ=38)
                                                                      tACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACg                                                                                                                                         >  6:51677/1‑92 (MQ=18)
                                                                                  cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                                                          >  3:331473/1‑95 (MQ=255)
                                                                                  cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                                                          <  4:331473/95‑1 (MQ=255)
                                                                                     aaTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGcc                                                                                                                                                   <  3:268537/67‑1 (MQ=25)
                                                                                     aaTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGcc                                                                                                                                                   >  4:268537/1‑67 (MQ=25)
                                                                                          ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGt                                                                                   >  7:202621/1‑126 (MQ=17)
                                                                                                 tAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGc                                                     >  7:19668/1‑149 (MQ=37)
                                                                                                   gCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGAt                                                                                                                                               >  5:296424/1‑57 (MQ=255)
                                                                                                   gCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGAt                                                                                                                                               <  6:296424/57‑1 (MQ=255)
                                                                                                    cGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACGCCAGACAGCGCAAt                                                                                                  >  7:313221/1‑101 (MQ=17)
                                                                                                    cGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGcc                                                                                              >  6:237152/1‑105 (MQ=17)
                                                                                                     gTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCa                                                 >  7:210905/1‑149 (MQ=21)
                                                                                                      tGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCAt                                                >  3:4626/1‑149 (MQ=37)
                                                                                                              tAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACAt                                         <  4:4626/148‑1 (MQ=37)
                                                                                                                  ccgatctaGGCGTTTACGCCGCATCCGACAGCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACg                                    <  3:272750/142‑1 (MQ=255)
                                                                                                                   cGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGAc                                           <  3:107780/141‑1 (MQ=37)
                                                                                                                   cGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGAc                                           >  4:107780/1‑141 (MQ=37)
                                                                                                                    ggATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTa                                                                                                                           >  5:247996/1‑60 (MQ=255)
                                                                                                                     gATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAAc                                     >  2:124841/1‑145 (MQ=37)
                                                                                                                     gATAAGGCGTTTACGCCGCATCCGACAGCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAAc                                     <  1:124841/145‑1 (MQ=37)
                                                                                                                          ggCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGAGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGGCTTCACCg                                    >  4:272750/1‑141 (MQ=255)
                                                                                                                             gTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGa                                                                                           >  8:47307/1‑83 (MQ=255)
                                                                                                                                   gccgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGAAAGACAGCGCAATAGCCGGATTTAGCGTGATTTTGTAGTTCGTATAAGGCGTTTAAGCCGCAACCGACATCAAAGACTGATGCGACGCTTac                   >  6:335801/1‑148 (MQ=11)
                                                                                                                                     cgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTc                                                                      <  5:214343/96‑1 (MQ=37)
                                                                                                                                     cgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTc                                                                      >  6:214343/1‑96 (MQ=37)
                                                                                                                                     cgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAAcg                 >  7:240692/1‑149 (MQ=34)
                                                                                                                                            gACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGAGTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTa          >  3:183166/1‑149 (MQ=21)
                                                                                                                                            gACATCAACGCCTCATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCTCAATATCCTGATTTAGCGTGATTTTTTATGTATGATAAGGCGTTTATGCCTCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTa          <  4:6697/149‑1 (MQ=11)
                                                                                                                                                 cAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGGTTATCGCGTCTTATCAgg     >  5:74252/1‑149 (MQ=17)
                                                                                                                                                  aaCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTaa                   <  8:318788/134‑1 (MQ=32)
                                                                                                                                                  aaCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTaa                   >  7:318788/1‑134 (MQ=32)
                                                                                                                                                    cGCCTTACGCTACGCTTGCCGCGTCTTATGAGGCCTACGCGAGACAGCGCAATATCCTGACTTAGCGTGATTTTGTAGGTCTGCTAAGGCGTTTATTCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCt  <  3:40313/149‑1 (MQ=11)
                                                                                                                                                    |                                                                                                                                                    
GTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCT  >  NC_000913/4295912‑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: