Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,094,489 G→A 7.0% intergenic (‑217/+234) rhaD ← / ← rhaA rhamnulose‑1‑phosphate aldolase/L‑rhamnose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,094,4890GA 7.0% 245.6 / 7.1 88intergenic (‑217/+234)rhaD/rhaArhamnulose‑1‑phosphate aldolase/L‑rhamnose isomerase
Reads supporting (aligned to +/- strand):  ref base G (30/49);  new base A (4/2);  total (36/52)
Fisher's exact test for biased strand distribution p-value = 2.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.97e-01

ATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCC  >  NC_000913/4094343‑4094635
                                                                                                                                                  |                                                                                                                                                  
aTCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGtt                                                                                                                                                  <  2:306791/149‑1 (MQ=32)
 tCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTg                                                                                                                                                 <  7:174610/149‑1 (MQ=32)
   cGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGc                                                                                                                                                <  7:80722/148‑1 (MQ=32)
   cGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGc                                                                                                                                                <  5:11998/148‑1 (MQ=32)
    ggCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCg                                                                                                                                              <  2:96523/149‑1 (MQ=32)
      caatcaatGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGAc                                                                                                                                            <  8:270899/149‑1 (MQ=32)
       aatcaatGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGACTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTCTTGCCGAcc                                                                                                                                           >  3:282533/1‑149 (MQ=17)
        atcaatGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCt                                                                                                                                          <  6:105623/149‑1 (MQ=32)
         tcaatGCCTGATGCGACGCTGTCGCGTCTTATCAGGCGTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTAGCAGGCGTACAAGTGTTGCCGACATg                                                                                                                                         >  5:252295/1‑149 (MQ=11)
             tGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAgg                                                                                                                                     <  3:303829/149‑1 (MQ=34)
                cTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCg                                                                                                                                  <  6:224127/149‑1 (MQ=34)
                cTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCg                                                                                                                                  <  5:335179/149‑1 (MQ=34)
                   aTGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGCAGGCCGGAt                                                                                                                               <  6:365145/149‑1 (MQ=21)
                   aTGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGa                                                                                                                                >  3:270005/1‑148 (MQ=17)
                      cGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCACTGCCTGATGCGACGCTGGCGCGTCTAATCCGGCCTACAACTGTTGCCGACCTGTAGGCCGGATaa                                                                                                                             <  5:295678/148‑1 (MQ=255)
                      cGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATaa                                                                                                                             <  8:372919/148‑1 (MQ=34)
                      cGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAg                                                                                                                            <  7:407365/149‑1 (MQ=34)
                       gACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGGTGCCGCCCTGTAGGCCGGATAAgg                                                                                                                           >  5:94319/1‑149 (MQ=21)
                       gACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATAATGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAgg                                                                                                                           >  4:131640/1‑149 (MQ=21)
                           cTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACt                                                                                                                       >  6:367104/1‑149 (MQ=17)
                              tCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCCGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACtt                                                                                                                      <  2:154211/147‑1 (MQ=37)
                                gcgTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGc                                                                                                                  <  3:262095/149‑1 (MQ=37)
                                 cgTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGcc                                                                                                                 <  4:170187/149‑1 (MQ=37)
                                 cgTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGcc                                                                                                                 <  3:66801/149‑1 (MQ=37)
                                   tCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTCATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGc                                                                                                               <  1:76343/149‑1 (MQ=37)
                                     ttATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCAT                                                                                                             >  1:148465‑M1/1‑149 (MQ=255)
                                     ttATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCA                                                                                                              >  1:390683‑M1/1‑148 (MQ=255)
                                      tATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTCGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATC                                                                                                            >  7:361948‑M1/1‑149 (MQ=255)
                                       aTCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGcc                                                                                                                 <  6:292467/143‑1 (MQ=37)
                                       aTCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGcc                                                                                                                 >  1:280847/1‑143 (MQ=37)
                                       aTCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGcc                                                                                                                 >  5:292467/1‑143 (MQ=37)
                                       aTCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGcc                                                                                                                 <  2:280847/143‑1 (MQ=37)
                                          aGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATGCGGCACAAGTGCCTTATCCGGC                                                                                                        <  5:273690‑M1/149‑1 (MQ=255)
                                           ggCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATGCGGCACAAGTGCCTTATCCGGCc                                                                                                       <  7:170586‑M1/149‑2 (MQ=255)
                                              cTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGcaat                                                                                                     >  2:6592/1‑148 (MQ=37)
                                                aCAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCTAGTCTTATCAGGCCTAAAACTGTTGACGAACTGTAGGCCGGATAAGGCACTTGTGCCGCGGCCGGcaatcaa                                                                                                  >  7:326574/1‑149 (MQ=11)
                                                 ctaCTATTGCCGACCTGTAGGCCTGATAAGGCTCTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGcaatcaat                                                                                                 <  5:311289/147‑1 (MQ=37)
                                                 cAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGcaatcaat                                                                                                 <  5:176639/149‑1 (MQ=37)
                                                   aCTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGc                                                                                               <  8:40385/149‑1 (MQ=37)
                                                     tATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGcaat                                                                                                     >  6:168530/1‑141 (MQ=35)
                                                     tATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGcaat                                                                                                     <  5:168530/141‑1 (MQ=35)
                                                     tATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCt                                                                                             >  3:2244/1‑149 (MQ=35)
                                                     tATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCt                                                                                             >  6:134114/1‑149 (MQ=35)
                                                          ccGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATg                                                                                         <  2:187030/148‑1 (MQ=34)
                                                          ccGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGc                                                                                        <  4:196421/149‑1 (MQ=34)
                                                             aCCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATGCGGCACAAGTGCCTTATCCGGCctacaggtcggcaacagtt                                                                                     <  7:84622‑M1/149‑20 (MQ=255)
                                                                tGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCt                                                                                  >  6:358508/1‑149 (MQ=34)
                                                                tGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCt                                                                                  >  5:84432/1‑149 (MQ=34)
                                                                    ggCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTTTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACTCTGTcg                                                                              >  2:305304/1‑149 (MQ=17)
                                                                     gCCTGATAAGGCACTTGTGCCTCATCGGGCAATCAATGCCTGTTGTGACGCTTTCGCGTCTTTTCAGGCCTACAACTGTTGCCTACATGTTGTCCGGATATGGCACTTGTGCAGCATCCGGCAATCAATGCCTGATGCGACTCTGTcgc                                                                             <  4:282533/149‑1 (MQ=11)
                                                                     gCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTcgc                                                                             <  6:222006/149‑1 (MQ=34)
                                                                       cTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGt                                                                           <  2:42536/149‑1 (MQ=33)
                                                                          aTAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTc                                                                          <  2:210252/147‑1 (MQ=32)
                                                                           tAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCtt                                                                        <  7:64305/148‑1 (MQ=32)
                                                                              ggCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCa                                                                    <  2:148465/149‑1 (MQ=32)
                                                                                cACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAgg                                                                  <  8:27053/149‑1 (MQ=32)
                                                                                      tgCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACa                                                            <  1:6592/149‑1 (MQ=32)
                                                                                      tgCCGCATCCGGCAATCAATGCCGGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTAc                                                             <  6:94319/148‑1 (MQ=17)
                                                                                       gCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACaa                                                           <  6:308169/149‑1 (MQ=32)
                                                                                         cGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAc                                                          >  6:230337/1‑148 (MQ=32)
                                                                                               cGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGc                                                    <  4:87436/148‑1 (MQ=34)
                                                                                               cGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTAAAAATGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAATATTGcc                                                   >  8:227951/1‑149 (MQ=255)
                                                                                                ggCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATc                                                                     >  1:198764/1‑130 (MQ=32)
                                                                                                ggCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATc                                                                     <  2:198764/130‑1 (MQ=32)
                                                                                                   aatcaatGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGAcc                                               <  3:131640/149‑1 (MQ=34)
                                                                                                   aatcaatGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATGCGGCACAAGTGCCTTATCAGGCctacaggtcgg                                                                                             <  1:282726‑M1/103‑12 (MQ=255)
                                                                                                   aatcaatGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATGCGGCACAAGTGCCTTATCAGGCctacaggtcgg                                                                                             >  2:282726‑M1/1‑92 (MQ=255)
                                                                                                    atcaatGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCt                                              >  2:121796/1‑149 (MQ=35)
                                                                                                     tcaatGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTg                                             <  5:134114/149‑1 (MQ=35)
                                                                                                           ccTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGcc                                       >  8:53063/1‑149 (MQ=37)
                                                                                                                 gCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGAGCCGCTAACGGCAGTCCATGCCTTATGCGAAGCTGTAGCTTCTTATCAGTCCTACAAGTATTGGCGCGATGTAGTACGGAt                                   >  6:219780/1‑147 (MQ=11)
                                                                                                                  cGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATaa                                 <  4:270005/148‑1 (MQ=17)
                                                                                                                   gACGCTGTCTCGTCTTATCAGGCCTACTACTTTTTCAGACCTTTAGGCCGGATATGGCACTTGTGCCGCATCCGTCAATCAATGCCTGATGCTACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAgg                               <  4:329047/149‑1 (MQ=11)
                                                                                                                        tGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAACAATTGCCGACCTGTAGGCGGGATAAGGCAAtt                          >  7:61939/1‑149 (MQ=255)
                                                                                                                            gcgTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGc                                                    >  7:36388/1‑119 (MQ=35)
                                                                                                                            gcgTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGc                                                    <  8:36388/119‑1 (MQ=35)
                                                                                                                            gcgTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGTTAAGGCACTTGTGc                      >  1:46856/1‑149 (MQ=37)
                                                                                                                                 ttATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCAt                 >  5:142769/1‑149 (MQ=37)
                                                                                                                                    tCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCg              >  8:16877/1‑149 (MQ=37)
                                                                                                                                        gCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGc                      <  6:333892/137‑1 (MQ=37)
                                                                                                                                        gCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGc                      >  5:333892/1‑137 (MQ=37)
                                                                                                                                         ccTACAACTGTTGCCGACATGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCATTAAATGCCTGATGCGACGCTGTCGAGTATTATATGGCCTACAACAAGTCCAGACCTGTAGGCCGGAAAAGGCACATGAGCCGCATCCGGgaat         >  6:164785/1‑149 (MQ=11)
                                                                                                                                              aaCTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATg    <  6:84432/149‑1 (MQ=37)
                                                                                                                                               aCTGTTGCCGACCTGTAGGGCGGATAAGGCACTTGTGGCGCTTCCGGCAATCAATGCCTGATGGGTCGCTGTCGCGTCTTATCAGGCGTACTACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGc   <  8:61939/149‑1 (MQ=11)
                                                                                                                                               aCTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGc   >  2:244885/1‑149 (MQ=37)
                                                                                                                                                cTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGcc  >  7:148660/1‑149 (MQ=37)
                                                                                                                                                 tGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAg                                                                   >  2:273902/1‑83 (MQ=32)
                                                                                                                                                 tGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAg                                                                   <  1:273902/83‑1 (MQ=32)
                                                                                                                                                  |                                                                                                                                                  
ATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCC  >  NC_000913/4094343‑4094635

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: