Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,049,519 A→C K1528Q (AAA→CAA)  yeeJ → putative adhesin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,049,5190AC100.0% 230.9 / NA 61K1528Q (AAA→CAA) yeeJputative adhesin
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (27/34);  total (27/34)

AGCAACAGTCGCTGATACAGAGGGCAATGCGATAGCCAACACTGAGGTAACATTTACTCTGCCGGAAGATGTGAAGGCGAACTTCACGCTGAGCGATGGCGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGAAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTGAGCAGTTGGTGGTGAACTTTATTGCGGATACGCTCACTGCGCAGGTTAATCTTAACGTTACCGAGGACAATTTTATCGCTAATAACGTC  >  NC_000913/2049374‑2049665
                                                                                                                                                 |                                                                                                                                                  
aGCAACAGTCGCTGATACAGAGGGCAATGCGATAGCCAACACTGAGGTAACATTTACTCTGCCGGAAGATGTGAAGGCGAACTTCACGCTGAGCGATGGCGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGt                                                                                                                                               >  2:104202/1‑151 (MQ=255)
    aCAGTCGCTGATACAGAGGGCAATGCGATAGCCAACACTGAGGTAACATTTACTCTGCCGGAAGATGTGAAGGCGAACTTCACGCTGAGCGATGGCGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACaa                                                                                                                                           <  1:828924/151‑1 (MQ=255)
      aGTCGCTGATACAGAGGGCAATGCGATAGCCAACACTGAGGTAACATTTACTCTGCCGGAAGATGTGAAGGCGAACTTCACGCTGAGCGATGGCGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACaaa                                                                                                                                          >  1:273216/1‑150 (MQ=255)
      aGTCGCTGATACAGAGGGCAATGCGATAGCCAACACTGAGGTAACATTTACTCTGCCGGAAGATGTGAAGGCGAACTTCACGCTGAGCGATGGCGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACaaa                                                                                                                                          >  2:591134/1‑150 (MQ=255)
       gTCGCTGATACAGAGGGCAATGCGATAGCCAACACTGAGGTAACATTTACTCTGCCGGAAGATGTGAAGGCGAACTTCACGCTGAGCGATGGCGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAg                                                                                                                                        <  2:4516/151‑1 (MQ=255)
       gTCGCTGATACAGAGGGCAATGCGATAGCCAACACTGAGGTAACATTTACTCTGCCGGAAGATGTGAAGGCGAACTTCACGCTGAGCGATGGCGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAg                                                                                                                                        <  2:417382/151‑1 (MQ=255)
          gCTGATACAGAGGGCAATGCGATAGCCAACACTGAGGTAACATTTACTCTGCCGGAAGATGTGAAGGCGAACTTCACGCTGAGCGATGGCGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGc                                                                                                                                       >  2:519876/1‑149 (MQ=255)
             gATACAGAGGGCAATGCGATAGCCAACACTGAGGTAACATTTACTCTGCCGGAAGATGTGAAGGCGAACTTCACGCTGAGCGATGGCGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGgcg                                                                                                                                  >  1:595292/1‑151 (MQ=255)
                aCAGAGGGCAATGCGATAGCCAACACTGAGGTAACATTTACTCTGCCGGAAGATGTGAAGGCGAACTTCACGCTGAGCGATGGCGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTc                                                                                                                               <  1:503470/151‑1 (MQ=255)
                 cAGAGGGCAATGCGATAGCCAACACTGAGGTAACATTTACTCTGCCGGAAGATGTGAAGGCGAACTTCACGCTGAGCGATGGCGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTc                                                                                                                               >  1:473460/1‑150 (MQ=255)
                    aaggCAATGCGATAGCCAACACTGAGGTAACATTTACTCTGCCGGAAGATGTGAAGGCGAACTTCACGCTGAGCGATGGCGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATAc                                                                                                                           <  2:315141/149‑1 (MQ=255)
                     gggCAATGCGATAGCCAACACTGAGGTAACATTTACTCTGCCGGAAGATGTGAAGGCGAACTTCACGCTGAGCGATGGCGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACt                                                                                                                          <  2:668845/151‑1 (MQ=255)
                       gCAATGCGATAGCCAACACTGAGGTAACATTTACTCTGCCGGAAGATGTGAAGGCGAACTTCACGCTGAGCGATGGCGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGt                                                                                                                        >  2:69361/1‑151 (MQ=255)
                             cGATAGCCAACACTGAGGTAACATTTACTCTGCCGGAAGATGTGAAGGCGAACTTCACGCTGAGCGATGGCGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGc                                                                                                                  >  2:785929/1‑151 (MQ=255)
                              gATAGCCAACACTGAGGTAACATTTACTCTGCCGGAAGATGTGAAGGCGAACTTCACGCTGAGCGATGGCGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCa                                                                                                                 <  1:384319/151‑1 (MQ=255)
                                 aGCCAACACTGAGGTAACATTTACTCTGCCGGAAGATGTGAAGGCGAACTTCACGCTGAGCGATGGCGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATc                                                                                                               <  2:455857/150‑1 (MQ=255)
                                 aGCCAACACTGAGGTAACATTTACTCTGCCGGAAGATGTGAAGGCGAACTTCACGCTGAGCGATGGCGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCg                                                                                                              <  2:704625/151‑1 (MQ=255)
                                                  cATTTACTCTGCCGGAAGATGTGAAGGCGAACTTCACGCTGAGCGATGGCGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATTGATGACTGGCGGTAagag                                                                                             <  1:369156/151‑1 (MQ=255)
                                                   aTTTACTCTGCCGGAAGATGTGAAGGCGAACTTCACGCTGAGCGATGGCGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAagag                                                                                             >  2:626481/1‑150 (MQ=255)
                                                     ttACTCTGCCGGAAGATGTGAAGGCGAACTTCACGCTGAGCGATGGCGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTg                                                                                           <  2:533818/150‑1 (MQ=255)
                                                     ttACTCTGCCGGAAGATGTGAAGGCGAACTTCACGCTGAGCGATGGCGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTg                                                                                           >  1:841951/1‑150 (MQ=255)
                                                             ccGGAAGATGTGAAGGCGAACTTCACGCTGAGCGATGGCGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTGAGCAGTtgg                                                                                  <  1:781592/151‑1 (MQ=255)
                                                               ggAAGATGTGAAGGCGAACTTCACGCTGAGCGATGGCGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGGGAGCAGTtggtg                                                                                >  1:651113/1‑151 (MQ=255)
                                                                   gATGTGAAGGCGAACTTCACGCTGAGCGATGGCGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTGAGCAGTTGGTGGTGa                                                                            <  2:595631/151‑1 (MQ=255)
                                                                      gtgAAGGCGAACTTCACGCTGAGCGATGGCGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTGAGCAGTTGGTGGTGAAc                                                                          >  1:586958/1‑150 (MQ=255)
                                                                        gAAGGCGAACTTCACGCTGAGCGATGGCGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTGAGCAGTTGGTGGTGAACttt                                                                       <  2:551861/151‑1 (MQ=255)
                                                                        gAAGGCGAACTTCACGCTGAGCGATGGCGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTGAGCAGTTGGTGGTGAACttt                                                                       <  2:478391/151‑1 (MQ=255)
                                                                         aaGGCGAACTTCACGCTGAGCGATGGCGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTGAGCAGTTGGTGGTGAACttt                                                                       >  2:80608/1‑150 (MQ=255)
                                                                          aGGCGAACTTCACGCTGAGCGATGGCGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTGAGCAGTTGGTGGTGAACTTTAt                                                                     <  1:677118/151‑1 (MQ=255)
                                                                             cGAACTTCACGCTGAGCGATGGCGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTGAGCAGTTGGTGGTGAACTTTATTGc                                                                  <  2:629753/151‑1 (MQ=255)
                                                                                    cACGCTGAGCGATGGCGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTGAGCAGTTGGTGGTGAACTTTATTGCGGATAc                                                            <  2:854009/150‑1 (MQ=255)
                                                                                    cACGCTGAGCGATGGCGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTGAGCAGTTGGTGGTGAACTTTATTGCGGATACg                                                           <  2:514274/151‑1 (MQ=255)
                                                                                     aCGCTGAGCGATGGCGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTGAGCAGTTGGTGGTGAACTTTATTGCGGATACGc                                                          <  1:764183/151‑1 (MQ=255)
                                                                                     aCGCTGAGCGATGGCGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTGAGCAGTTGGTGGTGAACTTTATTGCGGATACGc                                                          >  2:452070/1‑151 (MQ=255)
                                                                                       gCTGAGCGATGGCGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTGAGCAGTTGGTGGTGAACTTTATTGCGGATACGCTc                                                        >  2:776066/1‑151 (MQ=255)
                                                                                            gCGATGGCGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTGAGCAGTTGGTGGTGAACTTTATTGCGGATACGCTCACTgc                                                   <  1:702074/151‑1 (MQ=255)
                                                                                                 ggCGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTGAGCAGTTGGTGGTGAACTTTATTGCGGATACGCTCACTGCGCAgg                                              >  2:357275/1‑151 (MQ=255)
                                                                                                 ggCGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTGAGCAGTTGGTGGTGAACTTTATTGCGGATACGCTCACTGCGCAg                                               <  1:154393/150‑1 (MQ=255)
                                                                                                   cGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTGAGCAGTTGGTGGTGAACTTTATTGCGGATACGCTCACTGCGCAGGtt                                            >  2:239916/1‑151 (MQ=255)
                                                                                                            gATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTGAGCAGTTGGTGGTGAACTTTATTGCGGATACGCTCACTGCGCAGGTTAATCTTaa                                    <  1:252961/150‑1 (MQ=255)
                                                                                                            gATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTGAGCAGTTGGTGGTGAACTTTATTGCGGATACGCTCACTGCGCAGGTTAATCTTAAc                                   >  2:820295/1‑151 (MQ=255)
                                                                                                             aTTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTGAGCAGTTGGTGGTGAACTTTATTGCGGATACGCTCACTGCGCAGGTTAATCTTAACg                                  >  1:872475/1‑151 (MQ=255)
                                                                                                              ttACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTGAGCAGTTGGTGGTGAACTTTATTGCGGATACGCTCACTGCGCAGGTTAATCTTaa                                    >  1:82415/1‑148 (MQ=255)
                                                                                                               tACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTGAGCAGTTTGTGGTGAACTTTATTGCGGATACGCTCACTGCGCAGGTTAATCTTAACGtt                                >  2:335137/1‑151 (MQ=255)
                                                                                                                 cTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTGAGCAGTTGGTGGTGAACTTTATTGCGGATACGCTCACTGCGCAGGTTAATCTTAACGTTa                               >  2:805830/1‑150 (MQ=255)
                                                                                                                 cTGATGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTGAGCAGTTGGTGGTGAACTTTATTGCGGATACGCTCACTGCGCAGGTTAATCTTAACGTTAc                              >  1:399016/1‑151 (MQ=255)
                                                                                                                    aTGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTGAGCAGTTGGTGGTGAACTTTATTGCGGATACGCTCACTGCGCAGGTTAATCTTAACGTTACCg                            >  1:693090/1‑150 (MQ=255)
                                                                                                                    aTGCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTGAGCAGTTGGTGGTGAACTTTATTGCGGATACGCTCACTGCGCAGGTTAATCTTAACGTTACCGa                           <  2:359612/151‑1 (MQ=255)
                                                                                                                      gCTGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTGAGTAGTTGGTGGTGAACTTTATTGCGGATACGCTCACTGCGCAGGTTAATCTTAACGTTACCGAg                          <  1:428686/150‑1 (MQ=255)
                                                                                                                        tGAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTGAGCAGTTGGTGGTGAACTTTATTGCGGATACGCTCACTGCGCAGGTTAATCTTAACGTTACCGAGGAc                       >  2:664966/1‑151 (MQ=255)
                                                                                                                         gAAGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTGAGCAGTTGGTGGTGAACTTTATTGCGGATACGCTCACTGCGCAGGTTAATCTTAACGTTACCGAGGACa                      <  1:519859/151‑1 (MQ=255)
                                                                                                                          aaGGCAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTGAGCAGTTGGTGGTGAACTTTATTGCGGATACGCTCACTGCGCAGGTTAATCTTAACGTTACCGAGGACaa                     <  1:69359/151‑1 (MQ=255)
                                                                                                                              cAAAGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTGAGCAGTTGGTGGTGAACTTTATTGCGGATACGCTCACTGCGCAGGTTAATCTTAACGTTACCGAGGACAAtttt                 <  2:377730/151‑1 (MQ=255)
                                                                                                                               aaaGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTGAGCAGTTGGTGGTGAACTTTATTGCGGATACGCTCACTGCGCAGGTTAATCTTAACGTTACCGAGGACAAtttt                 >  2:337845/1‑150 (MQ=255)
                                                                                                                                aaGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTGAGCAGTTGGTGGTGAACTTTATTGCGGATACGCTCACTGCGCAGGTTAATCTTAACGTTACCGAGGACAATTTTAt               <  1:822571/151‑1 (MQ=255)
                                                                                                                                 aGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTGAGCAGTTGGTGGTGAACTTTATTGCGGATACGCTCACTGCGCAGGTTAATCTTAACGTTACCGAGGACAATTTTATc              <  2:679155/151‑1 (MQ=255)
                                                                                                                                 aGCGAAAGTCACGCTGCAAGGTACAAAAGCAGGCGATCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTGAGCAGTTGGTGGTGAACTTTATTGCGGATACGCTCACTGCGCAGGTTAATCTTAACGTTCCCGAGGACAATTTTATc              <  1:337498/151‑1 (MQ=255)
                                                                                                                                       aGTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGAAAGAGTGAGCAGTTGGTGGTGAACTTTATTGCGGATACGCTCACTGCGCAGGTTAATCTTAACGTTACCGAGGACAATTTTATCGCtaat        >  2:308875/1‑151 (MQ=255)
                                                                                                                                        gTCACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTGAGCAGTTGGTGGTGAACTTTATTGCGGATACGCTCACTGCGCAGGTTAATCTTAACGTTACCGAGGACAATTTTATCGCtaata       <  2:473475/151‑1 (MQ=255)
                                                                                                                                          cACGCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTGAGCAGTTGGTGGTGAACTTTATTGCGGATACGCTCACTGCGCAGGTTAATCTTAACGTTACCGAGGACAATTTTATCGCtaataa      <  1:791270/150‑1 (MQ=255)
                                                                                                                                             gCTGCAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTGAGCAGTTGGTGGTGAACTTTATTGCGGATACGCTCACTGCGCAGGTTAATCTTAACGTTACCGAGGACAATTTTATCGCTAATAACGTc  <  1:631331/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
AGCAACAGTCGCTGATACAGAGGGCAATGCGATAGCCAACACTGAGGTAACATTTACTCTGCCGGAAGATGTGAAGGCGAACTTCACGCTGAGCGATGGCGGTAAAGTGATTACTGATGCTGAAGGCAAAGCGAAAGTCACGCTGAAAGGTACAAAAGCAGGCGCTCATACTGTTACAGCATCGATGACTGGCGGTAAGAGTGAGCAGTTGGTGGTGAACTTTATTGCGGATACGCTCACTGCGCAGGTTAATCTTAACGTTACCGAGGACAATTTTATCGCTAATAACGTC  >  NC_000913/2049374‑2049665

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