Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,661,885 C→G V83L (GTT→CTT)  iscR ← DNA‑binding transcriptional dual regulator IscR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,661,8850CG97.4% 115.6 / NA 39V83L (GTT→CTT) iscRDNA‑binding transcriptional dual regulator IscR
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (14/24);  minor base T (0/1);  total (14/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

ATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTC  >  NC_000913/2661642‑2662129
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
attaAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAACCCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAAGGGCGCt                                                                                                                                                                                                                                                <  1:53006/250‑1 (MQ=255)
       gcgTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAAGGGCGCTAATTACtt                                                                                                                                                                                                                                        >  2:41032/1‑251 (MQ=255)
        cgTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAAGGGCGCTAATTACTTc                                                                                                                                                                                                                                       >  1:291585/1‑251 (MQ=255)
                              cTTGTGTGCGGATGGGTGCAGCGTCGTGAGTTTGCTGCCGACTAGAAACATCCAGCACTTCTTGGTTATTAACCAGTTCGCCTAAGGTACTGTTGTTGAGAAGACCGGTGCGACGGTCGCTCAAATCGCGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAAGGGCGCTAATTACTTCGCCACAGGCGATGCTGCTGGcc                                                                                                                                                                                                                 >  2:278440/1‑250 (MQ=255)
                                         gTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGGAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCAGCCGTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAAGGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTaa                                                                                                                                                                                                      <  1:199385/251‑1 (MQ=255)
                                               ggCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAAGGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATa                                                                                                                                                                                                <  1:263348/251‑1 (MQ=255)
                                                    gtcgtGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAAGGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAaccac                                                                                                                                                                                           <  2:101494/251‑1 (MQ=255)
                                                        tGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAAGGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGcc                                                                                                                                                                                       >  1:416505/1‑251 (MQ=255)
                                                           gTATGCTGACGACCAGACACATCCATCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAAGGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCTCCGCCTgg                                                                                                                                                                                    <  2:347792/251‑1 (MQ=255)
                                                             aTGCTGACGACCCGACACATCCAGCCCTTCCTGTATATTAACCAGTTCGCCTAAAGTAATGTTGTTGCGACACCCGGTGAGACGGTCGCTCAACTCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGTCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCCAGGGCGCTAATTACTTCGCCAACGGCGATGCTGCTCGCCTCTTTGCCTAACAGCTAACCACCGCCTGGTc                                                                                                                                                                                  >  2:26183/1‑251 (MQ=255)
                                                              tGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCACGGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATCACCACCGCCTGGTcc                                                                                                                                                                                 >  1:329333/1‑251 (MQ=255)
                                                               gCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGAAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTCCCCTGACAACGGGTGGCATCTACAGATTCGTCAAGGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTcc                                                                                                                                                                                 <  2:328351/250‑1 (MQ=255)
                                                                    cgacgaGACAAATCCAGCACTTCCTGGTTATTAACCATGTCGCCTAAGGGAATGTTGTGGAAAAAACCGGTGAGGCTGTCGCTCAAGTCACGCCACAGCGCGTGGGTAAGGCATTTATCGCCGCCCTGGCCGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGACAAGGTCGCTAATTTCTTCGCCAGCGGCGAGGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTAc                                                                                                                                                                           <  1:7374/251‑1 (MQ=255)
                                                                       ccAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAAGGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCt                                                                                                                                                                        <  2:25517/251‑1 (MQ=255)
                                                                                  aagaaCTTCCTGGTTATGAACCAGTTCGCCTAAAGTAAGGTGGTTGAGAAAACCGGTGAGACGGCCGCTAAAACCACGCCACAGCGCGTGGGTAAGGCATTTATCGCCGCCCTGGCAGCCGCCTGTACCCTGACAACGGGTGTCATCTACAGATTCGTCAAGGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGAc                                                                                                                                                             <  1:14466/247‑1 (MQ=255)
                                                                                          ccTGGTTATTAACCAGTGCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAAGGGCGCGAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTa                                                                                                                                                     <  2:407000/251‑1 (MQ=255)
                                                                                                              caaaaattaaggttgttGAGAAAACCGGTGAGACGTTCGCTCAAATCACGCCCCAGCGCGTGGGTCAGGCATTTACCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAAGGGCGCTAATTCCTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAATCCAGACCATTTTTACGCAGACGGGAAAACAgttg                                                                                                                                 <  2:160894/240‑1 (MQ=255)
                                                                                                                tAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAAGGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAgttgtt                                                                                                                               >  1:313486/1‑251 (MQ=255)
                                                                                                                      aaTGTTGTTGCGAAAACCGGTGAGACGGTCTCTCAAATCCCGCCCCAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTCCCCTGACAACGGGTGGCATCTACAGATTCGTCAAGGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGTAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGag                                                                                                                         <  2:2915/251‑2 (MQ=255)
                                                                                                                      aaTGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAAGGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGa                                                                                                                          <  1:107313/250‑1 (MQ=255)
                                                                                                                           tgttgAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAAGGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCGTTGCCTACCATATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGa                                                                                                                     <  1:388082/250‑1 (MQ=255)
                                                                                                                            gttgAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAAGGGCGCTAATTTCTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGaa                                                                                                                    <  1:51508/250‑1 (MQ=255)
                                                                                                                                agaAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAAGGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCa                                                                                                                                                           >  2:5129/1‑207 (MQ=255)
                                                                                                                                agaAAACCGGTGAGACGGTCGCTCAAACCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAAGGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCa                                                                                                                                                           <  1:5129/207‑1 (MQ=255)
                                                                                                                                      ccGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAAGGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGCTAAGACAGGGAAATTcc                                                                                                         >  2:295833/1‑251 (MQ=255)
                                                                                                                                                ggTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAAGGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCg                                                                                               <  2:281196/251‑1 (MQ=255)
                                                                                                                                                 gTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAAGGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCgg                                                                                              >  1:163416/1‑251 (MQ=255)
                                                                                                                                                           ccACGCCACGGCGCGTGGGTCAGGCATTTACCCCCGCCCTGGCAGCCGCCTTTTCCCTGACAATGGGTGGCCTCTACAGATTCGACAATGGCGCTAATTACTTCGCCAACGGCGATTCTGCTGGCATCTTTGCCTAACAGATAACCACCGCGTGGTACACGTACGCTGGAAACCAGACCAGTTTTACGCAGACGGGAAAACAGTTGTGCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATGTCAGc                                                                                    <  2:172826/250‑1 (MQ=255)
                                                                                                                                                                     gcgcGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTAAAGGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTAcc                                                                          <  2:7884/251‑1 (MQ=255)
                                                                                                                                                                                               cccTGGCAGCCGCCTTTACCCTGACACCGGGTGGCATCTACAGATTCGTCAAGGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTGGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAAc                                                <  1:414463/251‑1 (MQ=255)
                                                                                                                                                                                                           ccTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAAGGGCGCTACTTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGTCTAACAGATAACCACCGCCTGGTCCACGAACGCTGGAAACCAGACAATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAACTTCCCTGACGTTCGGAAATATCAGCCAAAGGTACCGGGCACGCTTCAGAGTTGAGCGCAACGTCAAGCATTGc                                    >  1:192556/1‑251 (MQ=255)
                                                                                                                                                                                                             tttACCCTGACAACGGGTGGCATCTACAGATTCGTCAAGGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCgg                                  >  1:20016/1‑251 (MQ=255)
                                                                                                                                                                                                                     gACAAGGGGTGGCATCTACAGAGTCGTCAGGGGCGCTATTTACTACGGCTAAGGCGATGCTGGTGGCATCTGTGCCTAACAGATAAGCACCGCCTGGTCCACGTACGCTGGAAACCAGGCCAGTTTTACGCAGACGGGAAAACAGGTGTTCCAGAGAAGAAAGTGAAATTTCCTGGCGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTTAGAGTTGTGCGCAACGTCAAGCATGTCGGTCACgg                            <  2:277372/249‑1 (MQ=255)
                                                                                                                                                                                                                       cAACGGGTGGCATCTACAGATCCGTCAAGGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAAACAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATa                        <  2:291585/251‑1 (MQ=255)
                                                                                                                                                                                                                         aCGGGTGGCGTCTACAGATTCGTCAAGGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCTCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACTGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAgc                      <  1:29297/251‑1 (MQ=255)
                                                                                                                                                                                                                         aCGGGTGGCATCTACAGATTCGGCAAGGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAgc                      <  1:366108/251‑1 (MQ=255)
                                                                                                                                                                                                                                ccATCTCCAGACTCGTCAAGGCCGCAAATTGCTTCGCCAACGGCGGTGCTGCTGGCATCTTTGCCGAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTTACGTTCGGAAACATCAGCCATCGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAATCATTGCGGTCACGGCATAGCGCCCttt               <  1:278440/250‑1 (MQ=255)
                                                                                                                                                                                                                                           ttCGTCAAGGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGtc    >  1:189463/1‑251 (MQ=255)
                                                                                                                                                                                                                                             cGTCAAGGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGtctc  >  1:410574/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
ATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTC  >  NC_000913/2661642‑2662129

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑