Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,667,370 G→A I159I (ATC→ATT hcaT ← putative 3‑phenylpropionate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,667,3700GA100.0% 714.3 / NA 174I159I (ATC→ATThcaTputative 3‑phenylpropionate transporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (93/81);  total (93/81)

Only 100 of 176 total aligned reads displayed.

TGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGTCGCACTTTGCCGTAATCAA  >  NC_000913/2667273‑2667467
                                                                                                 |                                                                                                 
tGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATc                                                                                                <  1:3450291/101‑1 (MQ=255)
tGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATc                                                                                                <  1:3450435/101‑1 (MQ=255)
 gCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCa                                                                                               >  2:3450830/1‑101 (MQ=255)
  ctctCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCAc                                                                                              >  1:3450834/1‑101 (MQ=255)
   tctcCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCAcc                                                                                             >  1:3450836/1‑101 (MQ=255)
    ctcCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCAcc                                                                                             <  2:3450477/100‑1 (MQ=255)
     tcCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCg                                                                                           >  1:3450837/1‑101 (MQ=255)
     tcCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCg                                                                                           >  2:3450838/1‑101 (MQ=255)
     tcCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCg                                                                                           <  2:3450546/101‑1 (MQ=255)
     tcCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCg                                                                                           <  2:3450448/101‑1 (MQ=255)
       ctgctgGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGAt                                                                                         >  2:3450841/1‑101 (MQ=255)
         gctgGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATaa                                                                                       <  2:3450600/101‑1 (MQ=255)
         gctgGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATaa                                                                                       >  1:3450843/1‑101 (MQ=255)
         gctgGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATaa                                                                                       <  1:3450670/101‑1 (MQ=255)
         gctgGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATaa                                                                                       <  2:3450396/101‑1 (MQ=255)
          ctgGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAAt                                                                                      <  1:3450538/101‑1 (MQ=255)
          ctgGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAAt                                                                                      <  2:3450563/101‑1 (MQ=255)
           tgGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATaa                                                                                       <  1:3450679/99‑1 (MQ=255)
               gACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAAc                                                                                 <  2:3450527/101‑1 (MQ=255)
               gACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAAc                                                                                 <  1:3450569/101‑1 (MQ=255)
               gACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAAc                                                                                 <  1:3450428/101‑1 (MQ=255)
                aCTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACa                                                                                <  2:3450578/101‑1 (MQ=255)
                  ttGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATa                                                                              <  2:3450664/101‑1 (MQ=255)
                     ccccTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTg                                                                           <  1:3450601/101‑1 (MQ=255)
                       ccTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGt                                                                            <  2:3450379/98‑1 (MQ=255)
                           gtgGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGt                                                                     <  1:3450508/101‑1 (MQ=255)
                             ggCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGttt                                                                   <  2:3450555/101‑1 (MQ=255)
                             ggCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGttt                                                                   <  1:3450370/101‑1 (MQ=255)
                             ggCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGttt                                                                   >  1:3450846/1‑101 (MQ=255)
                              gCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTg                                                                  >  2:3450847/1‑101 (MQ=255)
                               cTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGc                                                                 <  1:3450492/101‑1 (MQ=255)
                                tGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCg                                                                                           <  2:3450663/74‑1 (MQ=255)
                                tGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGcc                                                                >  2:3450850/1‑101 (MQ=255)
                                   aTCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGt                                                             >  1:3450851/1‑101 (MQ=255)
                                   aTCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGt                                                             <  2:3450515/101‑1 (MQ=255)
                                    tcgtcgGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTc                                                            >  1:3450852/1‑101 (MQ=255)
                                     cgtcgGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATa                                                                                        <  1:3450461/72‑1 (MQ=255)
                                     cgtcgGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCa                                                           <  2:3450712/101‑1 (MQ=255)
                                     cgtcgGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCa                                                           >  2:3450854/1‑101 (MQ=255)
                                      gtcgGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAg                                                          >  2:3450856/1‑101 (MQ=255)
                                       tcgGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCg                                                              >  2:3450859/1‑96 (MQ=255)
                                       tcgGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAgc                                                         <  2:3450686/101‑1 (MQ=255)
                                       tcgGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAgc                                                         >  1:3450858/1‑101 (MQ=255)
                                          gacggaTGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGcc                                                      <  2:3450523/101‑1 (MQ=255)
                                          gacggaTGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGcc                                                      >  2:3450861/1‑101 (MQ=255)
                                          gacggaTGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGcc                                                      <  1:3450789/101‑1 (MQ=255)
                                            cggaTGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGa                                                    >  2:3450863/1‑101 (MQ=255)
                                            cggaTGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGa                                                    <  1:3450525/101‑1 (MQ=255)
                                             ggaTGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAg                                                   >  1:3450864/1‑101 (MQ=255)
                                               aTGAGAAAGCCGAGCGGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCAGATAATCGAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGcc                                                 >  1:3450870/1‑101 (MQ=255)
                                               aTGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGcc                                                 >  1:3450868/1‑101 (MQ=255)
                                               aTGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGcc                                                 >  2:3450867/1‑101 (MQ=255)
                                                 gagaAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCaa                                               >  1:3450872/1‑101 (MQ=255)
                                                    aaaGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGAcc                                                                        >  2:3450877/1‑73 (MQ=255)
                                                    aaaGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGa                                            >  1:3450874/1‑101 (MQ=255)
                                                    aaaGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGa                                            >  2:3450875/1‑101 (MQ=255)
                                                      aGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACa                                          >  1:3450879/1‑101 (MQ=255)
                                                       gCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACaa                                         >  2:3450880/1‑101 (MQ=255)
                                                        ccGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACaaa                                        <  1:3450643/101‑1 (MQ=255)
                                                        ccGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACaa                                         >  2:3450882/1‑100 (MQ=255)
                                                         cGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAAc                                       >  1:3450883/1‑101 (MQ=255)
                                                         cGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTCGCCCGTCAGCGCCGAGCCAATGACAAAc                                       <  1:3450471/101‑1 (MQ=255)
                                                          gAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACg                                      <  1:3450528/101‑1 (MQ=255)
                                                          gAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACg                                      <  2:3450568/101‑1 (MQ=255)
                                                            gcagcaTGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGcc                                    <  1:3450456/101‑1 (MQ=255)
                                                             cagcaTGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCa                                   >  1:3450886/1‑101 (MQ=255)
                                                              agcaTGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCAc                                  <  1:3450580/101‑1 (MQ=255)
                                                              agcaTGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCAc                                  >  2:3450887/1‑101 (MQ=255)
                                                                caTGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCg                                >  2:3450891/1‑101 (MQ=255)
                                                                caTGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCg                                >  1:3450890/1‑101 (MQ=255)
                                                                 aTGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGa                               >  2:3450892/1‑101 (MQ=255)
                                                                 aTGGATGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGa                               <  1:3450255/101‑1 (MQ=255)
                                                                     aTGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGccc                           >  2:3450894/1‑101 (MQ=255)
                                                                      tGCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGcccc                          >  1:3450896/1‑101 (MQ=255)
                                                                       gCCACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCca                         <  1:3450502/101‑1 (MQ=255)
                                                                        ccACGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCcac                        >  2:3450898/1‑101 (MQ=255)
                                                                           cGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACaaa                                        <  2:3450813/82‑1 (MQ=255)
                                                                           cGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGt                     >  2:3450900/1‑101 (MQ=255)
                                                                           cGCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGt                     >  2:3450901/1‑101 (MQ=255)
                                                                            gCCCAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGTc                    >  1:3450902/1‑101 (MQ=255)
                                                                             cccAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGccc                           <  2:3450464/93‑1 (MQ=255)
                                                                             cccAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGTCg                   >  2:3450904/1‑101 (MQ=255)
                                                                              ccAACGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGTCGc                  >  2:3450906/1‑101 (MQ=255)
                                                                                aaCGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGTCGCAc                >  1:3450910/1‑101 (MQ=255)
                                                                                aaCGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGTCGCAc                >  2:3450909/1‑101 (MQ=255)
                                                                                 aCGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGTCGCACt               <  1:3450622/101‑1 (MQ=255)
                                                                                 aCGTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGTCGCACt               <  1:3450618/101‑1 (MQ=255)
                                                                                   gTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGTCGCACttt             >  1:3450912/1‑101 (MQ=255)
                                                                                   gTCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGTCGCACttt             <  1:3450559/101‑1 (MQ=255)
                                                                                    tCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAg                              <  1:3450488/83‑1 (MQ=255)
                                                                                    tCAACAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAg                              >  2:3450913/1‑83 (MQ=255)
                                                                                      aaCAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGTCGCACTTTGcc          <  1:3450567/101‑1 (MQ=255)
                                                                                       aCAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGTCGCACTTTGCCg         >  2:3450916/1‑101 (MQ=255)
                                                                                        cAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACATTCGCACTTTGCCGt        >  1:3450917/1‑101 (MQ=255)
                                                                                        cAGCGCGAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGTCGCACttt             <  1:3450634/96‑1 (MQ=255)
                                                                                            gcgAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGTCGCACTTTGCCGTAATc    >  2:3450920/1‑101 (MQ=255)
                                                                                             cgAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGTCGCACTTTGCCGTAATCa   >  2:3450923/1‑101 (MQ=255)
                                                                                             cgAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGTCGCACTTTGCCGTAATCa   >  1:3450924/1‑101 (MQ=255)
                                                                                             cgAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGTCGCACTTTGCCGTAATCa   >  1:3450922/1‑101 (MQ=255)
                                                                                             cgAGAATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGTCGCACTTTGCCGTAATCa   <  1:3450655/101‑1 (MQ=255)
                                                                                                 |                                                                                                 
TGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGTCGCACTTTGCCGTAATCAA  >  NC_000913/2667273‑2667467

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