Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,127,938 T→A 7.2% intergenic (‑44/+140) yiiX ← / ← metJ putative lipid binding hydrolase, DUF830 family protein/transcriptional repressor, S‑adenosylmethionine‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,127,9380TA 7.2% 290.7 / 11.3 83intergenic (‑44/+140)yiiX/metJputative lipid binding hydrolase, DUF830 family protein/transcriptional repressor, S‑adenosylmethionine‑binding
Reads supporting (aligned to +/- strand):  ref base T (31/46);  new base A (3/3);  total (34/49)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.65e-01

GATGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTA  >  NC_000913/4127794‑4128083
                                                                                                                                                |                                                                                                                                                 
gATGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTc                                                                                                                                               >  4:1252/1‑149 (MQ=255)
gATGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTc                                                                                                                                               <  3:452847/149‑1 (MQ=255)
 aTGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTc                                                                                                                                               >  2:382709/1‑148 (MQ=255)
 aTGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTc                                                                                                                                               >  2:266262/1‑148 (MQ=255)
 aTGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTc                                                                                                                                               >  1:756230/1‑148 (MQ=255)
  tGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCg                                                                                                                                              <  2:533272/148‑1 (MQ=255)
   gAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCgtt                                                                                                                                            <  3:147651/149‑1 (MQ=255)
   gAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCgt                                                                                                                                             >  1:634872/1‑148 (MQ=255)
     gCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCgttgt                                                                                                                                          <  3:44612/149‑1 (MQ=255)
     gCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCgttgt                                                                                                                                          <  4:372264/149‑1 (MQ=255)
      cGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACTAGCACgttgtc                                                                                                                                         <  3:725208/149‑2 (MQ=255)
      cGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCgttgtt                                                                                                                                         <  4:142832/149‑1 (MQ=255)
      cGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGAGCACgttgt                                                                                                                                          <  1:847989/148‑1 (MQ=255)
       gagagaTCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCgttgtt                                                                                                                                         <  1:382709/148‑1 (MQ=255)
       gagagaTCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGAGCACgttgtct                                                                                                                                        <  4:702190/149‑3 (MQ=255)
           gaTCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATg                                                                                                                                     >  1:842538/1‑148 (MQ=255)
           gaTCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATg                                                                                                                                     <  1:805588/148‑1 (MQ=255)
               tGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGa                                                                                                                                >  4:289452/1‑149 (MQ=255)
                gAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGa                                                                                                                                >  2:465936/1‑148 (MQ=255)
                gAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGa                                                                                                                                >  2:312850/1‑148 (MQ=255)
                  aaGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATg                                                                                                                              >  1:648736/1‑148 (MQ=255)
                   aGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCg                                                                                                                                  >  3:153679/1‑143 (MQ=255)
                   aGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCg                                                                                                                                  <  4:153679/143‑1 (MQ=255)
                     atgatgTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATgcgg                                                                                                                           <  2:434875/148‑1 (MQ=255)
                     atgatgTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGAGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATgcgg                                                                                                                           <  2:457428/148‑1 (MQ=255)
                       gatgTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATgcggcg                                                                                                                         <  2:36085/148‑1 (MQ=255)
                        atgTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGt                                                                                                                        <  2:423531/148‑1 (MQ=255)
                           tCGCCGGTTTGTGGCTTCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGaa                                                                                                                     <  2:733042/148‑1 (MQ=255)
                           tCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAAc                                                                                                                    <  4:745737/149‑1 (MQ=255)
                             gCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGc                                                                                                                  <  3:789741/149‑1 (MQ=255)
                              ccGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGc                                                                                                                  >  2:698129/1‑148 (MQ=255)
                              ccGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGc                                                                                                                  <  2:261404/148‑1 (MQ=255)
                                 gTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGACCGCCtt                                                                                                               <  2:115566/148‑1 (MQ=255)
                                             cAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACaa                                                                                                  <  4:270649/149‑1 (MQ=255)
                                             cAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACa                                                                                                   <  1:693461/148‑1 (MQ=255)
                                                   aaGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCg                                                                                             <  2:756230/148‑1 (MQ=255)
                                                      gCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCa                                                                                         >  4:819570/1‑149 (MQ=255)
                                                        aGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCaaa                                                                                       <  4:854921/149‑1 (MQ=255)
                                                          gTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAATCGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAAt                                                                                      <  2:642391/148‑1 (MQ=255)
                                                          gTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAAt                                                                                      <  1:465936/148‑1 (MQ=255)
                                                             cAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGcc                                                                                                                 >  3:652166/1‑118 (MQ=255)
                                                             cAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGcc                                                                                                                 <  4:652166/118‑1 (MQ=255)
                                                                  aCCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCTATaaa                                                                              <  2:648736/148‑1 (MQ=255)
                                                                       cagGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCaa                                                                        <  3:819570/149‑1 (MQ=255)
                                                                           cTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCaa                                                                                        >  2:846971/1‑129 (MQ=255)
                                                                           cTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCaa                                                                                        <  1:846971/129‑1 (MQ=255)
                                                                           cTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAatat                                                                     >  2:560998/1‑148 (MQ=255)
                                                                             gAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGAc                                                                  <  3:732712/149‑1 (MQ=255)
                                                                              aGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGAc                                                                  >  2:465252/1‑148 (MQ=255)
                                                                                gATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGt                                                                <  2:634872/148‑1 (MQ=255)
                                                                                  tCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAgg                                                             <  4:646608/149‑1 (MQ=255)
                                                                                   cagcagCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAgg                                                             >  1:765171/1‑148 (MQ=255)
                                                                                   cagcagCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAgg                                                             >  2:144632/1‑148 (MQ=255)
                                                                                        gcCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGa                                                        <  2:868232/148‑1 (MQ=255)
                                                                                        gcCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGa                                                        <  1:698129/148‑1 (MQ=255)
                                                                                        gcCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGAt                                                       >  4:10748/1‑149 (MQ=255)
                                                                                         ccTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGAt                                                       >  2:487876/1‑148 (MQ=255)
                                                                                          cTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATa                                                      <  1:487876/148‑1 (MQ=255)
                                                                                             tttttCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCg                                                  >  3:313937/1‑149 (MQ=255)
                                                                                                       ttAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATc                                         <  1:343427/148‑1 (MQ=255)
                                                                                                            cGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCg                                    <  1:266262/148‑1 (MQ=255)
                                                                                                             gTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGa                                   <  1:312850/148‑1 (MQ=255)
                                                                                                               ccATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGAtt                                 <  1:560998/148‑1 (MQ=255)
                                                                                                                 aTGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCa                              >  3:310029/1‑149 (MQ=255)
                                                                                                                    tataAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTa                                                >  3:459381/1‑128 (MQ=255)
                                                                                                                    tataAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTa                                                <  4:459381/128‑1 (MQ=255)
                                                                                                                     ataAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTc                           >  1:306963/1‑148 (MQ=255)
                                                                                                                         aaaGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCcgcg                       >  1:197001/1‑148 (MQ=255)
                                                                                                                          aaGCGGTGGGTCGCAGACAACGTGCTCGTGGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGcgcc                     >  4:698765/1‑149 (MQ=255)
                                                                                                                                 gggTCGCAGCCAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCttt              >  4:704634/1‑149 (MQ=255)
                                                                                                                                 gggTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCtt               <  1:465252/148‑1 (MQ=255)
                                                                                                                                   gTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCttttt            <  3:289452/149‑1 (MQ=255)
                                                                                                                                   gTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCtttt             <  2:765171/148‑1 (MQ=255)
                                                                                                                                       taaaCAACGAGCACGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGc         >  1:735848/4‑148 (MQ=255)
                                                                                                                                       cAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGc         >  2:393159/1‑148 (MQ=255)
                                                                                                                                        aGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCt        <  2:842538/148‑1 (MQ=255)
                                                                                                                                         gACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCttt      <  4:137796/149‑1 (MQ=255)
                                                                                                                                         gACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCttt      >  3:546902/1‑149 (MQ=255)
                                                                                                                                          aCAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCttt      <  1:393159/148‑1 (MQ=255)
                                                                                                                                          aCAACGAGCACGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCttt      >  2:394873/1‑148 (MQ=255)
                                                                                                                                          aCAACGAGCACGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCttt      >  1:389271/1‑148 (MQ=255)
                                                                                                                                           cAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTa     >  1:456353/1‑148 (MQ=255)
                                                                                                                                             aCGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTa  <  3:1252/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
GATGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTA  >  NC_000913/4127794‑4128083

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