Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,094,327 A→G 6.5% intergenic (‑55/+396) 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,3270AG 6.5% 256.2 / 5.8 77intergenic (‑55/+396)rhaD/rhaArhamnulose‑1‑phosphate aldolase/L‑rhamnose isomerase
Reads supporting (aligned to +/- strand):  ref base A (31/41);  new base G (3/2);  total (34/43)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

GCCGCCGTTGCGCTCATCCCAGCCTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTA  >  NC_000913/4094180‑4094474
                                                                                                                                                   |                                                                                                                                                   
gccgccGTTGCGCTCCTCCCAGCCTTTCAGCCAGGCGTCGGTGGTGGCTTTGTTCATTCCCTGGACAAACCAGGACTGAATAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATaa                                                                                                                                                    <  4:209189/149‑1 (MQ=255)
gccgccGTTGCGCTCATCCCAGCCTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATa                                                                                                                                                     <  1:502580/148‑1 (MQ=255)
 ccgccgTTGCGCTCATCCCAGCCTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAg                                                                                                                                                   <  4:862501/149‑1 (MQ=255)
   gccgTTGCGCTCATCCCAGCCTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAgg                                                                                                                                                  <  2:552217/148‑1 (MQ=255)
   gccgTTGCGCTCATCCCAGCCTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAgg                                                                                                                                                  >  1:731971/1‑148 (MQ=255)
     cgTTGCGCTCATCCCAGCCTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGgcgc                                                                                                                                               <  4:279805/149‑1 (MQ=255)
           gcTCATCCCAGCCTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCgcgc                                                                                                                                         <  4:361817/149‑1 (MQ=255)
               aTCCCAGCCTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGcgccgc                                                                                                                                      >  2:77747/1‑148 (MQ=255)
                     gCCTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCgg                                                                                                                                >  2:186365/1‑148 (MQ=255)
                       cTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAATGCGCTCGCGCCGCATCCGGca                                                                                                                              <  1:103590/148‑1 (MQ=255)
                            aGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGcaatcaa                                                                                                                         <  1:186365/148‑1 (MQ=255)
                              ccAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATg                                                                                                                       <  1:54766/148‑1 (MQ=255)
                              ccAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGc                                                                                                                      <  3:200605/149‑1 (MQ=255)
                              ccAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGc                                                                                                                      <  3:179684/149‑1 (MQ=255)
                                  gCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTGCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGa                                                                                                                  >  4:789267/1‑149 (MQ=255)
                                    gTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATg                                                                                                                <  4:711037/149‑1 (MQ=255)
                                         tggtggCTTTGATCATTCCCTGGACAAACCAGGCCTGAGTAATGCTTTGCATGTACTGTGTTCCTGTAAATTCGGAGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGAc                                                                                                            >  2:849900/1‑148 (MQ=255)
                                              gCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTc                                                                                                      <  3:407778/149‑1 (MQ=255)
                                                     tCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCt                                                                                                >  1:797814/1‑148 (MQ=255)
                                                        ttCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGcaatc                                                                                                                           <  3:849695/118‑1 (MQ=255)
                                                        ttCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGcaatc                                                                                                                           >  4:849695/1‑118 (MQ=255)
                                                         tCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATc                                                                                            <  2:797814/148‑1 (MQ=255)
                                                         tCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATc                                                                                            <  1:239810/148‑1 (MQ=255)
                                                            cTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAgg                                                                                         >  2:813154/1‑148 (MQ=255)
                                                             tGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCt                                                                                                <  1:28133/140‑1 (MQ=255)
                                                             tGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCt                                                                                                >  2:28133/1‑140 (MQ=255)
                                                               gACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGt                                                                                                       >  2:886605/1‑131 (MQ=255)
                                                               gACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGt                                                                                                       <  1:886605/131‑1 (MQ=255)
                                                                aCAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTa                                                                                     >  1:679137/1‑148 (MQ=255)
                                                                  aaaCCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAc                                                                                 >  3:781923/1‑149 (MQ=255)
                                                                     ccAGGACTGAGTAATGCTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTa                                                                               <  3:771905/149‑1 (MQ=255)
                                                                      cAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTa                                                                               <  2:679137/148‑1 (MQ=255)
                                                                      cAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTa                                                                               <  1:51146/148‑1 (MQ=255)
                                                                       aGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTAt                                                                              <  1:77747/148‑1 (MQ=255)
                                                                       aGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTAt                                                                              <  2:731971/148‑1 (MQ=255)
                                                                         gACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTg                                                                            <  1:813154/148‑1 (MQ=255)
                                                                                 aaTGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATGCGGCGCGAGcgcctta                                                                                                                                     >  1:782279/1‑80 (MQ=38)
                                                                                 aaTGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATGCGGCGCGAGcgcctta                                                                                                                                     <  2:782279/83‑4 (MQ=38)
                                                                                 aaTGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGt                                                                   >  3:130428/1‑149 (MQ=255)
                                                                                  aTGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGt                                                                   >  1:47409/1‑148 (MQ=255)
                                                                                   tGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATg                                                                                                                       >  2:282651/1‑95 (MQ=255)
                                                                                   tGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATg                                                                                                                       <  1:282651/95‑1 (MQ=255)
                                                                                       ttGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGcc                                                              <  2:650516/148‑1 (MQ=255)
                                                                                       ttGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGcc                                                              <  2:806694/148‑1 (MQ=255)
                                                                                       ttGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCt                                                             <  3:69611/149‑1 (MQ=255)
                                                                                          cATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATc                                                                                            >  1:491815/1‑115 (MQ=255)
                                                                                          cATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATc                                                                                            <  2:491815/115‑1 (MQ=255)
                                                                                          cATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTg                                                            <  1:135556/147‑1 (MQ=255)
                                                                                          cATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTg                                                            >  2:135556/1‑147 (MQ=255)
                                                                                            tGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATGCGGCGCGAGCGCCTTATCCgg                                                                                                                                >  1:172367/1‑77 (MQ=25)
                                                                                            tGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATGCGGCGCGAGCGCCTTATCCgg                                                                                                                                <  2:172367/77‑1 (MQ=21)
                                                                                            tGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATaa                                                        <  3:294799/149‑1 (MQ=255)
                                                                                            tGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATa                                                         <  2:47409/148‑1 (MQ=255)
                                                                                              ttCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAgg                                                      >  3:815247/1‑149 (MQ=255)
                                                                                              ttCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAg                                                       <  2:420140/148‑1 (MQ=255)
                                                                                               tCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGc                                                     <  4:815247/149‑1 (MQ=255)
                                                                                               tCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCa                                                    <  4:781923/149‑1 (MQ=255)
                                                                                                 tgtgTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCAc                                                   <  4:337906/149‑1 (MQ=255)
                                                                                                        cTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGcc                                            >  3:419430/1‑149 (MQ=255)
                                                                                                         tGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGc                                             >  3:230614/1‑147 (MQ=255)
                                                                                                         tGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGc                                             <  4:230614/147‑1 (MQ=255)
                                                                                                         tGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCg                                           >  4:63684/1‑149 (MQ=255)
                                                                                                          gTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTTGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGcc                                                                                       <  1:494184/104‑1 (MQ=255)
                                                                                                          gTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTTGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGcc                                                                                       >  2:494184/1‑104 (MQ=255)
                                                                                                               ttCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATcc                                      >  2:58731/1‑148 (MQ=255)
                                                                                                                         tCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATg                           >  3:392136/1‑149 (MQ=255)
                                                                                                                            gATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGcc                         >  1:424155/1‑148 (MQ=255)
                                                                                                                                 aCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGc                          >  4:197513/1‑142 (MQ=255)
                                                                                                                                 aCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGc                          <  3:197513/142‑1 (MQ=255)
                                                                                                                                 aCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGc                   >  4:291703/1‑149 (MQ=255)
                                                                                                                                  cGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGc                   >  1:693121/1‑148 (MQ=255)
                                                                                                                                  cGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCg                  >  3:151073/1‑149 (MQ=255)
                                                                                                                                         attgCCGGATGCGGCGCGAGCGCCTTATCCGGCCATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGt           >  3:298431/4‑149 (MQ=11)
                                                                                                                                            gCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTcg         <  1:58731/148‑1 (MQ=255)
                                                                                                                                             ccGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGt           <  1:88955/145‑1 (MQ=39)
                                                                                                                                             ccGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGt           >  2:88955/1‑145 (MQ=40)
                                                                                                                                                 aTAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGCCGCTGTCGCGTCtt   <  3:63684/149‑1 (MQ=255)
                                                                                                                                                 aTAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCt    <  1:849900/148‑1 (MQ=39)
                                                                                                                                                  tAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGTCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTa  >  3:865889/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCCGCCGTTGCGCTCATCCCAGCCTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTA  >  NC_000913/4094180‑4094474

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