Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 889,923 A→G intergenic (‑104/‑166) ybjL ← / → ybjM putative transport protein YbjL/putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913889,9230AG100.0% 253.0 / NA 67intergenic (‑104/‑166)ybjL/ybjMputative transport protein YbjL/putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/45);  total (22/45)

CCATTTAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAACGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACGCCACACTGGCGAAA  >  NC_000913/889788‑890053
                                                                                                                                        |                                                                                                                                  
ccATTTAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGCTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGAt                                                                                                                                <  1:155088/140‑1 (MQ=255)
ccATTTAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGCTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGAt                                                                                                                                >  2:155088/1‑140 (MQ=255)
 cATTTAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATa                                                                                                                               >  2:37357/1‑140 (MQ=255)
 cATTTAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATAGTAATTTACGATAATGTTTTTTACCAGAATTCAGGTCGCAGATTCATTCAGCGCTCCTAAGCGTTa                                                                                                                               <  2:214575/140‑1 (MQ=255)
  aTTTAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATag                                                                                                                              <  2:142074/140‑1 (MQ=255)
  aTTTAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATag                                                                                                                              >  2:166867/1‑140 (MQ=255)
   tttAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATagt                                                                                                                             >  2:552101/1‑140 (MQ=255)
    ttAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATagta                                                                                                                            <  1:837476/140‑1 (MQ=255)
        cAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATagtaaagt                                                                                                                        <  1:779103/140‑1 (MQ=255)
          aTTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATagtaaagtaa                                                                                                                      <  1:209503/140‑1 (MQ=255)
           ttCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAAc                                                                                                                     >  2:466265/1‑140 (MQ=255)
            tCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACa                                                                                                                    <  1:309995/140‑1 (MQ=255)
            tCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACa                                                                                                                    <  2:747818/140‑1 (MQ=255)
            tCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACa                                                                                                                    <  2:8797/140‑1 (MQ=255)
             cGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACaa                                                                                                                   <  2:243597/140‑1 (MQ=255)
             cGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACaa                                                                                                                   <  1:368251/140‑1 (MQ=255)
                 gACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAatata                                                                                                               <  1:157426/140‑1 (MQ=255)
                   cGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATAtt                                                                                                             >  2:441559/1‑140 (MQ=255)
                   cGTTTATATTCACGGAGGCTAACTTCTTGTTTACAGGAAGGTGGTGGAAGAAATGGGAATTTAAGAAAATGTTTTTTACCAAAATTCAGGGCGCAAGTTTATTTAACCGCACTTAGCCATAAGAAAGTAACAATATtttt                                                                                                            >  2:748910/1‑140 (MQ=255)
                       tataTTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTAc                                                                                                         <  2:916153/140‑1 (MQ=255)
                           ttCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAgt                                                                                                     <  1:166867/140‑1 (MQ=255)
                            tCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAgtg                                                                                                    <  2:588158/140‑1 (MQ=255)
                             cACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAgtgt                                                                                                   >  2:482745/1‑140 (MQ=255)
                              aCGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTa                                                                                                  >  2:495320/1‑140 (MQ=255)
                              aCGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTa                                                                                                  <  1:482745/140‑1 (MQ=255)
                               cGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAGTTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTaa                                                                                                 >  2:873808/1‑140 (MQ=255)
                               cGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTaa                                                                                                 <  1:638193/140‑1 (MQ=255)
                                 gAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATc                                                                                               <  2:241431/140‑1 (MQ=255)
                                   ggCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATcac                                                                                             >  2:828270/1‑140 (MQ=255)
                                    gCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATcaca                                                                                            <  2:651927/140‑1 (MQ=255)
                                    gCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATcaca                                                                                            <  1:839034/140‑1 (MQ=255)
                                        aCTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAg                                                                                        <  1:466265/140‑1 (MQ=255)
                                          ttcttGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGt                                                                                      <  2:366886/140‑1 (MQ=255)
                                               gTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCaa                                                                                 <  1:447930/140‑1 (MQ=255)
                                               gTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCaa                                                                                 >  2:604718/1‑140 (MQ=255)
                                                      aGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCtat                                                                          <  1:441559/140‑1 (MQ=255)
                                                          aGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGa                                                                      >  2:902139/1‑140 (MQ=255)
                                                           gCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGaa                                                                     <  2:164821/140‑1 (MQ=255)
                                                           gCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGaa                                                                     <  1:139432/140‑1 (MQ=255)
                                                               ttGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGTGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGc                                                                 >  1:403210/1‑140 (MQ=255)
                                                               ttGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGc                                                                 >  1:364161/1‑140 (MQ=255)
                                                               ttGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGc                                                                 >  2:313713/1‑140 (MQ=255)
                                                                tGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCg                                                                <  1:902139/140‑1 (MQ=255)
                                                                tGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCg                                                                <  2:483969/140‑1 (MQ=255)
                                                                 gaaagaaaTGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGt                                                               <  2:574070/140‑1 (MQ=255)
                                                                    agaaaTGGTAATTTACGATAATGTTTTTTACCAGAGTTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTgg                                                            <  1:873808/140‑1 (MQ=255)
                                                                    agaaaTGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTgg                                                            <  1:146621/140‑1 (MQ=255)
                                                                    agaaaTGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTgg                                                            <  1:748910/140‑1 (MQ=255)
                                                                    agaaaTGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCACACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTgg                                                            <  1:604718/140‑1 (MQ=255)
                                                                     gaaaTGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGc                                                           <  2:31400/140‑1 (MQ=255)
                                                                     gaaaTGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGc                                                           <  2:312257/140‑1 (MQ=255)
                                                                     gaaaTGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGc                                                           <  1:37357/140‑1 (MQ=255)
                                                                       aaTGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCt                                                         >  2:375688/1‑140 (MQ=255)
                                                                         tGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCtat                                                       >  1:416372/1‑140 (MQ=255)
                                                                          ggTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCtata                                                      <  2:416372/140‑1 (MQ=255)
                                                                                 tACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCAt                                               >  2:558608/1‑140 (MQ=255)
                                                                                           ttttttACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGaa                                     <  2:364161/140‑1 (MQ=255)
                                                                                           ttttttACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGaa                                     <  2:653394/140‑1 (MQ=255)
                                                                                           ttttttACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGaa                                     <  2:274164/140‑1 (MQ=255)
                                                                                               ttACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCgg                                 <  1:828270/140‑1 (MQ=255)
                                                                                                 aCCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCt                               <  2:647237/140‑1 (MQ=255)
                                                                                                 aCCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCt                               <  2:876562/140‑1 (MQ=255)
                                                                                                           aGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACt                     <  1:633436/140‑1 (MQ=255)
                                                                                                                gcAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACg                >  1:731338/1‑140 (MQ=255)
                                                                                                                   aGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACGcc              >  1:68303/1‑140 (MQ=255)
                                                                                                                      ttcattcaGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACGCcaca           <  2:574403/140‑1 (MQ=255)
                                                                                                                      ttcattcaGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACGCcaca           <  1:123857/140‑1 (MQ=255)
                                                                                                                               cgcACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACGCCACACTGGCGaaa  >  1:420416/1‑140 (MQ=255)
                                                                                                                               cgcACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACGCCACACTGGCGaaa  >  1:589926/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                  
CCATTTAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAACGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACGCCACACTGGCGAAA  >  NC_000913/889788‑890053

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