Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,094,361 G→A intergenic (‑89/+362) 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,3610GA98.8% 269.3 / ‑4.1 80intergenic (‑89/+362)rhaD/rhaArhamnulose‑1‑phosphate aldolase/L‑rhamnose isomerase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (31/48);  total (31/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

TGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCG  >  NC_000913/4094224‑4094495
                                                                                                                                         |                                                                                                                                      
tggCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAAt                                                                                                                                      >  2:564934/1‑140 (MQ=255)
tggCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAAt                                                                                                                                      >  2:664524/1‑140 (MQ=255)
      tGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACg                                                                                                                                <  2:23818/140‑1 (MQ=255)
       gATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGc                                                                                                                               >  1:860132/1‑140 (MQ=255)
          cATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGt                                                                                                                            <  2:956360/140‑1 (MQ=255)
          cATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGt                                                                                                                            <  1:1289246/140‑1 (MQ=255)
          cATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGt                                                                                                                            <  2:92382/140‑1 (MQ=255)
          cATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGt                                                                                                                            <  2:866730/140‑1 (MQ=255)
              cccTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTcgcg                                                                                                                        <  1:1002479/140‑1 (MQ=255)
               ccTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGt                                                                                                                       <  1:395455/140‑1 (MQ=255)
                     cAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATc                                                                                                                 <  2:860132/140‑1 (MQ=255)
                     cAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATc                                                                                                                 <  1:869161/140‑1 (MQ=255)
                     cAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATc                                                                                                                 <  1:733701/140‑1 (MQ=255)
                      aaaCCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCa                                                                                                                <  2:1203791/140‑1 (MQ=255)
                      aaaCCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCa                                                                                                                <  2:417680/140‑1 (MQ=255)
                        aCCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAgg                                                                                                              <  2:569867/140‑1 (MQ=255)
                        aCCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCACATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAgg                                                                                                              <  2:1088856/140‑1 (MQ=255)
                             gACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTAc                                                                                                         >  1:1221587/1‑140 (MQ=255)
                                  aGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTa                                                                                                    >  2:473922/1‑140 (MQ=255)
                                  aGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTa                                                                                                    >  1:467604/1‑140 (MQ=255)
                                  aGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTa                                                                                                    >  2:672064/1‑140 (MQ=255)
                                  aGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTa                                                                                                    >  1:1257672/1‑140 (MQ=255)
                                  aGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTa                                                                                                    >  1:625625/1‑140 (MQ=255)
                                  aGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTa                                                                                                    <  1:832099/140‑1 (MQ=255)
                                  aGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTa                                                                                                    >  1:800125/1‑140 (MQ=255)
                                   gTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTAt                                                                                                   <  1:1238955/140‑1 (MQ=255)
                                   gTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTAt                                                                                                   <  2:2088/140‑1 (MQ=255)
                                   gTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTAt                                                                                                   <  1:664524/140‑1 (MQ=255)
                                    tAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTAtt                                                                                                  >  2:25340/1‑140 (MQ=255)
                                    tAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTAtt                                                                                                  >  2:310043/1‑140 (MQ=255)
                                       tGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGcc                                                                                               <  1:672064/140‑1 (MQ=255)
                                       tGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGcc                                                                                               >  2:890659/1‑140 (MQ=255)
                                        gTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGTGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCg                                                                                              <  1:564934/140‑1 (MQ=255)
                                        gTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCg                                                                                              <  1:72869/140‑1 (MQ=255)
                                        gTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCg                                                                                              <  1:582643/140‑1 (MQ=255)
                                            tGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTAt                                                                                                   <  2:549119/131‑1 (MQ=255)
                                            tGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTAt                                                                                                   >  1:549119/1‑131 (MQ=255)
                                            tGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCt                                                                                          >  2:830501/1‑140 (MQ=255)
                                            tGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCt                                                                                          >  1:379115/1‑140 (MQ=255)
                                              cATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGt                                                                                        >  1:136769/1‑140 (MQ=255)
                                                 gTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAgg                                                                                     >  1:1225537/1‑140 (MQ=255)
                                                     tgtgTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTg                                                                                 >  2:59312/1‑140 (MQ=255)
                                                         ttCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATaa                                                                             <  2:800125/140‑1 (MQ=255)
                                                              gTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATaa                                                                             <  2:723003/135‑1 (MQ=255)
                                                              gTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATaa                                                                             >  1:723003/1‑135 (MQ=255)
                                                               tAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTg                                                                                         <  1:334875/122‑1 (MQ=255)
                                                               tAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTg                                                                                         >  2:334875/1‑122 (MQ=255)
                                                               tAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCt                                                                                  <  2:634175/129‑1 (MQ=255)
                                                               tAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCt                                                                                  >  1:634175/1‑129 (MQ=255)
                                                                aaaTTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACtt                                                                      <  1:890659/140‑1 (MQ=255)
                                                                 aaTTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTtg                                                                     >  2:45803/1‑140 (MQ=255)
                                                                             tCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCgg                                                         <  1:310043/140‑1 (MQ=255)
                                                                             tCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCgg                                                         <  1:25340/140‑1 (MQ=255)
                                                                             tCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCgg                                                         <  2:467604/140‑1 (MQ=255)
                                                                             tCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCgg                                                         <  2:455303/140‑1 (MQ=255)
                                                                             tCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCgg                                                         <  1:473922/140‑1 (MQ=255)
                                                                             tCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCgg                                                         >  1:1243040/1‑140 (MQ=255)
                                                                             tCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCgg                                                         <  2:1243040/140‑1 (MQ=255)
                                                                                  tGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGcaatc                                                    >  1:56908/1‑140 (MQ=255)
                                                                                  tGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGcaatc                                                    <  2:499023/140‑1 (MQ=255)
                                                                                      cGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATg                                                >  2:1290929/1‑140 (MQ=255)
                                                                                        aCCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGcc                                              <  2:1225537/140‑1 (MQ=255)
                                                                                                 ccGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGAc                                     <  2:136769/140‑1 (MQ=255)
                                                                                                   ggATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGc                                   >  2:22693/1‑140 (MQ=255)
                                                                                                    gATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCt                                  >  2:886539/1‑140 (MQ=255)
                                                                                                      tAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGt                                <  2:1112626/140‑1 (MQ=255)
                                                                                                      tAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGt                                <  2:950312/140‑1 (MQ=255)
                                                                                                      tAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGt                                <  2:228227/140‑1 (MQ=255)
                                                                                                      tAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGt                                <  1:1290929/140‑1 (MQ=255)
                                                                                                          gcgcTCGCGCCGCATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTcgcg                            <  2:56908/140‑1 (MQ=255)
                                                                                                                 cgccgcATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATc                     <  1:830501/140‑1 (MQ=18)
                                                                                                                    cgcATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAgg                  <  1:59312/140‑1 (MQ=17)
                                                                                                                     gcATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGt                                <  2:1082829/125‑1 (MQ=17)
                                                                                                                     gcATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGt                                >  1:1082829/1‑125 (MQ=17)
                                                                                                                     gcATCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGc                 >  2:388208/1‑140 (MQ=17)
                                                                                                                       aTCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCt               <  1:22693/140‑1 (MQ=17)
                                                                                                                        tCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTa              <  2:1257672/140‑1 (MQ=17)
                                                                                                                        tCCGGCAATCAATGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTa              >  1:461578/1‑140 (MQ=17)
                                                                                                                                 caatGCCTAATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTg     <  2:163165/140‑1 (MQ=17)
                                                                                                                                    tGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCg  <  1:1157614/140‑1 (MQ=32)
                                                                                                                                         |                                                                                                                                      
TGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCG  >  NC_000913/4094224‑4094495

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: