Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,585,517 G→A A436V (GCC→GTC)  ggt ← glutathione hydrolase proenzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,585,5170GA100.0% 132.5 / NA 40A436V (GCC→GTC) ggtglutathione hydrolase proenzyme
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/25);  total (15/25)

TGATGGTGGAAACGCGGCGCATTGGTCGCTTCGGCGACGTTCAAGCCATAATCGATGCTATTCACCACCATTTGCAGCACTGTAGTGATGATCCGGCTACCGCCTGGGCTACCGGTAACCAGCCAGGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTGCCGGGGCGAATTTCGCTGGATGGCTTCGCTTTATTGATATCAATTTGATCGGCAATAGAT  >  NC_000913/3585274‑3585741
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                
tGATGGTGGAAACGCGGCGCATTGGTCGCTTCGGCGACGTTCAAGCCATAATCGATGCTATTCACCACCATTTGCAGCACTGTAGTGATGATCCGGCTACCGCCTGGGCTACCGGTAACCAGCCAGGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGgaga                                                                                                                                                                                                                            >  1:269608/1‑250 (MQ=255)
           ccgcgGCGCATTGGTCGCTTCGGCGACGTTCAAGCCATAATCGATGCTATTCACCACCATTTGCAGCACTGTAGTGATGATCCGGCTACCGCCTGGGCTACCGGTAACCAGCCAGGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCt                                                                                                                                                                                                                <  2:65301/250‑1 (MQ=255)
            cgcgGCGCATTGTTCGCTGGGGCGTCGTTAAAGCCATAATCGATGCTATTCACCACCATTTGCAGCACTGTAGTGATGATCCGGCTACCGCCTGGGCTACCGGTAACCAGCCAGGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTg                                                                                                                                                                                                               <  2:110185/251‑1 (MQ=255)
                           gCTTCGGCGACGTTCAAGCCATAATCGATGCTATTCACCACCATTTGCAGCACTGTAGTGATGATCCGGCTACCGCCTGGGCTACCGGTAACCAGCCAGGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAAt                                                                                                                                                                                                >  1:87384/1‑251 (MQ=255)
                              tCGGGGACGTTAAAGCCATAATCGATGCTATTCACCACCATTTGCAGCACTGTAGTGATGATCCGGCTACCGCCTGGGCTACCGGTAACCAGCCAGGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATAc                                                                                                                                                                                              <  2:214809/250‑1 (MQ=255)
                                 gCGACTTTCAAGCCATAATCGATGCTATTCACCACCATTTGCAGCACTGTAGTGATGATCCGGCTACCGCCTGGGCTACCGGTAACCAGCCAGGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGct                                                                                                                                                                                          <  1:170767/251‑1 (MQ=255)
                                                  aTCGATGCTATTCACCACCATTTGCAGCACTGTAGTGATGATCCGGCTACCGCCTGGGCTACCGGTAACCAGCCAGGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCAACCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGc                                                                                                                                                                         >  2:296060/1‑251 (MQ=255)
                                                                caaccaTTTGCAGCACTGTAGTGATGATCCGGCTACCGCCTGGGCTACCGGTAACCAGCCAGGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTACACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCCTCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAgggg                                                                                                                                                           >  2:261305/3‑249 (MQ=255)
                                                                   ccagtggcagcaCTGGCGTGATGATCCGGCTACCGCCTGGGCTACCGGTACCCAGCCAGGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACAGCTCTCGCCCGCGACAATGCCCGTACCGAAGgtggtg                                                                                                                                                        <  2:131767/245‑1 (MQ=255)
                                                                            gcaCTGTAGTGATGATCCGGCTACCGCCTGGGTTACCGGTAACCAGCCAGGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAggtggagctaagcgt                                                                                                                                                >  2:166651/1‑240 (MQ=255)
                                                                                   aGTGATGATCCGGCTACCGCCTGGGCTACCGGTAACCAGCCAGGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTc                                                                                                                                         <  1:171848/250‑1 (MQ=255)
                                                                                       atgatCCGGCTACCGCCTGGGCTACCGGTAACCAGCCAGGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTcaccgc                                                                                                                                    >  1:10639/1‑251 (MQ=255)
                                                                                                  aCCGCCTGGGCTACCGGTAACCAGCCAGGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTAc                                                                                                                         <  2:112566/251‑1 (MQ=255)
                                                                                                     gCCTGGGCTACCGGTAACCAGCCAGGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCAt                                                                                                                      >  2:47612/1‑251 (MQ=255)
                                                                                                                cGGTAACCAGCCAGGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCAGCTTTATccacc                                                                                                           <  2:203624/251‑1 (MQ=255)
                                                                                                                       cagccagGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGt                                                                                                    <  2:248735/251‑1 (MQ=255)
                                                                                                                         gccagGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTaa                                                                                                  <  2:53869/251‑1 (MQ=255)
                                                                                                                         gccagGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTa                                                                                                   >  1:53869/1‑250 (MQ=255)
                                                                                                                          ccagGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAAt                                                                                                 >  1:69814/1‑251 (MQ=255)
                                                                                                                                tttACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTa                                                                                            <  1:17364/250‑1 (MQ=255)
                                                                                                                                  tACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGt                                                                                          <  1:12737/250‑1 (MQ=255)
                                                                                                                                            caccacAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGTTTGATTActct                                                                               <  2:10639/251‑1 (MQ=255)
                                                                                                                                                         ggCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATg                                                                                                <  1:235950/221‑1 (MQ=255)
                                                                                                                                                         ggCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGGGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATg                                                                                                >  2:235950/1‑221 (MQ=255)
                                                                                                                                                                           gggCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTGCCGGGGCGAATTTc                                                <  2:74413/251‑1 (MQ=255)
                                                                                                                                                                             gCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTGCCGGGGCGAATTTCGc                                              <  1:47612/251‑1 (MQ=255)
                                                                                                                                                                                      cGGCCCGACTGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATCCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTGACCTCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGGTTGATTACGCTCATAAGGCGCAAGCTTGCCGGGGCGAATTTCGCTGGATGGCt                                     <  1:283958/251‑1 (MQ=255)
                                                                                                                                                                                      cGGCCCGACGGCGTTGGCATCACCGCCCACCAGCACGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCTATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGTTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGATTAATGGGTAGTGTGATTACTCTCATAAGGCGCAAGCTTGCCGGGGCGCAGTTCGCTGGATGGc                                      <  2:271431/250‑1 (MQ=255)
                                                                                                                                                                                             aCGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGTTTGAGTACTCTCATAAGGCGCAAGCTTGCCGGGGCGAATTTCGCTGGATGGCTTCGCttt                              <  1:204986/251‑1 (MQ=255)
                                                                                                                                                                                               agCGTTGGCATTACCGCCCACCAGACCGTAAACGTTCGGTACGCCCGGTTTGACGGAGACATCCTCCATCTGGTTATTAAGCCGAAAACCGATCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGCATAGGTCCCCGCCACCGCGTTACCATATTTATCCACCACTGAGTAATGGGTAGTTTTATTACTCTCTTAAGGCGCAATCTTGCCGGGGCGAATTTCGAGGGATTACTTCGCTTTAt                            >  2:243052/2‑251 (MQ=255)
                                                                                                                                                                                                gCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTGCCGGGGCGAATTTCGCTGGATGGCTTCGCTTTAtt                           >  1:63614/1‑251 (MQ=255)
                                                                                                                                                                                                      gCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTGCCGGGGCGAATTTCGCTGGATGGCTTCGCTTTATTGATATc                     <  1:227135/251‑1 (MQ=255)
                                                                                                                                                                                                       cATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTGCCGGGGCGAATTTCGCTGGATGGCTTCGCTTTATTGATATCa                    <  1:196656/251‑1 (MQ=255)
                                                                                                                                                                                                          cACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTGCCGGGGCGAATTTCGCTGGATGGCTTCGCTTTATTGATATCAAtt                 >  1:132184/1‑251 (MQ=255)
                                                                                                                                                                                                            ccGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCAGCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTGCCGGGGCGAATTTCGCTGGATGGCTTCGCTTTATTGATATCAATTTg               <  2:103191/251‑1 (MQ=255)
                                                                                                                                                                                                                      gCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGTTTGATTCCTCTCATAAGGCGCAAGCTTGCCGGGGCGAATTTCGCTGGATGGCTTCGCTTTATTGATATCAATTTGATCGGCAAt      >  2:119057/1‑250 (MQ=255)
                                                                                                                                                                                                                      gCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTGCCGGGGCGAATTTCGCTGGATGGCTTCGCTTTATTGATATCAATTTGATCGGCAAt      >  1:86369/1‑250 (MQ=255)
                                                                                                                                                                                                                       cccGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTGCCGGGGCGAATTTCGCTGGATGGCTTCGCTTTATTGATATCAATTTGATCGGCAATAg    <  2:253476/251‑1 (MQ=255)
                                                                                                                                                                                                                        ccGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCGTTATCCACCACTGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTGCCGGGGCGAATTTCGCTGGATGGCTTCGCTTTATTGATATCAATTTGATCGGCAATAGa   <  1:13370/251‑1 (MQ=255)
                                                                                                                                                                                                                         cGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTGCCGGGGCGAATTTCGCTGGATGGCTTCGCTTTATTGATATCAATTTGATCGGCAATAGAt  <  2:87384/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                
TGATGGTGGAAACGCGGCGCATTGGTCGCTTCGGCGACGTTCAAGCCATAATCGATGCTATTCACCACCATTTGCAGCACTGTAGTGATGATCCGGCTACCGCCTGGGCTACCGGTAACCAGCCAGGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTGCCGGGGCGAATTTCGCTGGATGGCTTCGCTTTATTGATATCAATTTGATCGGCAATAGAT  >  NC_000913/3585274‑3585741

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