Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,094,335 C→A 6.7% intergenic (‑63/+388) 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,3350CA 6.7% 250.6 / 8.0 76intergenic (‑63/+388)rhaD/rhaArhamnulose‑1‑phosphate aldolase/L‑rhamnose isomerase
Reads supporting (aligned to +/- strand):  ref base C (30/39);  new base A (3/2);  total (34/42)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

GCTCATCCCAGCCTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCC  >  NC_000913/4094191‑4094481
                                                                                                                                                |                                                                                                                                                  
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                                                                                                                                  >  4:849695/1‑118 (MQ=255)
                                             ttCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGcaatc                                                                                                                                  <  3:849695/118‑1 (MQ=255)
                                              tCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATc                                                                                                   <  1:239810/148‑1 (MQ=255)
                                              tCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATc                                                                                                   <  2:797814/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                                                                                                                                            <  2:782279/83‑4 (MQ=38)
                                                                      aaTGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATGCGGCGCGAGcgcctta                                                                                                                                            >  1:782279/1‑80 (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:806694/148‑1 (MQ=255)
                                                                            ttGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGcc                                                                     <  2:650516/148‑1 (MQ=255)
                                                                            ttGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCt                                                                    <  3:69611/149‑1 (MQ=255)
                                                                               cATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATc                                                                                                   <  2:491815/115‑1 (MQ=255)
                                                                               cATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATc                                                                                                   >  1:491815/1‑115 (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                                                    <  4:230614/147‑1 (MQ=255)
                                                                                              tGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGc                                                    >  3:230614/1‑147 (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                                 <  3:197513/142‑1 (MQ=255)
                                                                                                                      aCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGc                                 >  4:197513/1‑142 (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)
                                                                                                                                            cgcTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAg     >  2:358790/1‑148 (MQ=38)
                                                                                                                                             gcTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGc   <  4:392136/149‑1 (MQ=37)
                                                                                                                                               tCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGcc  <  2:424155/148‑1 (MQ=35)
                                                                                                                                                |                                                                                                                                                  
GCTCATCCCAGCCTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCC  >  NC_000913/4094191‑4094481

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: