Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,838,231 A→T 11.1% intergenic (‑126/+23) hydN ← / ← ascG formate dehydrogenase‑H, [4Fe‑4S] ferredoxin subunit/asc operon transcriptional repressor; prpBC operon repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,838,2310AT11.1% 222.9 / 18.6 72intergenic (‑126/+23)hydN/ascGformate dehydrogenase‑H, [4Fe‑4S] ferredoxin subunit/asc operon transcriptional repressor; prpBC operon repressor
Reads supporting (aligned to +/- strand):  ref base A (37/27);  new base T (5/3);  total (42/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.49e-01

ATGGATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGAACGGTGTAGGGGGCAATAGCGATATCATCGAATCCGATAACTGACACCTGCTCTGGCACCGCTACGCCGCGCTCGTGTAACGCTTTCATCG  >  NC_000913/2837971‑2838501
                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                              
aTGGATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTGCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTTCAGCAATGACAGTTAATGCGCCGATACATTCACAAATTATGCCAGTTTTTAATTTACTTTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGTGGTCATCGACACGCGTGATGTAGACATGtgat                                                                                                                                                                                                                                                                  >  4:413481/1‑275 (MQ=255)
aTGGATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTACTGATTTTCCTGATGAGACACCACGCACGCTACTTCACAGGTCCGGCAACCAATACATTTACTCGCGTAAGCAATGATGAAACGGTTACTCAAATTCTCCCGCAATGAAAGTTAATGCGCCGCTACATTCACAACTCATGCCAGTTTTTAATTTACTTGTATTTAAGGAAATTAATTACTGTAATGCAGGAAAAAAGATGTTATCGACCGTAGTTACTAAGACAagttaa                                                                                                                                                                                                                                                                  >  4:532991/1‑272 (MQ=255)
                   gTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAg                                                                                                                                                                                                                                               <  2:413482/275‑1 (MQ=255)
                                                                                       aCGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGcc                                                                                                                                                                            >  2:518619/1‑274 (MQ=255)
                                                                                                       ttACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCa                                                                                                                                                           >  2:416774/1‑275 (MQ=255)
                                                                                                                             aaCGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACAtgtg                                                                                                                                                                                                                                                                    <  7:767920/148‑1 (MQ=255)
                                                                                                                              aCGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGCGAGAAATCCCCGCCATCGAGCATACAAATCAGCCGTCCAATAATTTCCTGAatc                                                                                                                                    >  2:402466/1‑275 (MQ=255)
                                                                                                                                               ccAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACg                                                                                                                                                                                                                                                                              >  8:527132/1‑120 (MQ=255)
                                                                                                                                               ccAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACg                                                                                                                                                                                                                                                                              <  7:527132/120‑1 (MQ=255)
                                                                                                                                                             tAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGtga                                                                                                                                                                                                                                                                   >  5:399396/1‑117 (MQ=255)
                                                                                                                                                             tAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGtga                                                                                                                                                                                                                                                                   <  6:399396/117‑1 (MQ=255)
                                                                                                                                                                      gATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGaa                                                                                            <  4:416773/275‑1 (MQ=255)
                                                                                                                                                                             tCACAAATCATGCCAGTTTTTCATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGAACGGTGTa                                                                                     <  2:532992/275‑1 (MQ=255)
                                                                                                                                                                                caAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGtgatg                                                                                                                                                                                                                                                                 >  5:912384/1‑100 (MQ=255)
                                                                                                                                                                                caAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGtgatg                                                                                                                                                                                                                                                                 <  6:912384/100‑1 (MQ=255)
                                                                                                                                                                                              ttttAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGa                                                                                                                                                                                                        >  5:135048/1‑143 (MQ=255)
                                                                                                                                                                                              ttttAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGa                                                                                                                                                                                                        <  6:135048/143‑1 (MQ=255)
                                                                                                                                                                                                          gTTATTTAATGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTAGCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTgag                                                                                                                                                                                       >  7:702450/1‑148 (MQ=255)
                                                                                                                                                                                                          gTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGa                                                                                                                                                                                                                                                 <  6:354922/90‑1 (MQ=255)
                                                                                                                                                                                                          gTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGa                                                                                                                                                                                                                                                 >  5:354922/1‑90 (MQ=255)
                                                                                                                                                                                                               ttAAGGAACTTAATCTCTGTAATGCAGGAAAAACGATGTCATCGACACTACTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTTGTGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATc                                                                                                                                                                                  <  8:702450/148‑1 (MQ=255)
                                                                                                                                                                                                                   ggAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATaaaaa                                                                                                                                                              >  4:387518/1‑164 (MQ=255)
                                                                                                                                                                                                                     aaaTTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCgcgg                                                                                                                                                                                                                           <  7:434675/101‑1 (MQ=255)
                                                                                                                                                                                                                     aaaTTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCgcgg                                                                                                                                                                                                                           >  8:434675/1‑101 (MQ=255)
                                                                                                                                                                                                                     aaaTTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGCAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTAAGTCNCCGGAATTTTTAAGCTGGAGAGCGCGGGAAGGGGGTAGGGGGCAATACGGAAAACATCCAATCCGaacaagaaaaac                                             >  3:119377/1‑264 (MQ=255)
                                                                                                                                                                                                                       attaattTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAg                                                                                                                                                          <  2:387519/164‑1 (MQ=255)
                                                                                                                                                                                                                        ttaattTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCATTCACCGGAATTTTTACGCTGGAGAGCGCCGGAACGGTGTAGGGGGCAATCGCGATATCATCGAATCCGATAACTGACACCTg                                           >  1:38294/1‑274 (MQ=255)
                                                                                                                                                                                                                         taattTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAtga                                                                                                                                                                                                                                      <  6:84418/86‑1 (MQ=255)
                                                                                                                                                                                                                         taattTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAtga                                                                                                                                                                                                                                      >  5:84418/1‑86 (MQ=255)
                                                                                                                                                                                                                             ttCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGa                                                                                                                                                                                                                                             >  6:835974/1‑75 (MQ=255)
                                                                                                                                                                                                                             ttCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGa                                                                                                                                                                                                                                             <  5:835974/75‑1 (MQ=255)
                                                                                                                                                                                                                              tCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGATAGCGCCGGAACGGTGTAGGGGGCAATAGCGATATCATCGAATCCGATAACTGACACCTGCCCTgg                                     >  2:84354/1‑274 (MQ=255)
                                                                                                                                                                                                                                 gTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGaa                                                                                                                                                                                     <  5:288681/127‑1 (MQ=255)
                                                                                                                                                                                                                                 gTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGaa                                                                                                                                                                                     >  6:288681/1‑127 (MQ=255)
                                                                                                                                                                                                                                  tAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGa                                                                                                                                                                                                        >  5:525959/1‑107 (MQ=255)
                                                                                                                                                                                                                                  tAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGa                                                                                                                                                                                                        <  6:525959/107‑1 (MQ=255)
                                                                                                                                                                                                                                     tGCAGGAAAAACGATGTCATCGACACTAGTGACGATGAc                                                                                                                                                                                                                                                                         >  4:426889/1‑39 (MQ=255)
                                                                                                                                                                                                                                         ggAAAAACGATGTCATCGACACTAGTGACGATGACAtgt                                                                                                                                                                                                                                                                     <  2:426890/39‑1 (MQ=255)
                                                                                                                                                                                                                                          gAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCgcggcg                                                                                                                                                                                                                         <  8:1002546/82‑1 (MQ=255)
                                                                                                                                                                                                                                          gAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCgcggcg                                                                                                                                                                                                                         >  7:1002546/1‑82 (MQ=255)
                                                                                                                                                                                                                                           aaaaaCGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGaa                                                                                                                                                                                                       <  7:521108/99‑1 (MQ=255)
                                                                                                                                                                                                                                           aaaaaCGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGaa                                                                                                                                                                                                       >  8:521108/1‑99 (MQ=255)
                                                                                                                                                                                                                                             aaaCGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTcc                                                                                                                                                    >  5:755006/1‑148 (MQ=255)
                                                                                                                                                                                                                                                cGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAtgag                                                                                                                                                                                                                                     >  6:164627/1‑64 (MQ=255)
                                                                                                                                                                                                                                                cGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAtgag                                                                                                                                                                                                                                     <  5:164627/64‑1 (MQ=255)
                                                                                                                                                                                                                                                 gATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATaaaaa                                                                                                                                                              >  7:805410/1‑134 (MQ=255)
                                                                                                                                                                                                                                                 gATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATaaaaa                                                                                                                                                              <  8:805410/134‑1 (MQ=255)
                                                                                                                                                                                                                                                    gTCATCGACACTATTGACGATGACATGTGATGACAATGTTTATCGCGACGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCACCGAGCATATAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGAACGGTGTAGGGGGCACTAGCGAAATCATCGAATCCGTTAACTGAAACCTGCTCTGGCACCCCTACGCCGCGCACgtgt               >  4:383872/1‑274 (MQ=255)
                                                                                                                                                                                                                                                      cATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCa                                                                                                                                                                                                                                          <  5:16217/53‑1 (MQ=255)
                                                                                                                                                                                                                                                      cATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCa                                                                                                                                                                                                                                          >  6:16217/1‑53 (MQ=255)
                                                                                                                                                                                                                                                       aTCGTCACTAGTGTCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTcc                                                                                                                                          >  7:387550/1‑148 (MQ=255)
                                                                                                                                                                                                                                                        tCGTCACTAGTGTCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAtgagtgag                                                                                                                                                                                                                                 >  8:106403/1‑60 (MQ=38)
                                                                                                                                                                                                                                                        tCGTCACTAGTGTCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGa                                                                                                                                                                                                        <  5:125757/85‑1 (MQ=255)
                                                                                                                                                                                                                                                        tCGTCACTAGTGTCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGa                                                                                                                                                                                                        >  6:125757/1‑85 (MQ=255)
                                                                                                                                                                                                                                                         cGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTg                                                                                                                                        >  5:250783/1‑148 (MQ=255)
                                                                                                                                                                                                                                                          gACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCaa                                                                                                                                                  <  6:530484/137‑1 (MQ=255)
                                                                                                                                                                                                                                                          gACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCaa                                                                                                                                                  >  5:530484/1‑137 (MQ=255)
                                                                                                                                                                                                                                                            cacTAGTGTCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGa                                                                                                                                                                                                        <  6:551226/81‑1 (MQ=255)
                                                                                                                                                                                                                                                            cacTAGTGTCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGa                                                                                                                                                                                                        >  5:551226/1‑81 (MQ=255)
                                                                                                                                                                                                                                                            cacTAGTGTCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTg                                                                                                                                                                                         >  6:288783/1‑96 (MQ=255)
                                                                                                                                                                                                                                                            cacTAGTGTCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTg                                                                                                                                                                                         <  5:288783/96‑1 (MQ=255)
                                                                                                                                                                                                                                                             acTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAatc                                                                                                                                    >  5:584647/1‑148 (MQ=255)
                                                                                                                                                                                                                                                             acTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGAACGGTGTAGGGGGCAATAGCGATATCATCGAATCCGATAACTGACACCTGCTCTGGCACCGCTACGCCGCGCTCGTGTAACGCTTTc      >  2:530426/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                aGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTcggcgg                                                                                                                                                                                           >  5:495153/1‑90 (MQ=255)
                                                                                                                                                                                                                                                                aGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTcggcgg                                                                                                                                                                                           <  6:495153/90‑1 (MQ=255)
                                                                                                                                                                                                                                                                aGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAatcatc                                                                                                                                 >  8:354812/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                 gTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGCACGGTGTAGGGGGCAATAGCGATCTCATCGAATCCGATAACTAACACCTGCCTTGGCCCCGCTACGCCCCGCCCGTTTAACGCTTTCATCg  >  4:479973/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                  tGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCAtctc                                                                                                                               <  6:250783/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                   gACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGt                                                                                                                                                                                                                  <  6:394908/64‑1 (MQ=255)
                                                                                                                                                                                                                                                                   gACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGt                                                                                                                                                                                                                  >  5:394908/1‑64 (MQ=255)
                                                                                                                                                                                                                                                                   gACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGtt                                                                                                                                                                                                   <  5:974589/79‑1 (MQ=255)
                                                                                                                                                                                                                                                                   gACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGtt                                                                                                                                                                                                   >  6:974589/1‑79 (MQ=255)
                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                              
ATGGATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGAACGGTGTAGGGGGCAATAGCGATATCATCGAATCCGATAACTGACACCTGCTCTGGCACCGCTACGCCGCGCTCGTGTAACGCTTTCATCG  >  NC_000913/2837971‑2838501

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