Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,661,908 G→T A75D (GCC→GAC)  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,9080GT98.3% 178.0 / NA 61A75D (GCC→GAC) iscRDNA‑binding transcriptional dual regulator IscR
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (29/30);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTC  >  NC_000913/2661664‑2662152
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
gATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGAtgc                                                                                                                                                                                                                                                >  1:31915/1‑251 (MQ=255)
   cgcgTCTTGTGTGCGGGGGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGGCGGTCGCTCTAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGAtgctg                                                                                                                                                                                                                                              <  2:164272/250‑1 (MQ=255)
      gTCTTGGGTGCGTGTTCGTGGCGCGTCTTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGTGTGGCATCTTCAGATTCGTCAACGGCTCTAATTACTTCGCCAACGTCGATGCTTCTgg                                                                                                                                                                                                                                          <  2:445606/251‑1 (MQ=255)
          tgtgggCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGc                                                                                                                                                                                                                                         <  2:177819/248‑1 (MQ=255)
           gtgtgAGGGTGCGGGGCGCGTCGTGAGTATGCTGACGTCCAGAAACTTCCAGCTCTTCCTGGTTTTTAACCAGTTCGCTAAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCGCAAATCACGCCACAGCGCGGGGGTCAGGCATTTATCGCCGCCCTGGTAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAAGTACTTCGCCAACGTCGATGCTGCTAGCagt                                                                                                                                                                                                                                      <  2:378277/250‑3 (MQ=255)
                cGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGcc                                                                                                                                                                                                                                <  1:276679/251‑1 (MQ=255)
                     gCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACa                                                                                                                                                                                                                           <  2:210785/251‑1 (MQ=255)
                         ggCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATa                                                                                                                                                                                                                       <  1:458827/251‑1 (MQ=255)
                            gcgTCGTTAGTATGCTGACGACAAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTGAACGTCGCTAATTACGTCGCCAACGTCGATGCTTCGGGCGTTTTTGCGTAACGGATAacc                                                                                                                                                                                                                    <  2:49407/251‑1 (MQ=255)
                              gtcgtTAGTATGCTGAAGACCAGACACCTACAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGCAAACCGGTGAGACGGTCGCTCAAGTCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCCGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCCACGCCGATGCTGCTGACATCTTTGCCTAACAGATcaccac                                                                                                                                                                                                                  >  2:334691/1‑251 (MQ=255)
                                cgtGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAaccacc                                                                                                                                                                                                                 <  2:21011/250‑1 (MQ=255)
                                   gAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCt                                                                                                                                                                                                             <  1:230200/251‑1 (MQ=255)
                                    aGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTg                                                                                                                                                                                                            <  1:501157/251‑1 (MQ=255)
                                     gTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTgg                                                                                                                                                                                                           >  1:191071/1‑251 (MQ=255)
                                               gacCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCCCGTACg                                                                                                                                                                                                 >  1:183081/1‑251 (MQ=255)
                                                acCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGc                                                                                                                                                                                                >  1:34540/1‑251 (MQ=255)
                                                         aTCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGCCGCTCCAATCACGCCACAGCGCGTGCGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATATTTGCCTAACAGATAACCACCGCCTGGTCCAAGTACGCTGGAAAcc                                                                                                                                                                                        >  1:218762/1‑250 (MQ=255)
                                                           ccAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACACCGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAAcc                                                                                                                                                                                        >  2:382069/1‑248 (MQ=255)
                                                                aCTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGAccctt                                                                                                                                                                                >  1:437543/1‑248 (MQ=255)
                                                                              aaCCAGTTCGCCTAAATTAAGGTTTTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTTACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCCCCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACgg                                                                                                                                                                   <  2:416720/250‑1 (MQ=255)
                                                                                ccAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGa                                                                                                                                                                 >  1:285314/1‑250 (MQ=255)
                                                                                    ttCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACg                                                                                                                                                                                                     >  2:269638/1‑210 (MQ=255)
                                                                                    ttCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACg                                                                                                                                                                                                     <  1:269638/210‑1 (MQ=255)
                                                                                     ccGCCTAAAGTAAGGTGGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCCCAGCGCGGGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTTCCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAg                                                                                                                                                           <  2:217408/250‑1 (MQ=255)
                                                                                      cGCCAAAAGTAATGTTGTTGAGAAAACCGGGGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAg                                                                                                                                                           <  2:54698/250‑1 (MQ=255)
                                                                                         cTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGTTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAgttgt                                                                                                                                                       >  2:355643/1‑251 (MQ=255)
                                                                                           aaaGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTc                                                                                                                                                     <  1:497684/251‑1 (MQ=255)
                                                                                               tAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTcc                                                                                                                                                    >  1:192980/1‑248 (MQ=255)
                                                                                                aaTGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTcc                                                                                                                                                    >  1:241561/1‑247 (MQ=255)
                                                                                                         gagaACACCGGTGGGACGGTCACTCAAATAACGCAACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCAGACTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACCGATAAACACCGACTGGTCCACGTACGCTGGCAACCAGACCATTTTTACACAGACGGGCAAAAAGTTGTGCCAGATAAGAAAggg                                                                                                                                       >  2:488257/1‑251 (MQ=255)
                                                                                                          cgaAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGCAAACAGTTGTTCCAGATAAGAAAggg                                                                                                                                       >  2:126221/2‑250 (MQ=255)
                                                                                                          agaAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAggg                                                                                                                                       >  2:2321/1‑250 (MQ=255)
                                                                                                           gaAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGaa                                                                                                                                     <  1:38154/251‑1 (MQ=255)
                                                                                                            aaaaCCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCTACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGaa                                                                                                                                     >  1:26238/1‑250 (MQ=255)
                                                                                                                  ggTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCTCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACCGTTGTTCCAGATAGGAACGGGAAATTCCCt                                                                                                                              >  2:22284/1‑251 (MQ=255)
                                                                                                                  ggTGAGACGGTCGCTCAAATCACGACACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGCAAGGGAAATTCCCt                                                                                                                              >  2:259557/1‑251 (MQ=255)
                                                                                                                     gcgaCGGCCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCGGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAGCAGTTGTTCCAGATAAGAAAGGGAAAGTCCCTGa                                                                                                                            <  1:136119/248‑1 (MQ=255)
                                                                                                                           gTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGAAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGGCCACGTACGCTGGGAACCAGACCATTTTTACGCAGACGGGAAAACCGTTGTTCCCGAAAAGAAAGGATCATTCCCTGACGATCgg                                                                                                                     >  2:10385/1‑251 (MQ=255)
                                                                                                                                  aaaTCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTCCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATc                                                                                                              >  1:300345/1‑251 (MQ=255)
                                                                                                                                          ccACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACg                                                                                                      >  1:370852/1‑251 (MQ=255)
                                                                                                                                              aGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTAc                                                                                                  <  1:15890/251‑1 (MQ=255)
                                                                                                                                                  cgtgtgtcagccgtgtATCGCCGCCCTGTCAGCAGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACTGCGCTAATTCCTTCGCCAACGTCGCTGCTGCTGGTATCTTTGCCTAAAAGATAACCACCGCCTGGTCCACGTACGCTGGAAACTAGACCATCTGTACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCggg                                                                                              <  2:45551/236‑1 (MQ=255)
                                                                                                                                                  cgTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCggg                                                                                              >  2:403761/1‑251 (MQ=255)
                                                                                                                                                                    cgccgcCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAgc                                                                            <  2:251405/251‑1 (MQ=255)
                                                                                                                                                                            gaatgtcagacTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTc                                                                    >  2:14958/11‑251 (MQ=255)
                                                                                                                                                                                              tgaaaaacggTGGCATCTACAAATTCGTACACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGAATCTTTGCCTAACCGATAACCCCCGCCTGTTCCACGTACACTGGAAACCAGAACATTTTTCCGCAGGCGGGACAACAGTTGTTCCAGATAAGAAAGGCAAATTCCCTGACGTTCGGAAATATCAGCCAAAGGTACCGGGACCGCTTCCAAGTTGAGCGCAACGTCACGCCTTGCGGTCACGGc                                                  >  2:248918/9‑251 (MQ=255)
                                                                                                                                                                                               gACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGc                                                  >  1:403522/1‑250 (MQ=255)
                                                                                                                                                                                                    cGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAgcg                                            <  2:370852/251‑1 (MQ=255)
                                                                                                                                                                                                          gCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGCTGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGCGTTGAACGCAACGTCAAGCATTGCGGTCACGGCATAGAGCCCttt                                      >  2:307694/1‑251 (MQ=255)
                                                                                                                                                                                                           cATCTAAAGATTCGGCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGGCAAGCATTGCGGTCACGGCATAGCGCCCTTTa                                     <  1:126221/251‑1 (MQ=255)
                                                                                                                                                                                                                   aaTTCTTCAACGGCGCTAATTACTGCGCCAAAGTCGAGGCTTCTGGCATCTTTGACTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGGTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAg                             <  1:182459/250‑1 (MQ=255)
                                                                                                                                                                                                                             cGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCt                   <  1:470494/251‑1 (MQ=255)
                                                                                                                                                                                                                                 gcTAATTGCTCCGCCAACGTCTATTCTGCTGACATCTGAGCCTTACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATGTTTCCGCAGACGTGAAAACAGTTGTTCCAGACAAGAAAGGGAAAAGCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGATTTGAGTGGAACGTTAAGAATGGCGGTCGCGGCATAGCGTCCTGTATATGTCAGTCGCATATCtatct               <  2:29448/251‑5 (MQ=255)
                                                                                                                                                                                                                                    aaTTCCTTCGCCAACGTCGAGGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCCGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACtt              <  1:423021/249‑1 (MQ=255)
                                                                                                                                                                                                                                     aTTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATCACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGACATTCCCTGACGTTCGGAAATATCAGCCAACGGTACAGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCAc           >  2:409931/1‑251 (MQ=255)
                                                                                                                                                                                                                                        actcgcccAAGGGCGATGCTGCTGGCAGCTTTGCCTAACAGATACCCACCGCCTGGTCCACGTACGCTGGAAACCAGACCAGTGTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGTCATAGCGCAATTTAGATGTCAGTCGCCTGTCTTACTTCTCCt         <  2:224530/244‑1 (MQ=255)
                                                                                                                                                                                                                                        aCTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACGGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATg                      >  1:445238/1‑236 (MQ=255)
                                                                                                                                                                                                                                        aCTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACGGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATg                      <  2:445238/236‑1 (MQ=255)
                                                                                                                                                                                                                                           tCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCaa      >  2:140983/1‑250 (MQ=255)
                                                                                                                                                                                                                                              ccAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTc  <  1:10385/251‑1 (MQ=255)
                                                                                                                                                                                                                                              ccAACGTCGAGGCTGCTGGCATCTTTGCCTAACAGATAATCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCTGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTATTACTTCACCTCAAACTc  <  1:88951/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
GATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTC  >  NC_000913/2661664‑2662152

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