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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913362,8840AG96.5% 274.1 / ‑0.6 86L99S (TTA→TCA) lacYlactose permease
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (49/34);  total (50/36)
Fisher's exact test for biased strand distribution p-value = 5.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

TTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTGTCAG  >  NC_000913/362787‑362980
                                                                                                 |                                                                                                
tttCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGAc                                                                                                <  1:1228031/100‑1 (MQ=255)
 ttCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACa                                                                                               >  1:351594/1‑100 (MQ=255)
 ttCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACa                                                                                               <  2:258537/100‑1 (MQ=255)
  tctcAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAg                                                                                              <  2:2481505/100‑1 (MQ=255)
   ctcAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGt                                                                                             <  2:93454/100‑1 (MQ=255)
   ctcAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGt                                                                                             >  1:1070419/1‑100 (MQ=255)
    tcAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTg                                                                                            >  2:2296651/1‑100 (MQ=255)
     cAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTgg                                                                                           >  2:242450/1‑100 (MQ=255)
      aaTAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGc                                                                                          >  2:302981/1‑100 (MQ=255)
      aaTAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGc                                                                                          >  2:189536/1‑100 (MQ=255)
        tAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAATCCTAGATAAATACAACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTgg                                                                                           >  2:891724/1‑97 (MQ=255)
          aaTGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGa                                                                                      >  2:2204996/1‑100 (MQ=255)
               ctctACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATa                                                                                 >  2:32734/1‑100 (MQ=255)
               ctctACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATa                                                                                 >  1:332158/1‑100 (MQ=255)
                 ctACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTAATAAAATGTTGTATTGTGACAGTGGCCCGAAGATaaa                                                                               >  2:1002357/1‑100 (MQ=255)
                  tACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATaaaa                                                                              >  2:405295/1‑100 (MQ=255)
                  tACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATaaaa                                                                              >  2:623125/1‑100 (MQ=255)
                  tACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGAGAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATaaaa                                                                              >  2:1364436/1‑100 (MQ=255)
                  tACTGATGGCGCACCGGTGTTAAAACAAAAGCCTAGATAATTACCACCAACTATCGAGCCCACTAAAGTGTTGTATTGTGAGAGTGGCCCGAAGa                                                                                   <  2:998581/95‑1 (MQ=255)
                   aCTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATaaaaa                                                                             >  1:1019015/1‑100 (MQ=255)
                   aCTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATaaaaa                                                                             >  2:646524/1‑100 (MQ=255)
                    ctgctgGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAAt                                                                            <  2:2481204/100‑1 (MQ=255)
                    ctgctgGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAAt                                                                            >  1:2273020/1‑100 (MQ=255)
                     tgctgGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATa                                                                           >  2:609666/1‑100 (MQ=255)
                     tgctgGCGCACCGGCGTGAAAACAAAAGCCTAGATAGATACCGCCAACAATCGATCCTACTAAAAGGTTGTATTGTGACAGTGGCCCGAAGATaaaaa                                                                             <  2:1133550/98‑1 (MQ=255)
                      gctgGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATaa                                                                          >  1:1387129/1‑100 (MQ=255)
                      gctgGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATaa                                                                          >  2:1222960/1‑100 (MQ=255)
                      gctgGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATaa                                                                          <  1:982387/100‑1 (MQ=255)
                      gctgGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATaa                                                                          <  1:110842/100‑1 (MQ=255)
                       ctgGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATgaa                                                                         <  2:868002/100‑1 (MQ=255)
                         ggCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAg                                                                        <  1:104684/99‑1 (MQ=255)
                         ggCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGa                                                                       >  2:696705/1‑100 (MQ=255)
                           cgcACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAAc                                                                     >  2:151694/1‑100 (MQ=255)
                            gcACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACg                                                                    <  1:1935231/100‑1 (MQ=255)
                             cACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACg                                                                    >  1:1304724/1‑99 (MQ=255)
                              aCCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGgc                                                                  >  1:1469941/1‑100 (MQ=255)
                                   cGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCaaa                                                             >  2:769556/1‑100 (MQ=255)
                                   cGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCaaa                                                             <  1:691915/100‑1 (MQ=255)
                                    gTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGCTCCTACTAAAATGTTATATTCTGACAATGTCAAAAAGATAAAAATAAAGAATGACTCAAAc                                                            >  2:1576935/1‑100 (MQ=255)
                                    gTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAAc                                                            <  1:2296665/100‑1 (MQ=255)
                                    gTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAAc                                                            >  2:2067163/1‑100 (MQ=255)
                                    gTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAAc                                                            <  2:1633900/100‑1 (MQ=255)
                                     ttAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACa                                                           <  1:1033150/100‑1 (MQ=255)
                                      tAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACAt                                                          >  2:2385707/1‑100 (MQ=255)
                                      tAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACAt                                                          >  1:1921512/1‑100 (MQ=255)
                                       aaaaCAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATc                                                         <  1:1482857/100‑1 (MQ=255)
                                            aaaaGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACt                                                      <  2:792695/98‑1 (MQ=255)
                                            aaaaGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTaa                                                    >  1:1637289/1‑100 (MQ=255)
                                             aaaGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAAc                                                   <  1:983640/100‑1 (MQ=255)
                                              aaGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACa                                                  >  2:230510/1‑100 (MQ=255)
                                                 ccTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGc                                               >  1:556614/1‑100 (MQ=255)
                                                 ccTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGc                                               <  2:904460/100‑1 (MQ=255)
                                                  cTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGcc                                              <  1:1002365/100‑1 (MQ=255)
                                                  cTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGcc                                              <  1:1648655/100‑1 (MQ=255)
                                                   tAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCg                                             >  2:1553092/1‑100 (MQ=255)
                                                    aGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCgg                                            <  1:646531/100‑1 (MQ=255)
                                                       tAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGtaa                                         <  1:2072598/100‑1 (MQ=255)
                                                        aaaTACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGtaat                                        >  2:2123998/1‑100 (MQ=255)
                                                         aaTACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGtaata                                       >  1:919524/1‑100 (MQ=255)
                                                          aTACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGtaataa                                      <  2:332154/100‑1 (MQ=255)
                                                           tACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTaataat                                     >  1:1827923/1‑100 (MQ=255)
                                                             ccaccaACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATcc                                   <  1:1772327/100‑1 (MQ=255)
                                                             ccaccaACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATcc                                   <  1:643402/100‑1 (MQ=255)
                                                              caccaACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCca                                  >  1:1806713/1‑100 (MQ=255)
                                                              caccaACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCca                                  >  2:251974/1‑100 (MQ=255)
                                                              caccaACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCca                                  >  1:731014/1‑100 (MQ=255)
                                                               accaACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATcccc                                 <  1:1722548/100‑3 (MQ=255)
                                                                   acaaTCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCAcagca                             <  2:24055/100‑1 (MQ=255)
                                                                     aaTCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAgg                           <  2:1122101/100‑1 (MQ=255)
                                                                      aTCGATCCTACTAAAATGTTGTATTGTGACAGGGGCCCGAAGATAAAAATAAAGAAAGGCGCAAACATCACTAACATGCCGGTAATAATATACAGGAGGt                          >  2:1152886/1‑100 (MQ=255)
                                                                         gATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTAtt                       >  1:2095849/1‑100 (MQ=255)
                                                                         gATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTAtt                       >  1:858193/1‑100 (MQ=255)
                                                                         gATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGACGGTAATAATGGACAGCAGGTAt                        >  2:2500171/1‑99 (MQ=255)
                                                                             ctactaAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTgcg                   <  2:1889663/100‑1 (MQ=255)
                                                                              tactaAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTgcgc                  <  2:2273006/100‑1 (MQ=255)
                                                                              tactaAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTgcgc                  >  2:1466181/1‑100 (MQ=255)
                                                                                 taAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGc               <  1:228930/100‑1 (MQ=255)
                                                                                  aaaaTGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGcc              >  1:1058258/1‑100 (MQ=255)
                                                                                  aaaaTGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGcc              >  1:752712/1‑100 (MQ=255)
                                                                                  aaaaTGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGcc              >  1:81024/1‑100 (MQ=255)
                                                                                   aaaTGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGccc             <  2:209787/100‑1 (MQ=255)
                                                                                    aaTGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCg            <  2:1827909/100‑1 (MQ=255)
                                                                                    aaTGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCg            <  1:181408/100‑1 (MQ=255)
                                                                                     aTGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCATGTATTTGCGCAGCCCGa           >  2:1105712/1‑100 (MQ=255)
                                                                                      tgttgtATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAg          >  1:142595/1‑100 (MQ=255)
                                                                                       gttgtATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGt         >  1:2556413/1‑100 (MQ=255)
                                                                                        ttgtattgtgACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGtt        <  2:556607/100‑1 (MQ=255)
                                                                                        ttgtattgagACAGTGGCCAGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGt         <  2:886241/99‑1 (MQ=255)
                                                                                         tgtattgtgACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGttt       <  2:1896743/100‑1 (MQ=255)
                                                                                          gtattgtgACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTg      >  1:1293609/1‑100 (MQ=255)
                                                                                            attgtgACAGTGGCCCGAAGATAAAAATAAAGAACGGCTCAAACATCACTAACATGCCGGTAATATTCCACAGCAGGTATTTGCGCAGCCCGAGTTTGTc    <  2:920982/100‑1 (MQ=255)
                                                                                            attgtgACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTGTc    <  1:2128817/100‑1 (MQ=255)
                                                                                             ttgtgACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTGTCa   >  1:2140020/1‑100 (MQ=255)
                                                                                             ttgtgACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTGTCa   <  1:2170697/100‑1 (MQ=255)
                                                                                             ttgtgACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTGTCa   >  1:2502586/1‑100 (MQ=255)
                                                                                             ttgtgACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTGTCa   <  1:1578891/100‑1 (MQ=255)
                                                                                              tgtgACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTGTCAg  >  1:1205095/1‑100 (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 < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: