Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,094,414 T→G 5.4% intergenic (‑143/+308) 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,4140TG 5.4% 250.5 / 2.5 100intergenic (‑143/+308)rhaD/rhaArhamnulose‑1‑phosphate aldolase/L‑rhamnose isomerase
Reads supporting (aligned to +/- strand):  ref base T (52/39);  new base G (5/4);  total (57/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 101 total aligned reads displayed.

CATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCG  >  NC_000913/4094269‑4094560
                                                                                                                                                 |                                                                                                                                                  
cATGTTCTGTGTTACTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGAGCTGTCAGCCTGAt                                                                                                                                                 >  1:825579/1‑149 (MQ=255)
   gTTCTTTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAg                                                                                                                                              <  1:1038648/149‑1 (MQ=255)
   gTTCTGTGTGCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAg                                                                                                                                              <  2:872266/149‑1 (MQ=255)
    ttCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGGAGAGTTGTAGGCGTGATTAgg                                                                                                                                             >  2:748948/1‑149 (MQ=21)
    ttCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAgg                                                                                                                                             <  1:224487/149‑1 (MQ=255)
    ttCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCGGTAGGCCTGATAAgg                                                                                                                                             >  2:943131/1‑149 (MQ=255)
    ttCTGTGTTCCTGTAAATTCGGTGTTGTAGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGCTACgg                                                                                                                                             >  1:371346/1‑149 (MQ=255)
     tgtgggtTCCTGTAAAGTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGc                                                                                                                                            <  2:521824/144‑1 (MQ=255)
     tCTGTGTTCCTGTAAATTCGGTGTTTGCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGc                                                                                                                                            <  1:1106370/149‑1 (MQ=255)
     tCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCACGTCTTATCAGACCTCCAACTATTACCGACCTGTAGGCATGATCAggg                                                                                                                                            >  1:808979/1‑148 (MQ=21)
       tgtgTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCAc                                                                                                                                          >  1:325582/1‑149 (MQ=255)
           ttCCTGTAAATTCGGTGTTGTCGGATGCACGACACGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTtgt                                                                                                                                      >  1:828915/1‑149 (MQ=255)
            tCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTtgtg                                                                                                                                     >  1:270306/1‑149 (MQ=255)
             ccagTAAATCCGGTCTTGTCGGATGCACGACCCGTAGGCCGGATACGGCGCCCGCGCCGCATCCGGCAATCACTGCCTGATGGGACGTTTTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGc                                                                                                                                    <  2:812998/146‑1 (MQ=17)
             ccTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGc                                                                                                                                    <  2:105095/149‑1 (MQ=255)
              cTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGcc                                                                                                                                   <  1:203131/149‑1 (MQ=255)
              cTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGcc                                                                                                                                   <  2:53365/149‑1 (MQ=255)
              cTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGcc                                                                                                                                   >  1:846824/1‑149 (MQ=255)
               tGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCg                                                                                                                                  >  1:415029/1‑149 (MQ=255)
                     ttCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCg                                                                                                                            <  1:965783/149‑1 (MQ=255)
                      tCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCgg                                                                                                                           >  1:1019784/1‑149 (MQ=255)
                         gtgtGGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGcaa                                                                                                                        <  2:658509/149‑1 (MQ=255)
                          tgttgtCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACCCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTATGCCGCAACCCGcgat                                                                                                                       >  2:929128/1‑149 (MQ=17)
                              gtCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCTCTTGTGCCGCATCCGGcaatgaat                                                                                                                   <  2:255365/149‑1 (MQ=255)
                              gtCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGGCGCATCCGGcaatcaat                                                                                                                   >  2:625522/1‑149 (MQ=255)
                              gtCGGATGCACGACCCGTAGGCCGGATAAGGCCAAAATACCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCCGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGcaatcaat                                                                                                                   >  1:562519/1‑149 (MQ=14)
                               tCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTAGAACTATTGCGGGCCTGTAGGGCGGGTGGGGGACTTGTGCCGCATTCGGcagtagatg                                                                                                                  >  2:472491/1‑142 (MQ=14)
                               tCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATg                                                                                                                  >  1:198503/1‑149 (MQ=255)
                                cGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGc                                                                                                                 <  1:713230/149‑1 (MQ=255)
                                 ggATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTAGCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGcc                                                                                                                <  1:722945/149‑1 (MQ=255)
                                  gATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAGTGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGGCCTGTTGGCCTGGTAAGGGTTTTTTGTCGCATGCGGCGATGAATTCCt                                                                                                               >  2:185717/1‑149 (MQ=14)
                                      cACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATg                                                                                                           <  1:212152/149‑1 (MQ=255)
                                       aCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGc                                                                                                          >  2:478368/1‑149 (MQ=255)
                                       aCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGc                                                                                                          >  1:432877/1‑149 (MQ=255)
                                        cGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCTCGCGGTGATGAATGACTGCTGCg                                                                                                         >  1:832052/1‑149 (MQ=14)
                                            ccGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGc                                                                                                     <  2:432877/149‑1 (MQ=25)
                                             cGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCt                                                                                                    <  1:651528/149‑1 (MQ=255)
                                                  gCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTcgc                                                                                               <  1:943131/149‑1 (MQ=40)
                                                   ccGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTCTCAGGCCTGCCACTATTGCCGACCTGTAGGCCTGATATGGCGCTTGTGCCGCATCCGGCAACCTATGCCTGATGCGACGCTGTcgcg                                                                                              >  2:707089/1‑149 (MQ=14)
                                                   ccGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTcgcg                                                                                              >  2:336563/1‑149 (MQ=255)
                                                    cGGATAAGGCGCTCGCGCCGCATCGGGCAATCACTTCCCGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATGGCCGACCCGTAGGCCGGATAAGGCCCTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGt                                                                                             <  2:247732/149‑1 (MQ=11)
                                                      gATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCt                                                                                           >  1:855862/1‑149 (MQ=255)
                                                        tAAGGCGCTCGCGACGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTa                                                                                         >  1:785763/1‑149 (MQ=255)
                                                          aGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATc                                                                                       <  1:388862/149‑1 (MQ=39)
                                                           ggCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTCTCa                                                                                      <  2:591202/149‑1 (MQ=255)
                                                            gcgcTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAg                                                                                     >  1:881551/1‑149 (MQ=39)
                                                            gcgcTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCATATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAg                                                                                     >  1:792800/1‑149 (MQ=255)
                                                             cgcTCGCGCCGCAGCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAgg                                                                                    <  1:791210/149‑1 (MQ=38)
                                                               cTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATa                                                                                                                                                >  2:1150166/1‑87 (MQ=40)
                                                               cTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATa                                                                                                                                                >  2:1150291/1‑87 (MQ=14)
                                                                tCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTTTTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCt                                                                                 <  1:840610/149‑1 (MQ=25)
                                                                tCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTTTCAGGCCt                                                                                 >  1:127475/1‑149 (MQ=25)
                                                                  gcgcCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTAc                                                                               <  2:828915/149‑1 (MQ=33)
                                                                         ctCCGGCAATCAGTGCCTGTTGCGGCGCTGTCGCGGCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGtt                                                                        <  2:934372/148‑1 (MQ=255)
                                                                          tCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTTTTg                                                                       <  2:198503/149‑1 (MQ=17)
                                                                           ccGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGTATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCtt                                                                                          >  1:599545/1‑129 (MQ=17)
                                                                           ccGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTCCAACTGTTGc                                                                      >  2:593582/1‑149 (MQ=17)
                                                                           ccGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACACCTGTTGc                                                                      >  1:761822/1‑149 (MQ=17)
                                                                           ccGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGc                                                                      >  1:1142232/1‑149 (MQ=32)
                                                                            cGGCAATCAATGCCGGGTGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGGGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGcc                                                                     <  2:181800/149‑1 (MQ=12)
                                                                             ggCAATCAATGCCTGATGCGACGGTGTCGCATCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTTCCg                                                                    >  1:604730/1‑149 (MQ=12)
                                                                                 atcgatGCCTGGTGCGACGCTGTCGCGTCGTATCAGTCCTACAACTATGGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGAGGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCt                                                                <  2:706856/149‑1 (MQ=11)
                                                                                   caatGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGAcccgt                                                              >  2:592350/1‑146 (MQ=1)
                                                                                   caatGCCCGATGCGACGCTGTCGCGTCTTACGAGGCCTACAACTATGGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCAGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGt                                                              <  1:773439/149‑1 (MQ=11)
                                                                                    aatGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTcgc                                                                                               >  1:568985/1‑115 (MQ=32)
                                                                                     atGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAg                                                            >  1:685908/1‑149 (MQ=35)
                                                                                       gCCGGATGCGGCGCTGTCGCGTCTTAGCAGGCCTACAACGCGTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGc                                                          <  2:417493/149‑1 (MQ=11)
                                                                                        ccTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGcc                                                         >  2:1030246/1‑149 (MQ=14)
                                                                                           gATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGa                                                      <  1:197537/149‑1 (MQ=35)
                                                                                              gCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATaa                                                   <  1:458956/149‑1 (MQ=34)
                                                                                              gCGACGCTGTCGCGTCTTATCAGGCCTACAACTAGTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATaa                                                   <  1:988072/149‑1 (MQ=25)
                                                                                               cGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTTATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAg                                                  <  2:400621/149‑1 (MQ=25)
                                                                                                tccggtGTCCCGTTGTATCAGGCCTACAACTATTGCTGACCGGTAGGCCTGATAGCGCACTTGTGCCCCATCCGCCAATCAATTCCTGATCCGACGCTGACGCGTCTTATCAGGCCTACAACTGTTGCCGACCAGTAGCCCGGATaacg                                                 <  2:322031/144‑3 (MQ=11)
                                                                                                gACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGAAAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTACGCCGGATAAgg                                                 >  1:898851/1‑149 (MQ=21)
                                                                                                   gCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCAc                                              >  2:384729/1‑149 (MQ=37)
                                                                                                     tGTCGCGTCTTATCAGGCCTACAACTATTGCCTACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACtt                                            <  2:808979/149‑1 (MQ=37)
                                                                                                     tGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACtt                                            <  2:299871/149‑1 (MQ=38)
                                                                                                         gcgTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTCTCAGGCCTACAACGGTTGCCGACCTGTCGGCCGGGTAAGGCACTTGTGc                                        >  1:402015/1‑149 (MQ=255)
                                                                                                          cgTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGcc                                       >  1:939816/1‑149 (MQ=38)
                                                                                                                aTCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGcc                                                                     >  2:189064/1‑113 (MQ=33)
                                                                                                                aTCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATcc                                 >  1:606270/1‑149 (MQ=38)
                                                                                                                 tCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCg                                <  1:748948/149‑1 (MQ=37)
                                                                                                                     gCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCGCTTGTGCCGCATCCGGcaa                            >  1:1147623/1‑149 (MQ=38)
                                                                                                                     gCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGcaa                            >  1:1147582/1‑149 (MQ=38)
                                                                                                                          cAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCACTTGCGCCGCATCCGGcagtcaat                       >  2:768911/1‑149 (MQ=16)
                                                                                                                            ccgaTTGCCGACCTGTAGGCCGGATAAGGCACTTGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGc                     <  2:415029/146‑1 (MQ=9)
                                                                                                                              tATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCt                   >  2:662394/1‑149 (MQ=37)
                                                                                                                                   ccGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGc              >  1:687292/1‑149 (MQ=35)
                                                                                                                                     gACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGCCCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCCATGCCCGATGCGa            >  2:420703/1‑149 (MQ=255)
                                                                                                                                      aCCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGAc           >  2:877107/1‑149 (MQ=35)
                                                                                                                                        cTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGc         >  1:138801/1‑149 (MQ=35)
                                                                                                                                             ggCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTcgc                                                                                               >  2:1067059/1‑58 (MQ=38)
                                                                                                                                              gCCTGATAAGGGACTTGTGCCCCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTcgc   <  2:243018/149‑1 (MQ=21)
                                                                                                                                              gCCTGATAAGGCACTTGTGCCGCATGCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGTCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTcgc   <  1:707089/149‑1 (MQ=25)
                                                                                                                                              gCCTGATAAGGCACTTGTGCCGCATCCGGCAATCGATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGAGAAGGCACTGATGCCGCATGCGGCAATCAATGCCTGATGCGACGCTGTcgc   >  1:439645/1‑149 (MQ=11)
                                                                                                                                              gCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTCGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTcgc   >  1:515039/1‑149 (MQ=25)
                                                                                                                                               ccTTATATGGCACTTGCGCCGCATCCGGCAATCAATGCCTGATGCGACGTTGTCGCGTCTTATCAGGCCAACAACGGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTcgcg  <  1:329871/149‑1 (MQ=11)
                                                                                                                                               ccTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTc                                                                                            >  1:65121/1‑59 (MQ=35)
                                                                                                                                                cTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCt                                                                >  1:222714/1‑86 (MQ=33)
                                                                                                                                                 tGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTa                                                             >  1:712436/1‑88 (MQ=35)
                                                                                                                                                 |                                                                                                                                                  
CATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCG  >  NC_000913/4094269‑4094560

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: