Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,654,330 T→C D271G (GAC→GGC)  rspA ← mandelate racemase/muconate lactonizing enzyme family protein RspA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,654,3300TC100.0% 274.2 / NA 77D271G (GAC→GGC) rspAmandelate racemase/muconate lactonizing enzyme family protein RspA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (43/34);  total (43/34)

CGTGTGAGCCAGTACGTACCTGATACAGCGAAGCAAAATCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGTCCCAGATGCTGTTGAAGACTTCACCCACTGCGATGGGTGTGACGGTATGTTGGCGAATGAGACGGAAGCATTCCTGGTTTTCCGCAGGCGTCGGGTCTTCCATCCAGAACATGCGATAATCTTCAATGCTTT  >  NC_000913/1654199‑1654462
                                                                                                                                   |                                                                                                                                    
cGTGTGAGCCAGTACGTACCTGATACAGCGAAGCAAAATCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATg                                                                                                                              <  2:215219/140‑1 (MQ=255)
cGTGTGAGCCAGTACGTACCTGATACAGCGAAGCAAAATCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATg                                                                                                                              >  1:181254/1‑140 (MQ=255)
  tgtgAGCCAGTACGTACCTGATACAGCGAAGCAAAATCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCt                                                                                                                            >  1:803864/1‑140 (MQ=255)
       gCCAGTACGTACCTGATACAGCGAAGCAAAATCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGa                                                                                                                       >  1:461611/1‑140 (MQ=255)
           gtacgtacCTGATACAGCGAAGCAAAATCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGAc                                                                                                                   <  2:579181/140‑1 (MQ=255)
            tacgtacCTGATACAGCGAAGCAAAATCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACt                                                                                                                  >  1:453538/1‑140 (MQ=255)
            tacgtacCTGATACAGCGAAGCAAAATCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACt                                                                                                                  >  2:305403/1‑140 (MQ=255)
            tacgtacCTGATACAGCGAAGCAAAATCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACt                                                                                                                  >  2:502015/1‑140 (MQ=255)
            tacgtacCTGATACAGCGAAGCAAAATCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACt                                                                                                                  >  1:943102/1‑140 (MQ=255)
            tacgtacCTGATACAGCGAAGCAAAATCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACt                                                                                                                  >  2:978841/1‑140 (MQ=255)
            tacgtacCTGATACAGCGAAGCAAAATCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACt                                                                                                                  <  2:803864/140‑1 (MQ=255)
              cgtacCTGATACAGCGAAGCAAAATCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTc                                                                                                                >  1:846362/1‑140 (MQ=255)
                tacCTGATACAGCGAAGCAAAATCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAg                                                                                                                     >  1:58223/1‑133 (MQ=255)
                tacCTGATACAGCGAAGCAAAATCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAg                                                                                                                     <  2:58223/133‑1 (MQ=255)
                tacCTGATACAGCGAAGCAAAATCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCAc                                                                                                              <  2:461611/140‑1 (MQ=255)
                tacCTGATACAGCGAAGCAAAATCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCAc                                                                                                              <  1:428300/140‑1 (MQ=255)
                 acCTGATACAGCGAAGCAAAATCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCAcc                                                                                                             >  2:760680/1‑140 (MQ=255)
                 acCTGATACAGCGAAGCAAAATCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCAcc                                                                                                             >  2:638246/1‑140 (MQ=255)
                 acCTGATACAGCGAAGCAAAATCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCAcc                                                                                                             <  1:978841/140‑1 (MQ=255)
                      aTACAGCGAAGCAAAATCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGa                                                                                                                       <  1:460190/125‑1 (MQ=255)
                      aTACAGCGAAGCAAAATCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGa                                                                                                                       >  2:460190/1‑125 (MQ=255)
                      aTACAGCGAAGCAAAATCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTg                                                                                                        >  1:246903/1‑140 (MQ=255)
                      aTACAGCGAAGCAAAATCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTg                                                                                                        >  1:666534/1‑140 (MQ=255)
                        aCAGCGAAGCAAAATCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCg                                                                                                      <  2:666534/140‑1 (MQ=255)
                        aCAGCGAAGCAAAATCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCg                                                                                                      <  2:658947/140‑1 (MQ=255)
                        aCAGCGAAGCAAAATCGGCAATACGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCg                                                                                                      >  2:932002/1‑140 (MQ=255)
                         cAGCGAAGCAAAATCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGa                                                                                                     >  2:582942/1‑140 (MQ=255)
                          aGCGAAGCAAAATCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGAt                                                                                                    >  1:62178/1‑140 (MQ=255)
                          aGCGAAGCAAAATCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGAt                                                                                                    >  1:129200/1‑140 (MQ=255)
                            cGAAGCAAAATCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATgg                                                                                                  >  2:129146/1‑140 (MQ=255)
                            cGAAGCAAAATCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATgg                                                                                                  >  1:760290/1‑140 (MQ=255)
                            cGAAGCAAAATCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATgg                                                                                                  >  2:525749/1‑140 (MQ=255)
                             gAAGCAAAATCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGAt                                                                                                    >  1:642904/1‑137 (MQ=255)
                             gAAGCAAAATCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGAt                                                                                                    <  2:642904/137‑1 (MQ=255)
                             gAAGCAAAATCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATggg                                                                                                 >  2:99141/1‑140 (MQ=255)
                              aaGCAAAATCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATGGgt                                                                                                <  1:129146/140‑1 (MQ=255)
                               accAAAATCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATGGgtg                                                                                               <  2:943102/138‑1 (MQ=255)
                                gCAAAATCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATGGgtgt                                                                                              >  1:284869/1‑140 (MQ=255)
                                  aaaaTCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATGGGTGTGa                                                                                            <  2:62178/140‑1 (MQ=255)
                                    aaTCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGAt                                                                                                    <  1:174345/130‑1 (MQ=255)
                                    aaTCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGAt                                                                                                    >  2:174345/1‑130 (MQ=255)
                                      tCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATGGGTGTGACGGt                                                                                        <  1:117846/140‑1 (MQ=255)
                                      tCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATGGGTGTGACGGt                                                                                        <  2:892044/140‑1 (MQ=255)
                                      tCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATGGGTGTGACGGt                                                                                        <  2:846362/140‑1 (MQ=255)
                                      tCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATGGGTGTGACGGt                                                                                        <  1:591521/140‑1 (MQ=255)
                                                   gcATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATGGGTGTGACGGTATGTTGGCGAATg                                                                           <  2:760290/140‑1 (MQ=255)
                                                       aCCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATGGGTGTGACGGTATGTTGGCGAATGAGAc                                                                       >  1:163817/1‑140 (MQ=255)
                                                           gTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATGGGTGTGACGGTATGTTGGCGAATGAGACGGaa                                                                   <  2:453538/140‑1 (MQ=255)
                                                            tAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATGGGTGTGACGGTATGTTGGCGAATGAGACGGAAg                                                                  <  2:163817/140‑1 (MQ=255)
                                                            tAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATGGGTGTGACGGTATGTTGGCGAATGAGACGGAAg                                                                  <  1:638246/140‑1 (MQ=255)
                                                                  cGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATGGGTGTGACGGTATGTTGGCGAATGAGACGGAAGCATTcc                                                            >  2:561262/1‑140 (MQ=255)
                                                                   gCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCATCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATGGGTGTGACGGTATGTTGGCGAATGAGACGGa                                                                    <  2:925850/131‑1 (MQ=255)
                                                                   gCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATGGGTGTGACGGTATGTTGGCGAATGAGACGGa                                                                    >  1:925850/1‑131 (MQ=255)
                                                                     cTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATGGGTGTGACGGTATGTTGGCGAATGAGACGGAAGCATTCCTgg                                                         >  2:222109/1‑140 (MQ=255)
                                                                     cTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATGGGTGTGACGGTATGTTGGCGAATGAGACGGAAGCATTCCTgg                                                         >  2:934327/1‑140 (MQ=255)
                                                                                      gtgCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATGGGTGTGACGGTATGTTGGCGAATGAGACGGAAGCATTCCTGGTTTTCCGCAGGCGTCgg                                        <  1:99141/140‑1 (MQ=255)
                                                                                       tgCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATGGGTGTGACGGTATGTTGGCGAATGAGACGGAAGCATTCCTGGTTTTCCGCAg                                               >  1:408547/1‑132 (MQ=255)
                                                                                       tgCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATGGGTGTGACGGTATGTTGGCGAATGAGACGGAAGCATTCCTGGTTTTCCGCAg                                               <  2:408547/132‑1 (MQ=255)
                                                                                       tgCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATGGGTGTGACGGTATGTTGGCGAATGAGACGGAAGCATTCCTGGTTTTCCGCAGGCGTCggg                                       >  1:459722/1‑140 (MQ=255)
                                                                                        gCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATGGGTGTGACGGTATGTTGGCGAATGAGACGGAAGCATTCCTGGTTTTCCGCAGGCGTCGGGt                                      <  1:305403/140‑1 (MQ=255)
                                                                                        gCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATGGGTGTGACGGTATGTTGGCGAATGAGACGGAAGCATTCCTGGTTTTCCGCAGGCGTCGGGt                                      <  1:753587/140‑1 (MQ=255)
                                                                                            atatAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATGGGTGTGACGGTATGTTGGCGAATGAGACGGAAGCATTCCTGGTTTTCCGCAGGCGTCGGGTCTtc                                  >  2:443154/1‑140 (MQ=255)
                                                                                                 aTCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATGGGTGTGACGGTATGTAGGCGAATGAGACGGAAGCATTCCTGGTTTTCCGCAGGCGTCGGGTCTtccatcc                             >  2:192490/1‑140 (MQ=255)
                                                                                                   cGATGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATGGGTGTGACGGTATGTTGGCGAATGAGACGGAAGCATTCCTGGTTTTCCGCAGGCGTCGGGTCTTCCATCCAg                           <  2:112712/140‑1 (MQ=255)
                                                                                                      tGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATGGGTGTGACGGTATGTTGGCGAATGAGACGGAAGCATTCCTGGTTTTCCGCAGGCGTCGGGTCtt                                   >  1:499215/1‑129 (MQ=255)
                                                                                                      tGAGTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATGGGTGTGACGGTATGTTGGCGAATGAGACGGAAGCATTCCTGGTTTTCCGCAGGCGTCGGGTCtt                                   <  2:499215/129‑1 (MQ=255)
                                                                                                         gTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCTCTGCGATGGGTGTGACGGTATGTTGGCGAATGAGACGGAAGCATTCCTGGTTTTCCGCAGGCGTCGGGTCTTCCATCCAGAACATg                     <  1:881607/140‑1 (MQ=255)
                                                                                                         gTTGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATGGGTGTGACGGTATGTAGGCGAATGAGACGGAAGCATTCCTGGTTTTCCGCAGGCGTCGGGTCTTCCATCCAGAACATg                     <  1:192490/140‑1 (MQ=255)
                                                                                                          ttGCTCTTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATGGGTGTGACGGTATGTTGGCGAATGAGACGGAAGCATTCCTGGTTTTCCGCAGGCGTCGGGTCTTCCATCCAGAACATGc                    >  1:176686/1‑140 (MQ=255)
                                                                                                             ctctTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATGGGTGTGACGGTATGTTGGCGAATGAGACGGAAGCATTCCTGGTTTTCCGCAGGCGTCGGGTCTTCCATCCAGAACATGCGAt                 >  1:344402/1‑140 (MQ=255)
                                                                                                               ctTCAATCAGTTGTTTGCAGCCCCAGATGCTGTTGTAGACTTCACCCACTGCGATGGGTGTGACGGTATGTTGGCGATTGAGACGGAAGCATTCCTGGTTTTCCGCAGGCGTCGGGTCTTCCATCCAGAACATGCGATaa               >  2:83339/1‑140 (MQ=255)
                                                                                                                  cAATCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATGGGTGTGACGGTATGTTGGCGAATGAGACGGAAGCATTCCTGGTTTTCCGCAGGCGTCGGGTCTTCCATCCAGAACATGCGATAATCt            >  2:75109/1‑140 (MQ=255)
                                                                                                                    aTCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATGGGTGTGACGGTATGTTGGCGAATGAGACGGAAGCATTCCTGGTTTTCCGCAGGCGTCGGGTCTTCCATCCAGAACATGCGATAATCTTc          >  1:449236/1‑140 (MQ=255)
                                                                                                                     tCAGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATGGGTGTGACGGTATGTTGGCGAATGAGACGGAAGCATTCCTGGTTTTCCGCAGGCGTCGGGTCTTCCATCCAGAACATGCGATAATCTTCa         <  2:129200/140‑1 (MQ=255)
                                                                                                                       aGTTGTTTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATGGGTGTGACGGTATGTTGGCGAATGAGACGGAAGCATTCCTGGTTTTCCGCAGGCGTCGGGTCTTCCATCCAGAACATGCGATAATCTTCAAt       <  2:176686/140‑1 (MQ=255)
                                                                                                                        gttgttTGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATGGGTGTGACGGTATGTTGGCGAATGAGACGGAAGCATTCCTGGTTTTCCGCAGGCGTCGGGTCTTCCATCCAGAACATGCGATAATCTTCAATg      >  2:879390/1‑140 (MQ=255)
                                                                                                                            tttGCAGCCCCAGATGCTGTTGAAGACTTCACCCACTGCGATGGGTGTGACGGTATGTTGGCGAATGAGACGGAAGCATTCCTGGTTTTCCGCAGGCGTCGGGTCTTCCATCCAGAACATGCGATAATCTTCAATGCttt  >  2:842617/1‑140 (MQ=255)
                                                                                                                                   |                                                                                                                                    
CGTGTGAGCCAGTACGTACCTGATACAGCGAAGCAAAATCGGCAATCCGGCGCATACCGGTAATTCCGCCTGCATGGGTCAGCGTGGTGCGGATATAATCGATGAGTTGCTCTTCAATCAGTTGTTTGCAGTCCCAGATGCTGTTGAAGACTTCACCCACTGCGATGGGTGTGACGGTATGTTGGCGAATGAGACGGAAGCATTCCTGGTTTTCCGCAGGCGTCGGGTCTTCCATCCAGAACATGCGATAATCTTCAATGCTTT  >  NC_000913/1654199‑1654462

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