Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,656,386 C→T 5.4% intergenic (‑28/+150) 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,3860CT 5.4% 367.5 / 6.9 111intergenic (‑28/+150)pepB/iscXaminopeptidase B/Fe(2+) donor and activity modulator for cysteine desulfurase
Reads supporting (aligned to +/- strand):  ref base C (59/46);  new base T (4/2);  total (63/48)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

Only 100 of 114 total aligned reads displayed.

GATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGACGCGGAA  >  NC_000913/2656238‑2656534
                                                                                                                                                    |                                                                                                                                                    
gATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGgcc                                                                                                                                                      >  4:114118/1‑149 (MQ=255)
     tCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAg                                                                                                                                                 >  1:10207/1‑149 (MQ=255)
       tGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGtt                                                                                                                                               <  3:360561/149‑1 (MQ=255)
        gCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTg                                                                                                                                              <  3:954306/149‑1 (MQ=255)
          cccGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGc                                                                                                                                            >  1:1029783/1‑149 (MQ=255)
           ccGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCa                                                                                                                                           >  3:644487/1‑149 (MQ=255)
           ccGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCa                                                                                                                                           >  1:437590/1‑149 (MQ=255)
            cGTTCAAATGCAGGGTAATGCCGTCAGTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAg                                                                                                                                          >  3:906971/1‑149 (MQ=255)
             gTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGc                                                                                                                                         >  2:1024073/1‑149 (MQ=255)
             gTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGc                                                                                                                                         <  1:702848/149‑1 (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)
                                  gTCATTATTAATGCTGTAAGTTGCTTTTTATCCCCAGCGCGCGTCGGCAGGTTTGGTGGAGAGGGTAAGCTTCATCGATTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCGCCAGGATTTTTCCGCTAAAt                                                                                                                    >  2:745369/1‑149 (MQ=255)
                                     attattAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGAt                                                                                                                 >  1:770232/1‑149 (MQ=255)
                                        attaATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTcg                                                                                                              >  3:1076694/1‑149 (MQ=255)
                                         ttaATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTcc                                                                                                                           >  1:553086/1‑135 (MQ=255)
                                         ttaATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTcc                                                                                                                           <  2:553086/135‑1 (MQ=255)
                                            aTGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGtt                                                                                                          >  1:1086622/1‑149 (MQ=255)
                                               cTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCa                                                                                                       <  3:346255/149‑1 (MQ=255)
                                               cTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCa                                                                                                       <  3:59766/149‑1 (MQ=255)
                                                tGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAg                                                                                                      >  2:382377/1‑149 (MQ=255)
                                                tGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAg                                                                                                      <  4:222741/149‑1 (MQ=255)
                                                      ttGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCCGGAAggc                                                                                                >  2:759706/1‑149 (MQ=255)
                                                      ttGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAggc                                                                                                >  3:740142/1‑149 (MQ=255)
                                                      ttGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAggc                                                                                                >  2:54125/1‑149 (MQ=255)
                                                      ttGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAggc                                                                                                <  3:67358/149‑1 (MQ=255)
                                                      ttGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAggc                                                                                                >  1:250457/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                                                                                              >  3:389948/1‑149 (MQ=255)
                                                        gCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAggcgg                                                                                              >  2:404969/1‑149 (MQ=255)
                                                         cTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAggcggc                                                                                             >  1:427805/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)
                                                                 ccccAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAagtgagt                                                                                     <  1:382377/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)
                                                                   ccAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTgg                                                                                                                                             <  1:975760/91‑1 (MQ=255)
                                                                   ccAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTgg                                                                                                                                             >  2:975760/1‑91 (MQ=255)
                                                                   ccAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAgtgagtga                                                                                   <  2:427805/149‑1 (MQ=255)
                                                                        gcgcgTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGccc                                                                              <  3:887635/149‑1 (MQ=255)
                                                                              cGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAg                                                                        <  4:381652/149‑1 (MQ=255)
                                                                                gCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCa                                                                      >  2:461923/1‑149 (MQ=255)
                                                                                gCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCa                                                                      <  2:1094629/149‑1 (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)
                                                                                         gTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTcgc                                                                                                             <  3:291363/101‑1 (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                                                                                                                                           >  2:770652/1‑44 (MQ=255)
                                                                                                                    gTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCa                                                                                                                                           <  1:770652/44‑1 (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                                                                                >  3:397212/1‑87 (MQ=255)
                                                                                                                                    ttcttTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGc                                                                                <  4:397212/87‑1 (MQ=255)
                                                                                                                                    ttcttTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCAt                                                                     <  3:907154/98‑1 (MQ=255)
                                                                                                                                    ttcttTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCAt                                                                     >  4:907154/1‑98 (MQ=255)
                                                                                                                                    ttcttTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGt                  >  4:249908/1‑149 (MQ=255)
                                                                                                                                    ttcttTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGt                  >  2:978920/1‑149 (MQ=255)
                                                                                                                                     tcttTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTg                 >  3:849932/1‑149 (MQ=255)
                                                                                                                                       tttGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGaa               <  2:10207/149‑1 (MQ=255)
                                                                                                                                           tAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCAt           <  1:15171/149‑1 (MQ=255)
                                                                                                                                             ctggggCTGCCAACTGGCGGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCTTGa         >  3:436674/4‑149 (MQ=37)
                                                                                                                                             aaaGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGa         >  3:196207/1‑149 (MQ=255)
                                                                                                                                             aaaGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGa         >  4:126679/1‑149 (MQ=255)
                                                                                                                                              tggggCTGCCAACTGGCGGCCCCAGGATTTTTCCGCTAAGTGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTAt                                                           <  4:183891/96‑1 (MQ=255)
                                                                                                                                              tggggCTGCCAACTGGCGGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTAt                                                           >  3:183891/3‑98 (MQ=38)
                                                                                                                                              aaGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGAc        >  2:937748/1‑149 (MQ=255)
                                                                                                                                                gggCTGCCAACTGGCGGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGAcgc      >  2:945980/1‑149 (MQ=255)
                                                                                                                                                   ccgccAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGACGCGGa   >  4:117807/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGACGCGGAA  >  NC_000913/2656238‑2656534

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