Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 362,884 A→G L99S (TTA→TCA)  lacY ← lactose permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913362,8840AG100.0% 488.4 / NA 139L99S (TTA→TCA) lacYlactose permease
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (78/61);  total (78/61)

Only 100 of 148 total aligned reads displayed.

TTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTGTCAG  >  NC_000913/362787‑362980
                                                                                                 |                                                                                                
tttCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGAc                                                                                                <  1:908791/100‑1 (MQ=255)
 ttCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACa                                                                                               >  2:3035991/1‑100 (MQ=255)
  tctcAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAg                                                                                              >  2:811520/1‑100 (MQ=255)
  tctcAATAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCCACTAAAATGTTGTATTGTGACAGt                                                                                             >  1:385154/1‑100 (MQ=255)
   ctcAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAATAAACGATTCTCCTAAAATGTTGTATTGTGACa                                                                                               >  1:1466048/1‑98 (MQ=255)
   ctcAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCGATATAAATACCACCAACAATCGATCCTACTAAAATTTTGTATTGTGACAGt                                                                                             >  2:1028321/1‑100 (MQ=255)
    tcAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTg                                                                                            >  1:3057040/1‑100 (MQ=255)
     cAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTgg                                                                                           <  1:953579/100‑1 (MQ=255)
       aTAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCGCCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGcc                                                                                         >  1:2222784/1‑100 (MQ=255)
       aTAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGcc                                                                                         <  1:191929/100‑1 (MQ=255)
       aTAAATGCCTCTACTGAAGGCGCAACAGCGAAAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGcc                                                                                         <  1:11568/100‑1 (MQ=255)
        tAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGccc                                                                                        >  2:2960718/1‑100 (MQ=255)
         aaaTGCCCCTACTGCTGGCGAACCGGAGTTAAAACAAAAGCCTACATCCATATCACTCACATTCGATCCTACTAAAATGTTGTATTGTGACAGTGTCCCg                                                                                       <  1:2896742/100‑1 (MQ=255)
           aTGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGaa                                                                                     <  1:3441613/100‑1 (MQ=255)
            tGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAg                                                                                    <  2:2226684/100‑1 (MQ=255)
               ctctACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATa                                                                                 >  2:1308059/1‑100 (MQ=255)
               ctctACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATa                                                                                 <  1:1746686/100‑1 (MQ=255)
                   aCTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATaaaaa                                                                             >  1:2536432/1‑100 (MQ=255)
                   aCTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATa                                                                                 <  1:533280/96‑1 (MQ=255)
                    ctgctgGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAAt                                                                            >  2:2412546/1‑100 (MQ=255)
                      gctgGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATaa                                                                          >  1:3252828/1‑100 (MQ=255)
                      gctgGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATaa                                                                          >  2:2325095/1‑100 (MQ=255)
                      gctgGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATaa                                                                          >  1:741747/1‑100 (MQ=255)
                       ctgGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATaaa                                                                         >  2:1413709/1‑100 (MQ=255)
                       ctgGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATaaa                                                                         >  2:1502818/1‑100 (MQ=255)
                        tgGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAg                                                                        >  2:932210/1‑100 (MQ=255)
                        tgGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAg                                                                        <  2:2131556/100‑1 (MQ=255)
                        tgGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAg                                                                        <  1:3159621/100‑1 (MQ=255)
                         ggCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGa                                                                       >  2:1504032/1‑100 (MQ=255)
                         ggCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGa                                                                       >  1:3411661/1‑100 (MQ=255)
                               ccGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGgcg                                                                 >  2:2672815/1‑100 (MQ=255)
                                 ggCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCa                                                               <  1:25568/100‑1 (MQ=255)
                                 ggCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCa                                                               <  2:629579/100‑1 (MQ=255)
                                   cGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCaaa                                                             >  2:1688353/1‑100 (MQ=255)
                                   cGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCaaa                                                             <  2:2211593/100‑1 (MQ=255)
                                     ttAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACa                                                           <  2:939689/100‑1 (MQ=255)
                                      tAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCGCGAAGATAAAAATAAAGAACGGCGCAAACAt                                                          >  2:3210638/1‑100 (MQ=255)
                                      tAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCaaa                                                             >  1:1333102/1‑97 (MQ=255)
                                      tAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACAt                                                          >  2:2278004/1‑100 (MQ=255)
                                       aaaaCAACAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAATGGCCCGAAGCTAAAAATAAAGAACGGCGCAAACAt                                                          >  1:1081344/1‑99 (MQ=255)
                                         aaCAAAAGCCTATATAAATACCACCAACAAACGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATc                                                         <  1:3210586/98‑1 (MQ=255)
                                         aaCAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCAc                                                       >  2:2705559/1‑100 (MQ=255)
                                          aCAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAACTGTTGTATTGTGACAGTGGCCCGAAGACAAAAATAAACAACGGCGCAAACATCACt                                                      >  2:3412212/1‑100 (MQ=255)
                                          aCAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACt                                                      <  1:1567897/100‑1 (MQ=255)
                                          aCAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACt                                                      >  2:3414331/1‑100 (MQ=255)
                                          aCAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGAGGACCGAAGATAAAAATAAAGAACGGCGAAAACATCaat                                                      >  1:1883105/1‑98 (MQ=255)
                                           cAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTa                                                     >  2:2668616/1‑100 (MQ=255)
                                              aaGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACa                                                  <  2:1350613/100‑1 (MQ=255)
                                              aaGCCTAGATAAATACCACCAACAATCGATCCAACTAAAATGTTGTATTGTGACAGTGGCCCGAAAATAAAAATAAAGAACGGCGCAAACATCACTAACa                                                  <  1:2193469/100‑1 (MQ=255)
                                               aGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACAt                                                 >  2:1125185/1‑100 (MQ=255)
                                               aGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACAt                                                 >  2:563284/1‑100 (MQ=255)
                                               aGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACAt                                                 >  2:2866259/1‑100 (MQ=255)
                                                gCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATg                                                >  2:647338/1‑100 (MQ=255)
                                                 ccTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCagatcggaag                                               >  1:1382757/1‑91 (MQ=255)
                                                  tctgATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGcc                                              <  2:1653584/97‑1 (MQ=255)
                                                  cTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGcc                                              >  1:3108828/1‑100 (MQ=255)
                                                    aGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCgg                                            <  2:2536466/100‑1 (MQ=255)
                                                    aGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCgg                                            <  2:1956605/100‑1 (MQ=255)
                                                     gATAAATACCACCAACAATCGATCCTATTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGccaga                                           >  2:2297288/1‑97 (MQ=255)
                                                     gATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGccaga                                           >  1:1653561/1‑97 (MQ=255)
                                                     gATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGt                                           >  1:63976/1‑100 (MQ=255)
                                                      aTAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGta                                          <  1:1522908/100‑1 (MQ=255)
                                                       tAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGtaa                                         >  2:3020855/1‑100 (MQ=255)
                                                       tAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGtaa                                         >  1:220564/1‑100 (MQ=255)
                                                        naaTACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGtaa                                         >  1:1847152/2‑99 (MQ=255)
                                                        aaaTACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGtaat                                        <  2:564504/100‑1 (MQ=255)
                                                        aaaTACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAAAGGCGAAAACATCACTAACATGCCGGaaa                                         >  1:3252200/1‑99 (MQ=255)
                                                           tACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTaataat                                     <  2:1394311/100‑1 (MQ=255)
                                                           tACCACCAACAATCGATCCTACTAAAATGGTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTaataat                                     <  1:2135723/100‑1 (MQ=255)
                                                            accaccAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCTGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATc                                    <  1:1835056/100‑1 (MQ=255)
                                                             ccaccaACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATcc                                   >  1:2597083/1‑100 (MQ=255)
                                                             ccaccaACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATcc                                   <  2:1201061/100‑1 (MQ=255)
                                                              caccaACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATcc                                   <  1:869408/99‑1 (MQ=255)
                                                               accaACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCcac                                 >  2:2300953/1‑100 (MQ=255)
                                                               accaACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCcac                                 <  1:225068/100‑1 (MQ=255)
                                                               accaACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCcac                                 <  2:2901923/100‑1 (MQ=255)
                                                                 caacaaTCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCAcag                               >  1:1645573/1‑100 (MQ=255)
                                                                  aacaaTCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCAcagc                              <  1:1095237/100‑1 (MQ=255)
                                                                  aacaaTCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCAcagc                              >  1:1588125/1‑100 (MQ=255)
                                                                  aacaaTCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCAcagc                              <  2:2811963/100‑1 (MQ=255)
                                                                    caaTCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCAcagcag                            <  2:1964753/100‑1 (MQ=255)
                                                                     aaTCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAgg                           <  1:1800577/100‑1 (MQ=255)
                                                                      aTCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAgg                           <  1:2410943/99‑1 (MQ=255)
                                                                         gATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTAtt                       >  1:341214/1‑100 (MQ=255)
                                                                            nctactaAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTgc                    >  2:1640795/2‑100 (MQ=255)
                                                                            ccTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTgc                    >  2:836157/1‑100 (MQ=255)
                                                                              tactaAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTgcgc                  <  1:2694835/100‑1 (MQ=255)
                                                                               actaAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCa                 <  1:3190750/100‑1 (MQ=255)
                                                                                ctaAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAg                >  1:2045238/1‑100 (MQ=255)
                                                                                  aaaaTGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGcc              <  1:540216/100‑1 (MQ=255)
                                                                                  aaaaTGTTGTATTGTGACAGTGACCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGcc              <  1:3020806/100‑1 (MQ=255)
                                                                                    aaTGTTGTATTGTGACAGTGGCCCGAAGATAAAAAGAAAGAACGGCGCAAACATCACTAGCCTGCCGGTAATAATCCACAGCAGGTATTTGCGCAGccc             <  1:2005251/99‑1 (MQ=255)
                                                                                      tgttgtATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAg          <  1:2427771/100‑1 (MQ=255)
                                                                                       gttgtATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAg          >  2:3094393/1‑99 (MQ=255)
                                                                                       gttgtATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGt         >  2:1641916/1‑100 (MQ=255)
                                                                                            attgtgACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTGTc    >  2:120357/1‑100 (MQ=255)
                                                                                            attgtgACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAAAATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTGTc    >  1:2038958/1‑100 (MQ=255)
                                                                                             ttgtgACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTGTCa   <  2:2689651/100‑1 (MQ=255)
                                                                                             ttgtgACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTGTCa   <  1:3204929/100‑1 (MQ=255)
                                                                                              tgtgACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTGTCAg  <  1:2971072/100‑1 (MQ=255)
                                                                                                 |                                                                                                
TTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTGTCAG  >  NC_000913/362787‑362980

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