Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,439,556 G→A 6.0% intergenic (‑294/‑31) ytfJ ← / → ytfK putative transcriptional regulator/DUF1107 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,439,5560GA 6.0% 264.8 / 4.1 82intergenic (‑294/‑31)ytfJ/ytfKputative transcriptional regulator/DUF1107 family protein
Reads supporting (aligned to +/- strand):  ref base G (35/42);  new base A (3/2);  total (38/44)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GACCGCGCTTCGGTAAAAAGATAATTCTGAATAATTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACATCAAGGACGTTGGCGCTTTTGAATTCGATAAGGGTAAGATT  >  NC_000913/4439410‑4439703
                                                                                                                                                  |                                                                                                                                                   
gACCGCGCTTCGGTAAAAAGATAATTCTGAATAATTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGc                                                                                                                                                    <  5:309623/148‑1 (MQ=255)
 aCCGCGCTTCGGTAAAAAGATAATTCTGAATAATTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGca                                                                                                                                                   <  2:350523/148‑1 (MQ=255)
  cctcgcTTCGGTAAAAAGATAATTCTGAATAATTGTAACCTTTATGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGcac                                                                                                                                                  <  4:99275/145‑1 (MQ=255)
  ccGCGCTTCGGTAAAAAGATAATTCTGAATAATTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATATTAGTTGGGCCGATTGTGGCACCGcac                                                                                                                                                  <  6:63251/148‑1 (MQ=255)
  ccGCGCTTCGGTAAAAAGATAATTCTGAATAATTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGcac                                                                                                                                                  <  2:82899/148‑1 (MQ=255)
  ccGCGCTTCGGTAAAAAGATAATTCTGAATAATTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGcac                                                                                                                                                  <  7:101295/148‑1 (MQ=255)
  ccGCGCTTCGGTAAAAAGATAATTCTGAATAATTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGcac                                                                                                                                                  <  6:177931/148‑1 (MQ=255)
  ccGCGCTTCGGTAAAAAGATAATTCTGAATAATTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGca                                                                                                                                                   <  8:96672/147‑1 (MQ=255)
  ccGCGCTTCGGTAAAAAGATAATTCTGAATAATTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGca                                                                                                                                                   >  5:334424/1‑147 (MQ=255)
  ccGCGCTTCGGTAAAAAGATAATTCTGAATAATTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGca                                                                                                                                                   >  6:104261/1‑147 (MQ=255)
  ccGCGCTTCGGTAAAAAGATAATTCTGAATAATTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGca                                                                                                                                                   >  1:357851/1‑147 (MQ=255)
  ccGCGCTTCGGTAAAAAGATAATTCTGAATAATTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGCGGCACCGcac                                                                                                                                                  <  1:306863/148‑1 (MQ=255)
          cGGTAAAAAGATAATTCTGAATAATTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTa                                                                                                                                           <  7:250495/147‑1 (MQ=255)
           ggTAAAAAGATAATTCTGAATAATTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGcac                                                                                                                                                  >  6:166690/1‑139 (MQ=255)
           ggTAAAAAGATAATTCTGAATAATTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGcac                                                                                                                                                  <  5:166690/139‑1 (MQ=255)
           ggTAAAAAGATAATTCTGAATAATTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGGCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAAt                                                                                                                                         <  6:337877/148‑1 (MQ=255)
             tAAAAAGATAATTCTGAATAATTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGcac                                                                                                                                                  <  6:358564/137‑1 (MQ=255)
             tAAAAAGATAATTCTGAATAATTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGcac                                                                                                                                                  >  5:358564/1‑137 (MQ=255)
                         tCTGAATAATTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGgaga                                                                                                                            >  1:332702/1‑147 (MQ=255)
                                aatTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCCCAGGCGTAATACTCAGCAGGAGATAACAAt                                                                                                                     >  5:380186/1‑147 (MQ=255)
                                     tAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAAt                                                                                                               <  1:137158/148‑1 (MQ=255)
                                     tAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAAt                                                                                                               <  3:320401/148‑1 (MQ=255)
                                           ttAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCACACAGGCGTAATACTcagc                                                                                                                                  >  6:301280/1‑123 (MQ=255)
                                           ttAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGCTGGGCCGATTGTGGCACCACACAGGCGTAATACTcagc                                                                                                                                  <  5:301280/123‑1 (MQ=255)
                                            tAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCa                                                                                                         <  7:206947/147‑1 (MQ=255)
                                              gttAAAAAAAGTTATACTCGGGGGAAACATTGCCCGTATAGTCTATAGTCACTAAGCATTAAAATTTTCGCCTCATAATAGTTGGGCCGATTGTGGCAACGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAAc                                                                                                       <  5:45117/145‑1 (MQ=255)
                                                   aaaaaGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACa                                                                                                 <  5:104261/148‑1 (MQ=255)
                                                   aaaaaGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGACTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACa                                                                                                 <  2:332702/148‑1 (MQ=255)
                                                                      aaaCATTGCCCGGATAGTCTATAGTCACTAAGCAATAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCg                                                                                >  8:198258/1‑146 (MQ=255)
                                                                       aaCATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAAt                                                                                                               <  4:131184/114‑1 (MQ=255)
                                                                       aaCATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAAt                                                                                                               >  3:131184/1‑114 (MQ=255)
                                                                                 ggATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCTCAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAGCGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGa                                                                   <  6:380186/148‑1 (MQ=255)
                                                                                   aTAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATc                                                                 >  8:154330/1‑148 (MQ=255)
                                                                                     aGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCt                                                               >  1:187024/1‑148 (MQ=255)
                                                                                        cTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCACACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGtt                                                            >  8:123853/1‑148 (MQ=255)
                                                                                          ataGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTcc                                                          >  5:223235/1‑148 (MQ=255)
                                                                                             gTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGAt                                                                  <  2:275788/137‑1 (MQ=255)
                                                                                             gTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGAt                                                                  >  1:275788/1‑137 (MQ=255)
                                                                                             gTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGt                                                        <  2:187024/147‑1 (MQ=255)
                                                                                               cACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACtt                                                                                         >  2:351438/1‑112 (MQ=255)
                                                                                               cACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACtt                                                                                         <  1:351438/112‑1 (MQ=255)
                                                                                                   aaGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACaa                                                                                                >  8:265604/1‑101 (MQ=255)
                                                                                                   aaGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACaa                                                                                                <  7:265604/101‑1 (MQ=255)
                                                                                                   aaGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGt                                                 >  4:121233/1‑148 (MQ=255)
                                                                                                   aaGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGt                                                 >  4:4337/1‑148 (MQ=255)
                                                                                                      cATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTtat                                              <  3:4337/148‑1 (MQ=255)
                                                                                                        ttAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATc                                                                 >  2:252702/1‑127 (MQ=255)
                                                                                                        ttAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATc                                                                 <  1:252702/127‑1 (MQ=255)
                                                                                                         tAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACa                                           >  8:310790/1‑148 (MQ=255)
                                                                                                         tAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCATAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACa                                           >  3:331748/1‑148 (MQ=255)
                                                                                                         tAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCACACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACa                                           >  1:207367/1‑148 (MQ=255)
                                                                                                          aaaaTTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAAGCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACAt                                          >  4:159906/1‑148 (MQ=255)
                                                                                                          aaaaTTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCg                                                         >  7:5804/1‑133 (MQ=255)
                                                                                                          aaaaTTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCg                                                         <  8:5804/133‑1 (MQ=255)
                                                                                                          aaaaTTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTAAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACAt                                          >  6:283902/1‑148 (MQ=255)
                                                                                                           aaaTTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGt                                                        <  8:33665/133‑1 (MQ=255)
                                                                                                           aaaTTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGt                                                        >  7:33665/1‑133 (MQ=255)
                                                                                                             aTTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACATCaa                                       <  2:357851/148‑1 (MQ=255)
                                                                                                             aTTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACATCaa                                       <  6:334424/148‑1 (MQ=255)
                                                                                                                      tCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACATCAAGGACGTTgg                              >  6:322487/1‑148 (MQ=255)
                                                                                                                       cATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCa                                                                                       <  5:129837/90‑1 (MQ=255)
                                                                                                                       cATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCa                                                                                       >  6:129837/1‑90 (MQ=255)
                                                                                                                       cATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTtata                                             >  7:201153/1‑132 (MQ=255)
                                                                                                                       cATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTtata                                             <  8:201153/132‑1 (MQ=255)
                                                                                                                       cATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACATCAAGGACGTTg                               >  8:5731/1‑146 (MQ=255)
                                                                                                                       cATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGGGTACGGTTATACATCAAGGACGTTGgc                             >  4:273729/1‑148 (MQ=255)
                                                                                                                       cATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATACTGTTCCGTGGACGGTTATACATCAAGGACGTTGgc                             >  2:304492/1‑148 (MQ=255)
                                                                                                                        ataataGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACATCAAGGACGTTgg                              <  7:5731/146‑1 (MQ=255)
                                                                                                                         taataGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAAttt                                                                                                             <  1:260759/66‑1 (MQ=255)
                                                                                                                         taataGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAAttt                                                                                                             >  2:260759/1‑66 (MQ=255)
                                                                                                                           ataGTTGGGCCGATTGTGGCACCACACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACATCAAGGACGTTGGCGCtt                         <  7:123853/148‑1 (MQ=255)
                                                                                                                               ttGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACATCAAGGACGTTGGCGCTTTTGa                     <  7:198258/148‑1 (MQ=255)
                                                                                                                                 gggCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAAGCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACATCAAGGACGTTGGCGCTTTTGaa                    <  3:159906/147‑1 (MQ=255)
                                                                                                                                  ggCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACATCAAGGACGTTGGCGCTTTTGAAtt                  <  7:154330/148‑1 (MQ=255)
                                                                                                                                    ccGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACATCAAGGACGTTGGCGCTTTTGAATTc                 >  8:61890/1‑147 (MQ=255)
                                                                                                                                    ccGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACATCAAGGACGTTGGCGCTTTTGAATTCg                <  7:310790/148‑1 (MQ=255)
                                                                                                                                      gATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATCACAATGAAAATTTTCCAACGGTACACCCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACATCAAGGACGTTGGCGCTTTTGAATTCGAt              <  5:239796/148‑1 (MQ=255)
                                                                                                                                       aTTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTc                                                                                        <  8:270962/73‑1 (MQ=255)
                                                                                                                                       aTTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTc                                                                                        >  7:270962/1‑73 (MQ=255)
                                                                                                                                        ttGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACATCAAGGACGTTGGCGCTTTTGAATTCGAt              >  1:286170/1‑146 (MQ=255)
                                                                                                                                          gtgGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACATCAAGGACGTTGgc                             >  7:318400/1‑129 (MQ=255)
                                                                                                                                          gtgGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACATCAAGGACGTTGgc                             <  8:318400/129‑1 (MQ=255)
                                                                                                                                                  gCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACATCAAGGACGTTGGCGCTTTTGAATTCGATAAGGGTAAGAtt  >  5:278353/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GACCGCGCTTCGGTAAAAAGATAATTCTGAATAATTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACATCAAGGACGTTGGCGCTTTTGAATTCGATAAGGGTAAGATT  >  NC_000913/4439410‑4439703

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