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,9080GT95.7% 102.0 / ‑5.4 46A75D (GCC→GAC) iscRDNA‑binding transcriptional dual regulator IscR
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (16/28);  total (17/29)
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.31e-01

GTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTC  >  NC_000913/2661661‑2662152
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
gTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCg                                                                                                                                                                                                                                                    >  2:21612/1‑250 (MQ=255)
 gcGATCGCGTTTTGGGTGCGGGTGCGTGGTGCGTCGTGAGTATGCGGACGTCCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTAGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTTCCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGAt                                                                                                                                                                                                                                                  <  2:374825/250‑1 (MQ=255)
   gATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGCTTCGTCAACGGCTCTACTTACTTCGCCAACGTCGAtgc                                                                                                                                                                                                                                                >  2:154983/1‑251 (MQ=255)
      cgcgTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGAtgctg                                                                                                                                                                                                                                              <  2:341202/250‑1 (MQ=255)
         gTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTgg                                                                                                                                                                                                                                          >  1:95139/1‑251 (MQ=255)
             tgtgtgCGGGTGCGTGGCGCGGCGTGAGTATGCTGACGACCAGGCACATCCAGCACGTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTTAGAAACCCTGTGAGACGGTCGCTCAAATCACGTCACAGCGCGTGGGTCAGGCATTCATCGCCGCCCTGGCATCCGCCGTTACCCTTACAACGGGTGGCATCTACAGAGTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCAt                                                                                                                                                                                                                                       <  1:466607/250‑1 (MQ=255)
                gtgCGGGTGCGTGGCTCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTTGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGTTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGTCAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCtt                                                                                                                                                                                                                                    <  1:318505/250‑1 (MQ=255)
                     ggTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCt                                                                                                                                                                                                                               >  1:459574/1‑250 (MQ=255)
                     ggTGCGTGGCGCGTCGTGAGGATGCTGCCGACCAGACACATTCAGCACTTCCTGGTCAGTAACCAGTTCGCCCAAAGTAAGCGTGGTGAGAAACCCGGTGAGACGGTCGCTCAAGTCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGTCGCCTTTGCCCGTACAACGAGGGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCTACAACGTCGATGCTGCTGGCATCTTTGCCtc                                                                                                                                                                                                                              <  2:69014/251‑2 (MQ=255)
                                            gCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCGCAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTcc                                                                                                                                                                                                        <  2:116878/250‑1 (MQ=255)
                                                acgacCAGACCCATCCAGCACTTCCTGGTTATGAACCAGTGCGCCTAAAGTAATGTGGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGAGTCGTCAACGGCGCTAATTACTTCGCCAACTTCGAGGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTa                                                                                                                                                                                                   <  1:316751/251‑1 (MQ=255)
                                                 cgcccAGACACATCCAGCACGTCCGGGTTATTAGCCAGTGCGCCTAAAGTAATGTTGGTGAGAAAACCGGTGAGACGGTCTCTCAAATCACGCCACAGCGCGTGGTTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCGGATAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCGTGCATCTTTGCTTAACAGATAACCACCGCCTGGTCCACGTAc                                                                                                                                                                                                  <  1:453329/248‑1 (MQ=255)
                                                   acCAGACTCATCAATCACTTCAGGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACTGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATACCCACCGCCTGGTCCACGTACGc                                                                                                                                                                                                <  2:160766/251‑1 (MQ=255)
                                                        acacATCCAGCACTTACTGGTTATTAAACAGTTCGCCTACACTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCCAATCACGCCACAGCGCGTGGGTCCGGCCTTTATCGCCGCCCTGGCCGCCGCATTTCCCCTGACAACGGGTGGCATCTCCAGCTTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAAACACCGCCTGTTCCACGTACGCTGGa                                                                                                                                                                                            >  2:164060/1‑250 (MQ=255)
                                                          acaTCCAGCACTTCCTGGTGATTAACCAGTTCCCCTAAATTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCGCGCCACAGCGCGTGGGTCAGGCAGTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGGGGCATCTACAGATGCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAAc                                                                                                                                                                                         <  2:420792/251‑1 (MQ=255)
                                                                aGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGAcc                                                                                                                                                                                   <  1:435783/251‑1 (MQ=255)
                                                                 ccACTTCCTGGTTATTACCCAGTTCGCCTAAAGGAATGTTGTGGAGAAAACCGGTGAGAGGGTCGCTCAAAGCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTCCCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCa                                                                                                                                                                                  <  1:455580/250‑1 (MQ=255)
                                                                      tCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAACCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCAttt                                                                                                                                                                               <  2:459574/249‑1 (MQ=255)
                                                                       ccTGGTTATTACCCAGTTCGCCTAAAGTAATGTTGTTGAGAAACCCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTa                                                                                                                                                                            <  1:253767/251‑1 (MQ=255)
                                                                               ttaACCAGTTCGCCTAAAGTAATGTTGTTGAGCAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACg                                                                                                                                                                    >  2:471562/1‑251 (MQ=255)
                                                                                taACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACgg                                                                                                                                                                   >  1:451433/1‑251 (MQ=255)
                                                                                    cAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGaaa                                                                                                                                                               <  1:476788/251‑1 (MQ=255)
                                                                                     aGTTCTCCTAAAGTAATGTTGTTGAGAAAACCGGGGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCACCGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGaaaa                                                                                                                                                              <  2:315892/251‑1 (MQ=255)
                                                                                         cGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAgt                                                                                                                                                          >  2:182894/1‑251 (MQ=255)
                                                                                                    aggggtgtgagaatcccGGTGAGACGGTCGCGCAAATGCCGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTTCCCTGGCAGCGGGTGGCATCTACTGAGTCGTCAACGGGGCGAATTACTTCGCCGACGTCGTGGCTGCTGGCATCTTGTCCTAACAGATATCCTGCGCCTGGTCCACGTAGGCTGTAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATa                                                                                                                                               <  1:285715/236‑1 (MQ=255)
                                                                                                     gtttgtagAGAAAACCGGTGAGTCGGTCGCTCAAGGCACGCCACAGCGCGTGGTTCAGGCATTGATCGCCGCCCTTGCAGCCGCCTTTCCCCTGACAACGGGTGGCATCGACAGATTCGTCAACGGCGCTAATGACTGCGCCACCGGCGATGCTTCTGGCATCTTTGCCTAACAGATTACCACCGGCTGTGCCACGTACGCTTGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTGCCAGATac                                                                                                                                              <  2:422615/249‑2 (MQ=255)
                                                                                                          ttgAGAAAACCGGAGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGACAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGaa                                                                                                                                           >  1:401280/1‑249 (MQ=255)
                                                                                                           tgAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAgg                                                                                                                                        >  1:466868/1‑251 (MQ=255)
                                                                                                             agaAAACCGGTGAGACGGTCGCTCAAATCACGCCCCAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGAGTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGa                                                                                                                                      <  2:172382/251‑1 (MQ=255)
                                                                                                                aaaCCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAAt                                                                                                                                   >  1:193759/1‑251 (MQ=255)
                                                                                                                       ggaaaCGGTCGCTAAAATCCCGCCACAGCGCTTGGGTCAGGCATTTATCTCCGCCCTGGCAGCCGCCTTTACCCGGACAACGGGGGGCATCGACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCGGCTGGCTTCTTTGCCTAAGAGATAACCGCCGCCTGGTCCACGTACGCTGGAAACCAGGCCATTTTTACGCAGACGTGAAAACAGTTGTTCCAGATAAGAAACGGAAATTCCCTGa                                                                                                                            <  2:492860/247‑1 (MQ=255)
                                                                                                                          gaCGGTCGCTCAAATCCCGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCCGCCGCCTTTACCCTGACAACGGGTGGCATCTCCAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGt                                                                                                                         <  1:83790/251‑1 (MQ=255)
                                                                                                                                             ccACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACg                                                                                                      <  2:401280/251‑1 (MQ=255)
                                                                                                                                                            cAGGCATTTATCGCCGCCCGGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCtt                                                                                       <  1:271444/251‑1 (MQ=255)
                                                                                                                                                                                                         ggTGGCATTTACAGATCTGTCTACGGCGCTAATTACTTTGCCAACGTCGTTGCTGCTGGCATCTTTGCCTATCAGATACCCACCGCCGGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCACCCAACGGTCCCGGGCCCGCTTCAGAGTTGAGCGCAACGTAAAGCATTGCTGTCACGGCATAGCGcc                                          <  2:278669/251‑1 (MQ=255)
                                                                                                                                                                                                         ggTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGcc                                          >  1:163193/1‑251 (MQ=255)
                                                                                                                                                                                                                   aCAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAg                                    <  1:4498/247‑1 (MQ=255)
                                                                                                                                                                                                                   aCAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATg                                 <  1:154983/250‑1 (MQ=255)
                                                                                                                                                                                                                   aCAGATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGt                                >  1:327523/1‑251 (MQ=255)
                                                                                                                                                                                                                      gATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTc                               >  1:438939/1‑249 (MQ=255)
                                                                                                                                                                                                                      gATTCGTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAg                             <  1:158880/251‑1 (MQ=255)
                                                                                                                                                                                                                           gTCAACGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTCACAGATAACCACCGCCTGGTCCACGTACGCTGGAACCCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCACTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGtctcc                        >  1:304897/1‑250 (MQ=255)
                                                                                                                                                                                                                               aCGGCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTc                    <  1:12085/251‑1 (MQ=255)
                                                                                                                                                                                                                                 ggCGCTAATTACTTCGCCAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCtt                  >  1:67227/1‑251 (MQ=255)
                                                                                                                                                                                                                                      tAATTACTTCGCCAAAGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAACGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGAGGTCACGGCATACCGACCTTTAGATGTCAGTCTCATGTCTTACTTc             >  2:288661/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 ccAACGTCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGGCAGTCTCATGTCTTACTTCACCTCAAACTc  <  2:82862/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
GTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTC  >  NC_000913/2661661‑2662152

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