Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,296,190 A→G 11.0% intergenic (+396/+246) gltP → / ← yjcO glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0009134,296,1900A→G11.0% 3.0 ‑0.3 0.059intergenic (+396/+246)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  new base (3/4):  ref base (26/24):  total (31/28)
Fisher's exact test strand distribution p-value = 7.06e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACC  >  NC_000913/4295892‑4296482
                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                      
gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGCTTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAATGCCCGATGCGACGCTTaa                                                                                                                                                                                                                                                                                                       >  2:459731/1‑300 (MQ=25)
 ttttAGCGTGATTTTGTAGTTGGGATAAGGCGTTTCTGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAAcg                                                                                                                                                                                                                                                                                                     <  2:473689/300‑1 (MQ=25)
  tttGGCGTCGTTTTGTAGGTCGGATAAGCCCGGTATTCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCTTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTGGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAAcgc                                                                                                                                                                                                                                                                                                    <  2:255780/301‑1 (MQ=255)
      gCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATCGCCTGCTTTAGCGTGATTTTGTAGGTCGGCTAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGTGACTCTTACAGCGt                                                                                                                                                                                                                                                                                                  >  2:595584/1‑299 (MQ=255)
           aTTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTAtga                                                                                                                                                                                                                                                                                           >  2:509719/1‑299 (MQ=25)
              ttGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCCGGCCTCCGCCAGACAGCGCACGAGCCTGATTTACCACGACTTTGTAGGTCGGCTAAGCAGTTTTTGCCGCCTCCCCCACGAACCCCCTATGCTCCGCTTTACGCGACTTaacccga                                                                                                                                                                                                                                                                                        >  1:806136/1‑295 (MQ=11)
              ttGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACCGCGACAACGCTTGATTTCGCTTGATTTTGGCGATCGGATCCGGCCTTTCTGCTCCATTCCAAACCAAACCCTCATGCCCCGCTTacccctccgtatccgt                                                                                                                                                                                                                                                                                         >  1:831113/1‑285 (MQ=11)
                  aGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACACCAGACAGCGCAATAGCCTGATTCAGCGTGATTTTGTGAGTCGGCTAAGGCGCTTATGCCACATCCGACATCAACTCCTGATCCGACGCATAACGCGTCttcttcgacctc                                                                                                                                                                                                                                                                                     >  1:550229/1‑290 (MQ=11)
                   ggTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGAGTCTTATCAGGCCTAc                                                                                                                                                                                                                                                                                    >  2:13927/1‑300 (MQ=25)
                      cGGAGAAGGCGTGTATGCCGCATCCGACATCAACGCCTGAGGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGcc                                                                                                                                                                                                                                                                                 <  2:189617/300‑1 (MQ=25)
                         aaacgtcGTTTATCCCGCATCCGAAATAAACGCCTGTGGCGACGCTTAACGGTTCTTAGCAGGCCTACGCCAGACAGTGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCGGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGAc                                                                                                                                                                                                                                                                             <  2:497781/295‑1 (MQ=11)
                         aTAAGGCGTTTATGCCGCGTCCGACATCAACGCCTGATGCGACGCTGAACGCGTCTGAGCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTAATTTTGTAGGGCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGAc                                                                                                                                                                                                                                                                             <  2:269840/301‑1 (MQ=255)
                                 tttcgtccgctgccGACACGAACGCCTTATGCGACTCTTAGCGCTCTTTATCAGGCCTACGCCAGACAGCGCAAGAGCCTTATTTATCGTTATTTGGTAGGTCGGATAAGGGCTTGACGCCGGATCCTACATCAACGTCTGATGCGACGCTTGACGCGTCTTATCAGGCTTACTCCAGCCAGCGAAATAGCCTGATTTATCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGTCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAgcgc                                                                                                                                                                                                                                                                        <  2:628444/286‑1 (MQ=11)
                                   tatgcccactCCGCCAGCCACGACTGATGCGACGCTTAACGCGTCTTTTGAGGCCTACGCCAGACAGCGCAATTGCCTGATTTCGCGTGATTTGGTAGGTCGGATAAGGCGTTTACGCCCCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTGTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAg                                                                                                                                                                                                                                                                   <  2:305154/292‑1 (MQ=11)
                                     tGCCGCACCCGCCATCAACGCCTTATGCGACGCTTAACGCGTCTTATCATGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATGTTGTAGGTCTGATAAGGCGTTTACGCCGCATCCGACATCACCGCCGGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGcc                                                                                                                                                                                                                                                                 <  2:946137/301‑1 (MQ=255)
                                      ggggctcggggatgtgtataagagccagACGCCTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGGGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGcc                                                                                                                                                                                                                                                                 <  2:418228/273‑1 (MQ=11)
                                           aacgagaacgaaCGCCTGATGAGAAGCTTACGGCGGATGGTCGGGCACACGCCAGGCGGGGCAATGGCCTGATGTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGAttt                                                                                                                                                                                                                                                           <  1:459717/291‑1 (MQ=11)
                                           aTCAGACATCAACGCCTGATGCAACGCTTAACGCGTATTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTGAGCGAGATTGTGTAGGTCGGATAAGGCGGTTATGCCTGATGCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGAttt                                                                                                                                                                                                                                                           <  2:690055/301‑1 (MQ=11)
                                            tCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGCTGCGACGCTTGCCGCGTGTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACACCGCATCCGACATCAACTCCTGCTGCGACGCTTAACGCGTCTTATCAGCCCTACGCCCGACACCCCACTACCCTGAttca                                                                                                                                                                                                                                                          >  1:215062/1‑299 (MQ=11)
                                               aataTCAACGCCTGATGCGACGCTGAACGCGTCTTAGCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCGGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCg                                                                                                                                                                                                                                                       <  2:519121/298‑1 (MQ=21)
                                                                            cgTCTTATCAGGCTTACGCCAGACAGCGCAATAGCCGGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTac                                                                                                                                                                                                                          <  2:215069/301‑2 (MQ=17)
                                                                                      ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGTTTTACCGTGCTTTTGTAGCTCGGATACGGCGTTGAGTCCGCATCCGACATCAACGTCTGATGCGACGCTTACCGAGTCCTCTCAGGCCTAAGCCTGCCCTCGCAATCGCATGCTTTACCGTGTTTCTGTCGTTTGGATAAGGCTtttcactcaatccc                                                                                                                                                                                                                 >  1:761330/1‑289 (MQ=11)
                                                                                      ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGCTAAGCCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACCCCAGACCGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGCTAAGGCGTTTACGCCGCATCCg                                                                                                                                                                                                                >  2:251371/1‑301 (MQ=255)
                                                                                      ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCCGACCGCGCAATAGCCTGATTTAGCGTGATTTTGTCGGTCGGATACGGCGTTTACGCCCCATCCg                                                                                                                                                                                                                >  2:747296/1‑301 (MQ=19)
                                                                                      ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATcc                                                                                                                                                                                                                 >  2:400041/1‑300 (MQ=25)
                                                                                           aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACGGCGCAATAGCCTGATTTAGCGTGATTTTGTCGGTCGGATAAGGCGTTTATGCCGCCTCCGACATCAACGCCTGATTCGACGCTTAACGCGTCTTATCAGTCCTCAGCACGACCGCGCAATAGCCTGATTTAGCGTTGTTTTGTAGTCCAGCTACGGCGTTTaccacccatccgccca                                                                                                                                                                                                            >  1:649248/1‑285 (MQ=11)
                                                                                           aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGTCGTTTAGGCCGCATCCGCCATCACAGCCTGATCCGACGCTTCACGCGTCTTATCAGTCCTACCCCAGCCAGCGCACTATCCTGCTTTCTCGTCATTCTGTGTGTCAGAtcctgcttttacgccccactcactcaa                                                                                                                                                                                                           >  1:40543/1‑275 (MQ=11)
                                                                                           aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCAAGACAGCGCAATAGCCGGATTTAGCGTGATTTGGTAGGGCGGATAAGGCGTTTACGCCGCATTCTACAt                                                                                                                                                                                                            >  2:420870/1‑300 (MQ=255)
                                                                                           aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTCTCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACAt                                                                                                                                                                                                            >  2:412363/1‑300 (MQ=25)
                                                                                           aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTCCGCCAGACAGCGCCATCGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGcagataccgactt                                                                                                                                                                                                            >  2:841319/1‑288 (MQ=255)
                                                                                           aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAGTCCTCATACGACATCAACGCCTGATGCGCCGCTTAACGCGTGTTGTAAGGCCTATGCCAGACAGCGGTATAGACGGCATTCGCGTCATTATGTCGGTCAGAGACGGGTTTTACGCGCCATCCGCCAt                                                                                                                                                                                                            >  1:193054/1‑300 (MQ=11)
                                                                                           aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGGTTACGCCGCCTCCGACATCAACGCCTGCTGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAAcgc                                                                                                                                                                                                                                                                                                    <  2:903517/212‑1 (MQ=21)
                                                                                               gggtaagatgtgtatgagtgtcgtgCGTTATTTGGTAGGTAGGAGAAGGCGTGTCGGCTGCGTTGTACATGACCGGGTGTGGCGGCGGTTGGCTCGCGTTATGAGGGCTACGCCAGACAGAGCAATAGACTGTTTTGGCGGGATTTTGTAGGGCGGGTAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCaat                                                                                                                                                                                                        <  1:835773/276‑2 (MQ=2)
                                                                                                aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATACGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTCTCAGGCCTACGCCAGCCAGCGCCATCGCCTGATTTAGCGTGATTTTGTCGGTCGGCTCAGGCTTTTATGCCGCATCCGACATCCCCGCCTGATTCGCCTCTTACCGCGTCTTATCCGGACTCCTCCCGCCACCGACATAGCTTCCTCTAGCTTGATTTTGCACGTCGGCACAGGCGTATaccgcgacttccccatca                                                                                                                                                                                                          >  1:279302/1‑280 (MQ=2)
                                                                                                aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTATGTCGGATAAGGCGTTTATGCCGCATCCGACATCTACGCCTGATGCGCCGCTTGACGCGTCTTATCAGCCCTACGCCATACCGCGCCATAGCCTGACTTCGCTTGATTTTGTAGATCTGATAAGACGTTTacaccccctccgccaactatc                                                                                                                                                                                                       >  1:350851/1‑280 (MQ=11)
                                                                                                  aCAGCGCAATAGCCGAATTTCGCGGGATGTTGTAGGTCGGAGAGGCGGTTTACGCCGTATCCGCCATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGAAATAGCCTGAGTGAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCaatg                                                                                                                                                                                                       <  2:698309/298‑3 (MQ=11)
                                                                                                     gcgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGACTGATTTAGCGTGATGTTGTAGGTCGGATAAGGCGATTACGCCGCATCATACATAAATGCATGAt                                                                                                                                                                                                 >  2:180813/1‑301 (MQ=255)
                                                                                                                 gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCGGCGTCTTAGCAGGCCTACGCCAGACAGAGCATTAGCCTGATTTAGGGTGTGTTTGTGGGTGGGATAAGGCGTTTATGCCGCATCCAACATCAAGGCCGGATGCGACGCTTAACGCGTCTTGTCAGGCCTGCGCCCACCGGCGCAATAGCCGAATTCAGCGTGCTTTTGTAGTGCGGCTACGGCGGTTAGACCTCCTCCGCCATCTATGCCTGATGCgcagagtgcg                                                                                                                                                                                     >  2:859703/1‑292 (MQ=11)
                                                                                                                 gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGGCGGATACGGCGTTTGCGCGCCATCCGCCTTCAATCCCTGTTggtcgcattgc                                                                                                                                                                                      >  2:378154/1‑290 (MQ=11)
                                                                                                                       gCGTGATTTTGTATGTCGGATAAGGCGTTTACGCCGCAGCCTACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTc                                                                                                                                                                                <  2:550245/300‑1 (MQ=12)
                                                                                                                       gCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCACTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGGCATCact                                                                                                                                                                                                        >  2:835801/1‑274 (MQ=25)
                                                                                                                       gCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATCGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGAATCCGACATAAATGCCTGATGCGAccattgccgcgt                                                                                                                                                                                 >  2:465221/1‑288 (MQ=255)
                                                                                                                           taTTTTGTTGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCTGGATGCGACGCTTGCCGAGTCTGATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATc                                                                                                                                                                           <  2:806165/300‑1 (MQ=255)
                                                                                                                            aTTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTCGCGTGATTTTGTCGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATc                                                                                                                                                                           >  2:801289/1‑300 (MQ=12)
                                                                                                                            aTTTTGTAGGTAGGAAAAGGCTTTTACGCCGCATCCGTCATCAACGCCTGAGGCGACTGTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGGCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATc                                                                                                                                                                           <  2:326542/300‑1 (MQ=11)
                                                                                                                                tggaaaaagggacaGGCGGTTCCCCCGCATCCGACATCAACGCCTGATGGGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGAGAAGGCGTTGACGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGAGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                                                                      <  1:668888/288‑1 (MQ=11)
                                                                                                                                tGTAGGGCGGATAAGGCGTTTATGGCGCATCCGACATACACGCCGGATGCGACGCTTGCCCGGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                                                                      <  2:345578/301‑1 (MQ=11)
                                                                                                                                aagagcagtgATAATGGGTTTACGCCGCATCCGACATTAACGGTGGATGCGCGGCTGGCCGGGCTTTAGCAGGCCAACGCCAGACAGCGCAATAGCCTGAGTTAGCGTGATTTTGGAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGTCGCTTGCCGCGTCTTATCAGGcc                                                                                                                                                                      <  1:32214/292‑1 (MQ=11)
                                                                                                                                        ggATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTACCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCCTTCCACATCAATGCCTGCGGCGACCCTCTCCGCTTCTTATCAGCCCTAcattaa                                                                                                                                                              >  1:604458/1‑296 (MQ=11)
                                                                                                                                         gATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGTCGCTTGCCGCGTCTTATCAGGcc                                                                                                                                                                      >  2:32214/1‑292 (MQ=17)
                                                                                                                                               gCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGCCGCTTGCCGCGTCTTATCAGGcc                                                                                                                                                                      >  2:668916/1‑286 (MQ=12)
                                                                                                                                                             tCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTGTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATc                                                                                                                                                                           <  2:709863/267‑1 (MQ=9)
                                                                                                                                                                             gATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTACGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACTCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATCGCCGGATTTACCGTCATTTTTTAGCTCGGATAAGGCGGACATGCCCCATCCAGCATCAACCCCTGCTTCGCACCGTCACGCCGATTATCCTGCATACGCTCGACCGCGCAATCGccacctttcaccttatttactagat                                                                                                                         >  1:690027/1‑278 (MQ=2)
                                                                                                                                                                                                                                                          gATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTAGGCCAGCCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTaa                                             >  2:712284/1‑300 (MQ=255)
                                                                                                                                                                                                                                                            aaaGGAGTTTACGCCGCATCCGAAATCAACGCCGGAGGCGACGCTGGCAGCGTCTGATCAGGCCTACGCCAGACATCGCAATAGCCTGAGTGAGCGTGATTTTGTAGGTCGGATAAGTCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCg                                          <  2:279309/300‑1 (MQ=11)
                                                                                                                                                                                                                                                             gcggcaggtgagagggttgtaaagggaaagCCTGATGCGACGCTTGCGGCGTCTTAGCAGGCGTACGCCAGACATCGAAATAGCCTGATGTAGCGTGATTTGGTAGGTCGGATAAGGCGGTGATGCCGCGTCCGACATAAACGCCTGATGCGACGCTTGACGCGTCTGAGCAGGCCTACGCCAGACAGCGCAATAGCCTGATGTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCg                                          <  1:691143/271‑1 (MQ=11)
                                                                                                                                                                                                                                                                gCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCTGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCCTCCGACCTCCATTCCTGATGCGACGCTTGCCGCGTCTTATCCGGCCTATCTTACCCGTTgg                                      >  2:686699/1‑301 (MQ=11)
                                                                                                                                                                                                                                                                   tttACGCCGCATCCGCCAACGTCGCCGGATGCAACGCTTACCGCGGCTTAGCGGGCCTACGCCAGACAGCGCAATGGCCTGATTTAGCGTAATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGtt                                    <  2:761358/300‑1 (MQ=11)
                                                                                                                                                                                                                                                                                              gATGGGGCGCTTGCCGCGTCTTTGCAGGCCTACGCCAGACAGCGCAATGGCCTGATTGAGCGGGATGTGTTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTa         <  2:235238/300‑1 (MQ=14)
                                                                                                                                                                                                                                                                                               aTGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGTTTAATTTTTCAAACTCTTCACACCCGGTAt        >  2:671934/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                    aCGCTTGCCCCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAcc  <  2:649276/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                      
GATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACC  >  NC_000913/4295892‑4296482

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.