Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,656,397 A→G 6.6% intergenic (‑39/+139) pepB ← / ← iscX aminopeptidase B/Fe(2+) donor and activity modulator for cysteine desulfurase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,656,3970AG 6.6% 224.3 / 4.2 71intergenic (‑39/+139)pepB/iscXaminopeptidase B/Fe(2+) donor and activity modulator for cysteine desulfurase
Reads supporting (aligned to +/- strand):  ref base A (37/29);  new base G (2/3);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

GTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGACGCGGAAATTATTCGGCC  >  NC_000913/2656251‑2656545
                                                                                                                                                  |                                                                                                                                                    
gTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGc                                                                                                                                                    >  1:119141/1‑149 (MQ=255)
  tCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCTGCCAACTGGCGGccc                                                                                                                                                  <  1:451122/149‑1 (MQ=255)
  tCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCTGCCAACTGGCGGccc                                                                                                                                                  <  3:440068/149‑1 (MQ=255)
        gCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAg                                                                                                                                               <  1:382019/146‑1 (MQ=255)
        gCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAg                                                                                                                                               >  2:382019/1‑146 (MQ=255)
         cAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGAtt                                                                                                                                           <  3:140904/149‑1 (MQ=255)
          aGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTTATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGAttt                                                                                                                                          <  4:403903/149‑1 (MQ=255)
                    cGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTaaa                                                                                                                                <  1:128872/149‑1 (MQ=255)
                     gTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAAt                                                                                                                               >  2:591650/1‑149 (MQ=255)
                     gTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAAt                                                                                                                               >  2:82055/1‑149 (MQ=255)
                        antattAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGAt                                                                                                                            >  3:113192/1‑149 (MQ=255)
                         ttattaATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGAtt                                                                                                                           <  1:136009/149‑1 (MQ=255)
                               aTGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGtt                                                                                                                     <  2:71013/149‑1 (MQ=255)
                                  cTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCa                                                                                                                  <  1:229490/149‑1 (MQ=255)
                                  cTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGAGTTGCa                                                                                                                  <  2:422647/149‑1 (MQ=255)
                                        gTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAgg                                                                                                            <  2:299373/149‑1 (MQ=255)
                                         ttGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAggc                                                                                                           >  2:104682/1‑149 (MQ=255)
                                         ttGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAggc                                                                                                           <  4:26154/149‑1 (MQ=255)
                                         ttGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAggc                                                                                                           >  4:6455/1‑149 (MQ=255)
                                          tGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAggcg                                                                                                          <  1:553465/149‑1 (MQ=255)
                                            cTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAggcggc                                                                                                        >  3:339435/1‑149 (MQ=255)
                                             tttttCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCa                                                                                                       >  4:425940/1‑149 (MQ=255)
                                             tttttCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCa                                                                                                       >  3:379543/1‑149 (MQ=255)
                                              ttttCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCaa                                                                                                      <  4:32770/149‑1 (MQ=255)
                                              ttttCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCaa                                                                                                      <  3:126172/149‑1 (MQ=255)
                                                    ccccAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAagtgagt                                                                                                >  2:99188/1‑149 (MQ=255)
                                                     cccAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAagtgagtg                                                                                               <  1:104682/149‑1 (MQ=255)
                                                        aGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAg                                                                                            <  2:455919/149‑1 (MQ=255)
                                                             gcgTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCa                                                                                       <  4:293877/149‑1 (MQ=255)
                                                              cgTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAg                                                                                      <  1:275567/149‑1 (MQ=255)
                                                                tCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGa                                                                                    <  2:119141/149‑1 (MQ=255)
                                                                   gCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGGGCa                                                                                 >  1:591083/1‑149 (MQ=255)
                                                                     aGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATa                                                                               >  3:34599/1‑149 (MQ=255)
                                                                       gTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGAt                                                                                                                            <  3:258029/102‑1 (MQ=255)
                                                                       gTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGAt                                                                                                                            >  4:258029/1‑102 (MQ=255)
                                                                          gggTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATaa                                                                          <  1:266995/149‑1 (MQ=255)
                                                                            gTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACt                                                                        >  1:445952/1‑149 (MQ=255)
                                                                                    ggTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGccc                                                                                                                                                  <  3:133476/67‑1 (MQ=255)
                                                                                    ggTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGccc                                                                                                                                                  >  4:133476/1‑67 (MQ=255)
                                                                                    ggTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGAAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACt                                                                >  2:211513/1‑149 (MQ=255)
                                                                                     gTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTg                                                               >  4:558628/1‑149 (MQ=255)
                                                                                      tAATCTTCATCGCTTCTGTCATCTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTgg                                                              >  2:287591/1‑149 (MQ=255)
                                                                                       aaTCTTCATCGCTTCTGTCATTCTAGTTATCCTTCTTTGTAAACGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGAGACTggg                                                             >  2:219223/1‑149 (MQ=255)
                                                                                       aaTCTTCATCGCTTCTGTCATCTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTggg                                                             <  1:287591/149‑1 (MQ=255)
                                                                                               tCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGa                                                     >  2:57404/1‑149 (MQ=255)
                                                                                               tCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGa                                                     <  1:57404/149‑1 (MQ=255)
                                                                                                       gTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGcc                                             >  2:309052/1‑149 (MQ=255)
                                                                                                       gTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACAATGTGACTGGGGTGAACGAGCGCAGcc                                             <  3:391730/149‑1 (MQ=255)
                                                                                                           ttttAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACg                                         >  2:302023/1‑149 (MQ=255)
                                                                                                           ttttAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACg                                         >  1:503014/1‑149 (MQ=255)
                                                                                                               aGTTAACCTGCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTAAGTGACTGGGGTGAACGAGCGCAGCCAACGCATc                                     <  1:99188/149‑1 (MQ=255)
                                                                                                                     ccTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTgcggc                               >  3:313553/1‑149 (MQ=255)
                                                                                                                       ttcttTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGt                             >  3:509280/1‑149 (MQ=255)
                                                                                                                       ttcttTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGt                             >  1:169001/1‑149 (MQ=255)
                                                                                                                       ttcttTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGt                             >  1:215872/1‑149 (MQ=255)
                                                                                                                       ttcttTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGt                             >  2:571488/1‑149 (MQ=255)
                                                                                                                          tttGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGaa                          >  2:72/1‑149 (MQ=255)
                                                                                                                           ttGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAg                         >  2:349510/1‑149 (MQ=255)
                                                                                                                             gTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCa                       >  1:496994/1‑149 (MQ=255)
                                                                                                                               aaaaGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGa                                                     >  1:112236/1‑117 (MQ=255)
                                                                                                                               aaaaGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGa                                                     <  2:112236/117‑1 (MQ=255)
                                                                                                                                  ggggCTGCCAACTGGCGGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGAcg                  >  2:45203/2‑149 (MQ=255)
                                                                                                                                   gggCTGCCAACTGGCGGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGaa                          >  2:355809/1‑140 (MQ=255)
                                                                                                                                   gggCTGCCAACTGGCGGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGaa                          <  1:355809/140‑1 (MQ=37)
                                                                                                                                   gggCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGAcgc                 >  4:491486/1‑149 (MQ=255)
                                                                                                                                        gccAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGACGCGGaaa            >  3:575217/1‑149 (MQ=255)
                                                                                                                                            gTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGaa                                                        >  4:256535/1‑101 (MQ=255)
                                                                                                                                            gTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGaa                                                        <  3:256535/101‑1 (MQ=255)
                                                                                                                                              tGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGACGCGGAAATTATTc      <  1:571488/149‑1 (MQ=255)
                                                                                                                                               ggCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGACGCGGAAATTATTCg     <  3:558628/149‑1 (MQ=255)
                                                                                                                                               ggCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGACGCGGAAATTATTCg     >  2:332158/1‑149 (MQ=255)
                                                                                                                                                  aGCCCCAGGATTTTTACGCTAAATGATTCGTGTTGCAGGAACGCGGCAAGTGAGTGAAGCCCCAGGAGCGTAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGACGCGGACATTATTCGGcc  >  4:490049/1‑149 (MQ=38)
                                                                                                                                                  |                                                                                                                                                    
GTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGACGCGGAAATTATTCGGCC  >  NC_000913/2656251‑2656545

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