Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,656,397 A→G 5.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 5.6% 342.2 / 6.2 106intergenic (‑39/+139)pepB/iscXaminopeptidase B/Fe(2+) donor and activity modulator for cysteine desulfurase
Reads supporting (aligned to +/- strand):  ref base A (56/44);  new base G (4/2);  total (60/46)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

Only 100 of 111 total aligned reads displayed.

CCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGACGCGGAAATTATTCGG  >  NC_000913/2656249‑2656543
                                                                                                                                                    |                                                                                                                                                  
ccGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCa                                                                                                                                                    >  1:437590/1‑149 (MQ=255)
ccGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCa                                                                                                                                                    >  3:644487/1‑149 (MQ=255)
ccGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCa                                                                                                                                                    >  4:794020/1‑149 (MQ=255)
 cGTTCAAATGCAGGGTAATGCCGTCAGTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAg                                                                                                                                                   >  3:906971/1‑149 (MQ=255)
  gTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGc                                                                                                                                                  >  2:1024073/1‑149 (MQ=255)
         tGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGa                                                                                                                                           <  4:38516/149‑1 (MQ=255)
          gCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGAt                                                                                                                                          >  1:551100/1‑149 (MQ=255)
                tAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGTGTAGAGCGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTcc                                                                                                                                    <  2:919736/149‑1 (MQ=255)
                   tGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCt                                                                                                                                 >  1:446104/1‑149 (MQ=255)
                   tGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCAATTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCt                                                                                                                                 >  4:928034/1‑149 (MQ=255)
                       gTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAAt                                                                                                                             >  3:243492/1‑149 (MQ=255)
                          attattAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGAt                                                                                                                          >  1:770232/1‑149 (MQ=255)
                             attaATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTcg                                                                                                                       >  3:1076694/1‑149 (MQ=255)
                             attaATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATACTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTcg                                                                                                                       >  3:1024074/1‑149 (MQ=255)
                              ttaATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTcc                                                                                                                                    <  2:553086/135‑1 (MQ=255)
                              ttaATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTcc                                                                                                                                    >  1:553086/1‑135 (MQ=255)
                                aacgCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAATGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGt                                                                                                                    >  4:354809/4‑149 (MQ=255)
                                 aTGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGtt                                                                                                                   >  1:1086622/1‑149 (MQ=255)
                                    cTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCa                                                                                                                <  3:59766/149‑1 (MQ=255)
                                    cTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCa                                                                                                                <  3:346255/149‑1 (MQ=255)
                                     tGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAg                                                                                                               <  4:222741/149‑1 (MQ=255)
                                     tGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAg                                                                                                               >  2:382377/1‑149 (MQ=255)
                                        aaGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGaa                                                                                                            <  2:642246/149‑1 (MQ=255)
                                           ttGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCCGGAAggc                                                                                                         >  2:759706/1‑149 (MQ=255)
                                           ttGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAggc                                                                                                         >  3:740142/1‑149 (MQ=255)
                                           ttGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAggc                                                                                                         <  3:67358/149‑1 (MQ=255)
                                           ttGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAggc                                                                                                         >  1:250457/1‑149 (MQ=255)
                                           ttGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAggc                                                                                                         >  2:54125/1‑149 (MQ=255)
                                           ttGCTTTTACTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAggc                                                                                                         >  3:996624/1‑149 (MQ=255)
                                            tggtttttCTTCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAggcg                                                                                                        <  1:747647/146‑1 (MQ=255)
                                            tGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAggcg                                                                                                        >  3:510387/1‑149 (MQ=255)
                                             gCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAggcgg                                                                                                       >  2:404969/1‑149 (MQ=255)
                                             gCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAggcgg                                                                                                       >  3:389948/1‑149 (MQ=255)
                                              cTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAggcggc                                                                                                      >  1:427805/1‑149 (MQ=255)
                                               tttttCTCCCCAGCGCGCGTCGGCAGGTTGGGTTGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCa                                                                                                     >  3:879331/1‑149 (MQ=255)
                                                  ttCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAagt                                                                                                  <  1:470699/149‑1 (MQ=255)
                                                     tcCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAagtgag                                                                                               <  1:404969/149‑1 (MQ=255)
                                                     tcCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGCTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAagtgag                                                                                               <  3:882493/149‑1 (MQ=255)
                                                      ccccAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAagtgagtg                                                                                             >  2:78652/1‑149 (MQ=255)
                                                       cccAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAagtgagtg                                                                                             <  3:1036355/149‑1 (MQ=255)
                                                       cccAGCGCGCGTCGGCAGGTTGAGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAATGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGt                                                                                                                    >  4:1014131/1‑126 (MQ=255)
                                                       cccAGCGCGCGTCGGCAGGTTGAGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGt                                                                                                                    <  3:1014131/126‑1 (MQ=255)
                                                        ccggcgcgcGTCGGCAGGTTGGGTAGAGAGGGTACTCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAgagagcga                                                                                            >  3:719462/4‑149 (MQ=255)
                                                        ccAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAgtgagtga                                                                                            <  2:427805/149‑1 (MQ=255)
                                                                   cGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAg                                                                                 <  4:381652/149‑1 (MQ=255)
                                                                     gCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCa                                                                               <  2:1094629/149‑1 (MQ=255)
                                                                     gCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCa                                                                               >  2:461923/1‑149 (MQ=255)
                                                                     gCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCa                                                                               >  1:27688/1‑149 (MQ=255)
                                                                      cAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCAt                                                                              <  4:412175/149‑1 (MQ=255)
                                                                       aGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATa                                                                             <  3:794020/149‑1 (MQ=255)
                                                                          ttGGGTTGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGAt                                                                          <  4:879331/149‑1 (MQ=255)
                                                                          ttGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGAt                                                                          >  1:172491/1‑149 (MQ=255)
                                                                              gTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTcgc                                                                                                                      >  4:291363/1‑101 (MQ=255)
                                                                              gTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGACGGCGGCAagt                                                                                                  >  4:84593/1‑121 (MQ=255)
                                                                              gTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAagt                                                                                                  <  3:84593/121‑1 (MQ=255)
                                                                                agagagGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTAt                                                                    <  3:463542/149‑1 (MQ=255)
                                                                                 gagagGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTAtg                                                                   >  1:376378/1‑149 (MQ=255)
                                                                                    agGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGa                                                                <  2:317942/149‑1 (MQ=255)
                                                                                    agGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGa                                                                <  4:996624/149‑1 (MQ=255)
                                                                                     gggTAAGCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGCCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCATGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGAc                                                               <  2:437590/149‑1 (MQ=255)
                                                                                      ggTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACt                                                              >  3:555721/1‑149 (MQ=255)
                                                                                       gTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTg                                                             >  2:836806/1‑149 (MQ=255)
                                                                                        tAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTgg                                                            <  2:782607/149‑1 (MQ=255)
                                                                                              tCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGaa                                                      >  2:896092/1‑149 (MQ=255)
                                                                                                aTCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAACCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACg                                                    >  4:870945/1‑149 (MQ=255)
                                                                                                 tCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGa                                                   <  4:969522/149‑1 (MQ=255)
                                                                                                 tCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGa                                                   <  1:38484/149‑1 (MQ=255)
                                                                                                   gCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAgc                                                 >  4:408409/1‑149 (MQ=255)
                                                                                                       cTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAg                                             <  2:172491/149‑1 (MQ=255)
                                                                                                         gTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCa                                                                                                                                                    <  1:770652/44‑1 (MQ=255)
                                                                                                         gTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCa                                                                                                                                                    >  2:770652/1‑44 (MQ=255)
                                                                                                          tCATTTTAGTTATCGTGCTTGGTAAAACGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTCACTGATTCGCGTTGCACGAAGGCGGTAAGTGAGTGACGCCCCAGAAGAATAGGTAACTATGTGACGTGGGTGAACGAGCGCAGCCa                                          >  2:790648/1‑149 (MQ=255)
                                                                                                             ttttAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACg                                       >  4:727136/1‑149 (MQ=255)
                                                                                                             ttttAGTTATCCTGCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACg                                       <  1:280035/149‑1 (MQ=255)
                                                                                                                    tATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTgc                                <  4:510387/149‑1 (MQ=255)
                                                                                                                     aTCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTgcg                               >  1:966501/1‑149 (MQ=255)
                                                                                                                         ttcttTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGc                                                                                         <  4:397212/87‑1 (MQ=255)
                                                                                                                         ttcttTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGc                                                                                         >  3:397212/1‑87 (MQ=255)
                                                                                                                         ttcttTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCAt                                                                              <  3:907154/98‑1 (MQ=255)
                                                                                                                         ttcttTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCAt                                                                              >  4:907154/1‑98 (MQ=255)
                                                                                                                         ttcttTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGt                           >  2:978920/1‑149 (MQ=255)
                                                                                                                          tcttTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTg                          >  3:849932/1‑149 (MQ=255)
                                                                                                                            tttGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGaa                        <  2:10207/149‑1 (MQ=255)
                                                                                                                            tttGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGaa                        <  4:1076694/149‑1 (MQ=255)
                                                                                                                               gTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCa                     <  1:464582/149‑1 (MQ=255)
                                                                                                                                tAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCAt                    <  1:15171/149‑1 (MQ=255)
                                                                                                                                  ctggggCTGCCAACTGGCGGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCTTGa                  >  3:436674/4‑149 (MQ=37)
                                                                                                                                  aaaGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGa                  >  4:126679/1‑149 (MQ=255)
                                                                                                                                  aaaGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGa                  >  3:196207/1‑149 (MQ=255)
                                                                                                                                   tggggCTGCCAACTGGCGGCCCCAGGATTTTTCCGCTAAGTGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTAt                                                                    <  4:183891/96‑1 (MQ=255)
                                                                                                                                   tggggCTGCCAACTGGCGGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTAt                                                                    >  3:183891/3‑98 (MQ=38)
                                                                                                                                   aaGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGCGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGAc                 >  1:356221/1‑149 (MQ=255)
                                                                                                                                     gggCTGCCAACTGGCGGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGAcgc               >  2:945980/1‑149 (MQ=255)
                                                                                                                                     gggCTGCCAACTGGCGGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGAcgc               >  3:368385/1‑149 (MQ=255)
                                                                                                                                        ccgccAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGACGCGGa            >  4:117807/1‑149 (MQ=255)
                                                                                                                                         cgccAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGACGCGGaa           >  2:837069/1‑149 (MQ=255)
                                                                                                                                          gccAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGACGCGGaaa          >  1:1040904/1‑149 (MQ=255)
                                                                                                                                              gTTGGCAGCCCCAGGATTTTTCCGCTAAATGACTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGACGCGGAAattat      >  2:70627/1‑149 (MQ=255)
                                                                                                                                                tGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGACGCGGAAATTATTc    <  1:1092977/149‑1 (MQ=255)
                                                                                                                                                  gCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGACGCGGAAATTATTCgg  >  4:316597/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGACGCGGAAATTATTCGG  >  NC_000913/2656249‑2656543

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