Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,361,854 A→T 6.0% S120T (TCA→ACA)  rhmR ← putative DNA‑binding transcriptional regulator for the rhm operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,361,8540AT 6.0% 212.9 / 3.9 67S120T (TCA→ACA) rhmRputative DNA‑binding transcriptional regulator for the rhm operon
Reads supporting (aligned to +/- strand):  ref base A (29/34);  new base T (2/2);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.16e-01

TTGTGGGTGTCGCCTGCTCCCAGACCAGCCCTTCGATAATACTTTGCTGTACCGCTGCAGGTTGCCAGGCAAGCAGGCATTTACCTATTCCGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCAATAATCCCCAGGTGACATAACAGCCCGGTGGTATCCATCAGTTGCGTCAAACGCGGTCGAGCCAGTTCCCGGAGATCCATCTTGCTTAACGCATGACCGGAAAGCT  >  NC_000913/2361707‑2362000
                                                                                                                                                   |                                                                                                                                                  
ttGTGGGTGTCGCCTGCTCCCAGACCAGCCCTTCGATAATACTTTGCTGTACCGCTGCAGGTTGCCAGGCAAGCAGGCATTTACCTATTCCGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCtgat                                                                                                                                                   <  2:206425/149‑1 (MQ=255)
ttGTGGGTGTCGCCTGCTCCCAGACCAGCCCTTCGATAATACTTTGCTGTACCGCTGCAGGTTGCAAGGCAAGCAGGCATTTACCTATTCCGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCtgat                                                                                                                                                   >  3:370577/1‑149 (MQ=255)
      gtgtCGCCTGCTCCCAGACCAGCCCTTCGATAATACTTTGCTGTACCGCTGCAGGTTGCCAGGCAAGCAGGCATTTACCTATTCCGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTc                                                                                                                                              >  4:104596/1‑148 (MQ=255)
           gCCTGCTCCCAGACCAGCCCTTCGATAATACTTTGCTGTACCGCTGCAGGTTGCCAGGCCAGCAGGCATGTACCTATTCCGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCCCGCTGATAGTGGCTGATGACTCAACttt                                                                                                                                        >  1:185223/1‑149 (MQ=255)
             cTGCTCCCAGACCAGACCTTCGATAATACTTTGCTGTACCGCTGCAGGTTGCCAGGCAAGCAGGCATTTACCTATTCCGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTc                                                                                                                                       >  4:34115/1‑148 (MQ=255)
                 tCCCAGACCAGCCCTTCGATAATACTTTGCTGTACCGCTGCAGGTTGCCAGGCAAGCAGGCATTTACCTATTCCGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAAt                                                                                                                                  >  6:125435/1‑149 (MQ=255)
                    cAGACCAGCCCTTCGATAATACTTTGCTGTACCGCTGCAGGTTGCCAGGCAAGCAGGCATTTACCTATTCCGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGt                                                                                                                                >  8:377207/1‑148 (MQ=255)
                     aGACCAGCCCTTCGATAATACTTTGCTGTACCGCTGCAGGTTGCCAGGCAAGCAGGCATTTACCTATTCCGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGtaa                                                                                                                              <  3:91674/149‑1 (MQ=255)
                      gACCAGCCCTTCGTTAATACTTTGCTGTAGCTCTTCAGGTTGCCTGGCAAGCAGTCATTTACCTATTCCGGAGCGATAAATCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGtaat                                                                                                                             <  6:64675/149‑1 (MQ=255)
                      gACCAGCCATGCGATAATACTTTGCTGTCCCGCTGCAGGTTGCCAGGCAAGAAGGCATTTACCTATTCCGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACTCTGATAGTGGCTGATGACTCCACTTTCAAAATGtaat                                                                                                                             <  6:70069/149‑1 (MQ=255)
                          aGCCCTTCGATAATACTTTGCTGTACCGCTGCAGGTTGCCAGGCAAGCAGGCATTTACCTATTCCGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAgc                                                                                                                         <  5:125435/149‑1 (MQ=255)
                                  gATAATACTTTGCTGTACCGCTGCAGGTTGCCAGGCAAGCAGGCATTTACCTATTCCGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGt                                                                                                                 >  4:165995/1‑149 (MQ=255)
                                    taataCTTTGCTGTACCGCTGCAGGTTGCCAGGCAAGCAGGCATTTACCTATTCCGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTa                                                                                                               <  1:347933/149‑1 (MQ=255)
                                      ataCTTTGCTGTACCGCTGCAGGTTGCCAGGCAATCAGGCATTTACCTATTCCGGAGCGATAAAGCTAAAGGCTTTTTCCTTCATTGGAACGCACGCTGATAGTGGCTGATGACTCAACTTTCAAAATGTAATAAGCGCTGCCGTTAt                                                                                                              <  5:123128/148‑1 (MQ=255)
                                      ataCTTTGCTGTACCGCTGCAGGTTGCCAGGCAAGCAGGCATTTACCTATTCCGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTAt                                                                                                              >  7:240019/1‑148 (MQ=255)
                                        aCTTTGCTGTACCGCTGCAGGTTGCCAGGCAAGCAGGCATTTACCTATTCCGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGTTGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCa                                                                                                            >  7:262862/1‑148 (MQ=255)
                                        aCTTTGCTGTACCGCTGCAGGTTGCCAGGCAAGCAGGCATTTACCTATTCCGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGTTGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCa                                                                                                            >  7:256897/1‑148 (MQ=255)
                                        aCTTTGCTGTACCGCTGCAGGTTGCCAGGCAAGCAGGCATTTACCTATTCCGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGTTGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCa                                                                                                            <  8:262862/148‑1 (MQ=255)
                                        aCTTTGCTGTACCGCTGCAGGTTGCCAGGCAAGCAGGCATTTACCTATTCCGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGTTGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCa                                                                                                            <  8:256897/148‑1 (MQ=255)
                                         cTTTGCTGTACCGCTGCAGGTTGCCAGGCAAGCAGGCATTTACCTATTCCGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCaa                                                                                                           <  6:236771/148‑1 (MQ=255)
                                         cTTTGCTGTACCGCTGCAGGTTGCCAGGCAAGCAGGCATTTACCTATTCCGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCaa                                                                                                           <  7:377207/148‑1 (MQ=255)
                                            tGCTATACCGCTGCAGGTTGCCAGGCAAGCAGGCATTTACCTATTCCGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCaataa                                                                                                        <  6:212663/148‑1 (MQ=255)
                                              cTGTACCGCTGCAGGTTGCCAGGCAAGCAGGCATTTACCTATTCCGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCAATAATc                                                                                                      <  6:222086/148‑1 (MQ=255)
                                               tGTACCGCTGCAGGTTGCCAGGCAAGCAGGCATTTACCTATTCCGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCa                                                                                                            <  8:306317/141‑1 (MQ=255)
                                               tGTACCGCTGCAGGTTGCCAGGCAAGCAGGCATTTACCTATTCCGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCa                                                                                                            >  7:306317/1‑141 (MQ=255)
                                               tGTACCGCTGCAGGTTGCCAGGCAAGCAGGCATTTACCTATTCCGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCAATAATccc                                                                                                    >  1:119750/1‑149 (MQ=255)
                                                 tacggatGCAGGTTGCCAGGCATGCAGGGATTTACCTATTCCGGAGCGATAAAGCGAAAGTGTTTTTCGTTAATGGGAACGTACGCTGTTAGTGGGTGATGACTGCACTTTCTAAATGTATTAAGCGCTGCGGTTATAAATAATccc                                                                                                    <  8:240019/141‑1 (MQ=255)
                                                 tACCGCTGCAGGTTGCCAGGCAAGCAGGCATTTACCTATTCCGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCAATAATcccc                                                                                                   <  2:119750/148‑1 (MQ=255)
                                                 tACCGCTGCAGGTTGCCAGGCAAGCAGGCATTTACCTATTCCGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCAATAATccc                                                                                                    <  7:184705/147‑1 (MQ=255)
                                                   ccGCTGCAGGTTGCCAGGCAAGCAGGCATTTACCTATTCCGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCAATAATCCCCAgg                                                                                                >  5:244493/1‑149 (MQ=255)
                                                    cGCTGCAGGTTACCAGGCAAGCAGGCATTTACCTATTCCGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGtaa                                                                                                                              <  6:388129/118‑1 (MQ=255)
                                                    cGCTGCAGGTTACCAGGCAAGCAGGCATTTACCTATTCCGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGtaa                                                                                                                              >  5:388129/1‑118 (MQ=255)
                                                     gCTGCAGGTTGCCAGGCAAGCAGGCAATTACGTATTCCGGAGCGATAAAGCGATAGGCTTTTTCCTTCATGGGAACGCACGGTGATAGTGTCTGATGACTCCCGTTTCAAAATGTAATAAGCGCTGCCGTTATAAATAATCCCCAGGt                                                                                               >  5:245674/1‑148 (MQ=255)
                                                      cTGCAGGTTGCCAGGCAAGCAGGCATTTACCTATTCCGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCAATAATCCACAGGTGa                                                                                             >  6:227787/1‑149 (MQ=255)
                                                         cAGGTTTCCAGGCAAGCAGGCATTTACCTATTCCGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGCGAACGCACGCTGATAGAGGCTGATGACTCCACTTTCAAAATGTAAATAGCGCTGCAGTTATCAATAATCCCCAGGGGAc                                                                                            <  2:150579/147‑1 (MQ=255)
                                                         cAGGTTGCCAGGCAAGCAGGCATTTACCTATTCCGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGt                                                                                                                                >  2:224503/1‑111 (MQ=255)
                                                             tggcaaGGCTAGCAGGCATTTACCTATTCCGGAGCGATAAAGCGAAAGCCTTTTTCCTTCATGGGAACGCCCGCTGATTTTGGCTGATGACTCCACTTTCAAAATCTAATAAGCGCTGGCTTTATCAATTATGCCAACGTGAGTTAACa                                                                                      <  4:319055/144‑1 (MQ=255)
                                                               gCCAGGCAAGCAGTCATTTACATATTCCGGATCGATTAAGCGAAAGGCTTTTTCCTTCATGGGAACGGACGCTGATAGTTGCTGATGACTCCACTTTCAACATGTAATAAGCTCTGCCGTTATCAATAATCCCCAGGTGACATAACAg                                                                                     <  6:245674/148‑1 (MQ=255)
                                                                     cAAGCAGGCATTTACCTATTCCGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCAATAATCCCCAGGTGACATAACAGCCCggtg                                                                              <  7:148639/149‑1 (MQ=255)
                                                                         cAGGCATTTACCTATTCCGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCAATAATCCCCAGGTGACATAACAGCCCGGTGGTAt                                                                          >  4:277803/1‑149 (MQ=255)
                                                                               tttACCTATTCCGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCAATAATCCCCAGGTGACATAACAGCCCGGTGGTATCCATCa                                                                    <  3:34115/149‑1 (MQ=255)
                                                                                 tACCTATTCCGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCAATAATCCCCAGGTGACATAACAGCCCGGTGGTATCCATCAGt                                                                  >  8:198459/1‑149 (MQ=255)
                                                                                   ccTATTCCGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCAATAATCCCCAGGTGACATAACAGCCCGGTGGTATCCATCAGtt                                                                 <  7:85407/148‑1 (MQ=255)
                                                                                          cGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCAATAATCCCCAg                                                                                                 <  8:221975/109‑1 (MQ=255)
                                                                                          cGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCAATAATCCCCAg                                                                                                 >  7:221975/1‑109 (MQ=255)
                                                                                          cGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCAATAATCCCCAGGTGACATAACAGCCCGGTGGTATCCATCAGTTGCGTCaac                                                         >  1:246784/1‑148 (MQ=255)
                                                                                            gAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCAATAATCCCCAGGTGACATAACAGCCCGGTGGTATCCATCAGTTGCGTCAAAcg                                                       >  7:202551/1‑149 (MQ=255)
                                                                                              gCGATAAAGCGAAGGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCAATAATCCCCAGGTGACATAACAGCCCGGTGGTATCCATCAGTTGCGTCAAAcgcg                                                     >  4:43466/1‑149 (MQ=255)
                                                                                                    aaGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCAATAATCCCCAGGTGACATAACAGCCCGGTGGTATCCATCAGTTGCGTCAAACGCGGTCGAg                                               >  3:203739/1‑149 (MQ=255)
                                                                                                      gCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGGTATCAATAATCGCCAGGTGACATAACAGCCCGGTGGTATCCAGCAGTTGCGTCAAACGCGGTCGAGcc                                             >  2:292659/1‑149 (MQ=255)
                                                                                                        gAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCAATAATCCCCAGGTGACATAACAGCCCGGTGGTATCCATCAGTTGCGTCAAACGCGGTCGAGCCAg                                           <  3:104596/149‑1 (MQ=255)
                                                                                                                ttttCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCAATAATCCCCAGGTGACATAACAGCCCGGTGGTATCCATCAGTTGCGTCAAACGCGGTCGAGCCAGTTCCCg                                     >  3:175946/1‑147 (MQ=255)
                                                                                                                 tttCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCAATAATCCCCAGGTGACATAACAGCCCGGTGGTATCCATCAGTTGCGTCAAACGCAGTCGAGCCAGTTCCCg                                     >  1:174743/1‑146 (MQ=255)
                                                                                                                  ttCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCAATAATCCCCAGGTGACATAACAGCCCGGTGGTATCCATCAGTTGCGTCAAACGCGGTCGAGCCAGTTCCCGgag                                  <  8:202551/148‑1 (MQ=255)
                                                                                                                   tCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCAATAATCCCCAGGTGACATAACAGCCCGGTGGTATCCATCAGTTGCGTCAAACGCGGTCGAGCCAGTTCCCGgaga                                 <  3:43466/148‑1 (MQ=255)
                                                                                                                     cTTCATTGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCAATAATCCCCAGGTGACATaa                                                                                        <  8:135584/91‑1 (MQ=255)
                                                                                                                     cTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCAATAATCCCCAGGTGACATAACAGCCCGGTGGTATCCATCAGTTGCGTCAAACGCGGTCGAGCCAGTTCCCGGAGATc                               <  7:198459/148‑1 (MQ=255)
                                                                                                                     cTTCATAGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCAATAATCCCCAGGTGACATaa                                                                                        >  7:135584/1‑91 (MQ=255)
                                                                                                                        cATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCAATAATCCCCAGGTGACATAACAGCCCggtg                                                                              <  2:176603/98‑1 (MQ=255)
                                                                                                                        cATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCAATAATCCCCAGGTGACATAACAGCCCggtg                                                                              >  1:176603/1‑98 (MQ=255)
                                                                                                                        cATGCGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCAATAATCCCCAGGTGACATAACAGCCCGGTGGTATCCATCAGTTGCGTCAAACGCGGTCGAGCCAGTTCCCGGAGATCCATc                           >  6:107255/1‑149 (MQ=255)
                                                                                                                            ggAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCAATAATCCCCAGGTGACATAACAGCCCGGTGGTATCCATCAGTTGCGTCAAACGCGGTCGAGCCAGTTCCCGGAGATCCATCTTg                        <  3:156354/148‑1 (MQ=255)
                                                                                                                               acgcacgcTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCAATAATCCCCAGGTGACATAACAGCCCGGTGGTATCCATCAGTTGCGTCAAACGCGGTCGAGCCAGTTCCCGGAGATCCATCTTGCTTa                    <  3:266493/149‑1 (MQ=255)
                                                                                                                                cgcacgcTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCAATAATCCCCAGGTGACATAACAGCCCGGTGGTATCCATCAGTTGCGTCAAACGCGGTCGAGCCAGTTCCCGGAGATCCATCTTGCTTaa                   <  5:107255/149‑1 (MQ=255)
                                                                                                                                 gcacgcTGATAGTGGCTGATGACTCGACTTTCAAAATGTAATCAGCGCTGCCGTTATCAATAATCCCCAGGTGACATAACAGCCCGGTGGTATCCATCAGTTGCGTCAAACGCGGTCGAGCCAGTTCCCGGAGATCCATCTTGCTTAAc                  <  5:288933/149‑1 (MQ=255)
                                                                                                                                     gcTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCAATAATCCCCACGTGACATAACAGCCCGGTGTTATCCATCAGTTGCGTCAAACGCGGTCGAGCCAGTTCCCGGAGATCCATCTTGCTTAACGCAt              <  4:370577/149‑1 (MQ=255)
                                                                                                                                       tGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCAATAATCCCCAGGTGACATAACAGCCCGGTGGTATCCATCAGTTGCGTCAAACGCGGTCGAGCCAGTTCCCGGAGATCCATCTTGCTTAACGCATGa            >  2:148654/1‑149 (MQ=255)
                                                                                                                                             tGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCAATAATCCCCAGGTGACATAACAGCCCGGTGGTATCCATCAGTTGCGTCAAACGCGGTCGAGCCAGTTCCCGGAGATCCATCTTGCTTAACGCATGACCGGaa      <  3:277803/149‑1 (MQ=255)
                                                                                                                                             tGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCAATAATCCCCAGGTGACATAACAGCCCGGTGGTATCCATCAGTTGCGTCAAACGCGGTCGAGCCAGTTCCCGGAGATCCATCTTGCTTAACGCATGACCGGaa      >  6:29005/1‑149 (MQ=255)
                                                                                                                                                 tgatgaCTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCAATAATCCCCAGGTGACATAACAGCCCGGTGGTATCCATCAGTTGCGTCAAACGCGGTCGAGCCAGTTCCCGGAGATCCATCTTGCTTAACGCATGACCGGAAAGCt  >  7:3888/1‑149 (MQ=255)
                                                                                                                                                 tgatgaCTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCAATAATCCCCAGGTGACATAACAGCCCGGTGGTATCCATCAGTTGCGTCAAACGCGGTCGAGCCAGTTCCCGGAGATCCATCTTGCGTAACGCATGACCGGAAAGCt  >  5:208219/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTGTGGGTGTCGCCTGCTCCCAGACCAGCCCTTCGATAATACTTTGCTGTACCGCTGCAGGTTGCCAGGCAAGCAGGCATTTACCTATTCCGGAGCGATAAAGCGAAAGGCTTTTTCCTTCATGGGAACGCACGCTGATAGTGGCTGATGACTCCACTTTCAAAATGTAATAAGCGCTGCCGTTATCAATAATCCCCAGGTGACATAACAGCCCGGTGGTATCCATCAGTTGCGTCAAACGCGGTCGAGCCAGTTCCCGGAGATCCATCTTGCTTAACGCATGACCGGAAAGCT  >  NC_000913/2361707‑2362000

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