Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 470,643 G→A *591* (TAG→TAA) 
S3N (AGT→AAT) 		 mdlA →
mdlB →		 putative multidrug ABC transporter ATPase
putative multidrug ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913470,6430GA100.0% 250.4 / NA 77*591* (TAG→TAA) 
S3N (AGT→AAT) 		mdlA
mdlB		putative multidrug ABC transporter ATPase
putative multidrug ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (43/34);  total (43/34)

GACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTGCGGCGGCGGCAGAAGTCAGTGGGCCGCTGCTTATCAGCT  >  NC_000913/470507‑470777
                                                                                                                                        |                                                                                                                                      
gACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAAtt                                                                                                                                      >  2:791907/1‑139 (MQ=255)
     atCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGc                                                                                                                                 <  2:939893/139‑1 (MQ=255)
      tCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGcc                                                                                                                                >  2:52685/1‑139 (MQ=255)
      tCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGAAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGcc                                                                                                                                >  2:1108660/1‑139 (MQ=255)
          ccAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACt                                                                                                                            >  2:1561101/1‑139 (MQ=255)
           cAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACtg                                                                                                                           >  1:1324448/1‑139 (MQ=255)
              cGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTgg                                                                                                                        >  1:1180295/1‑139 (MQ=255)
              cGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTgg                                                                                                                        >  1:887311/1‑139 (MQ=255)
                   cAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGc                                                                                                                       >  2:458741/1‑135 (MQ=255)
                   cAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGc                                                                                                                       <  1:458741/135‑1 (MQ=255)
                                cTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTg                                                                                                      >  1:425545/1‑139 (MQ=255)
                                  ggCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGtt                                                                                                    >  1:1395282/1‑139 (MQ=255)
                                    cacaACAAAGCGGCTGGTATCGCGATATGTATCCCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAg                                                                                                  >  2:1208754/1‑139 (MQ=255)
                                        acaaAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTa                                                                                              <  2:1731539/139‑1 (MQ=255)
                                        acaaAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTa                                                                                              >  1:1731539/1‑139 (MQ=255)
                                            aGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGt                                                                                          <  2:1180295/139‑1 (MQ=255)
                                                    gTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTg                                                                                  <  2:1324448/139‑1 (MQ=255)
                                                          cgATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCATCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTaa                                                                            >  1:923542/1‑139 (MQ=255)
                                                           gATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTaaa                                                                           <  1:1598174/139‑1 (MQ=255)
                                                              atGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCaa                                                                                                             <  2:744826/102‑1 (MQ=255)
                                                              atGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCaa                                                                                                             >  1:744826/1‑102 (MQ=255)
                                                              atGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCg                                                                        >  1:1650323/1‑139 (MQ=255)
                                                               tGTATCGCTATCAACAAGTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTaaa                                                                           >  2:588420/1‑135 (MQ=255)
                                                               tGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTaaa                                                                           <  1:588420/135‑1 (MQ=255)
                                                                   tCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTa                                                                                              >  2:1073413/1‑112 (MQ=255)
                                                                   tCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTa                                                                                              <  1:1073413/112‑1 (MQ=255)
                                                                   tCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTggg                                                                   <  1:52685/139‑1 (MQ=255)
                                                                       tATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGAtt                                                               >  2:1749121/1‑139 (MQ=255)
                                                                        aTCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGc                                                                                                                       >  2:611072/1‑82 (MQ=255)
                                                                        aTCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGc                                                                                                                       <  1:611072/82‑1 (MQ=255)
                                                                        aTCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTg                                                              <  2:156218/139‑1 (MQ=255)
                                                                        aTCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTg                                                              <  2:875142/139‑1 (MQ=255)
                                                                         tCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGc                                                             >  1:571950/1‑139 (MQ=255)
                                                                         tCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGc                                                             >  1:373426/1‑139 (MQ=255)
                                                                           aacaacTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCgg                                                           >  2:91545/1‑139 (MQ=255)
                                                                             caacTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTc                                                         <  1:1208754/139‑1 (MQ=255)
                                                                                cTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCtg                                                      <  1:1749121/139‑1 (MQ=255)
                                                                                   gAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAAttt                                                                                                                                     <  1:1689510/57‑1 (MQ=255)
                                                                                   gAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAAttt                                                                                                                                     >  2:1689510/1‑57 (MQ=255)
                                                                                   gAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTg                                                                                                      >  2:1178252/1‑88 (MQ=255)
                                                                                   gAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTg                                                                                                      <  1:1178252/88‑1 (MQ=255)
                                                                                   gAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCtgatg                                                   >  2:267555/1‑139 (MQ=255)
                                                                                   gAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCtgatg                                                   >  1:386590/1‑139 (MQ=255)
                                                                                   gAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCtgatg                                                   >  2:1687814/1‑139 (MQ=255)
                                                                                     ggcggcGCTCGACGACGCTCCGGACAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGcc                                                                                                        <  1:1186176/84‑1 (MQ=255)
                                                                                     ggcggcGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGcc                                                                                                        >  2:1186176/1‑84 (MQ=255)
                                                                                     ggcggcGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCgg                                                           >  1:728943/1‑129 (MQ=255)
                                                                                     ggcggcGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCgg                                                           <  2:728943/129‑1 (MQ=255)
                                                                                     ggcggcGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTgatgat                                                 >  2:413151/1‑139 (MQ=255)
                                                                                     ggcggcGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTgatgat                                                 >  1:1473761/1‑139 (MQ=255)
                                                                                     ggcggcGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTgatgat                                                 <  2:571950/139‑1 (MQ=255)
                                                                                          cgcTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAg                                                                                                                                  <  2:311929/53‑1 (MQ=255)
                                                                                          cgcTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAg                                                                                                                                  >  1:311929/1‑53 (MQ=255)
                                                                                          cgcTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTa                                                                             <  1:192749/106‑1 (MQ=255)
                                                                                          cgcTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTa                                                                             >  2:192749/1‑106 (MQ=255)
                                                                                           gcTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGc                                                                                                                                 <  1:498148/53‑1 (MQ=255)
                                                                                           gcTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGc                                                                                                                                 >  2:498148/1‑53 (MQ=255)
                                                                                           gcTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGt                                           >  2:57699/1‑139 (MQ=255)
                                                                                             tCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTg                                         <  1:1624139/139‑1 (MQ=255)
                                                                                               gacgacGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTgcg                                       >  1:224925/1‑139 (MQ=255)
                                                                                                 cgacgCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTgcggc                                     <  1:791907/139‑1 (MQ=255)
                                                                                                       tCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGc                                                             <  2:1161641/109‑1 (MQ=255)
                                                                                                       tCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGc                                                             >  1:1161641/1‑109 (MQ=255)
                                                                                                                     aggaggCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGATTGCGGCGGCGGCAGAAGTCAGTGGGc                 >  1:1035417/1‑139 (MQ=255)
                                                                                                                       gaggCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTGCGGCGGCGGCAGAAGTCAGTGGGCCg               <  1:57699/139‑1 (MQ=255)
                                                                                                                        aggCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCgg                                                           <  2:477000/94‑1 (MQ=255)
                                                                                                                        aggCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCgg                                                           >  1:477000/1‑94 (MQ=255)
                                                                                                                          gCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGa                                                                                                                    <  2:300721/35‑1 (MQ=255)
                                                                                                                          gCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGa                                                                                                                    >  1:300721/1‑35 (MQ=255)
                                                                                                                          gCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTa                                                                             >  2:1722074/1‑74 (MQ=255)
                                                                                                                          gCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTa                                                                             <  1:1722074/74‑1 (MQ=255)
                                                                                                                          gCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTGCGGCGGCGGCAGAAGTCAGTGGGCCgctg            >  1:195377/1‑139 (MQ=255)
                                                                                                                                gATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTTCTGATGATGTGGGTTGCGGCGGCGGCAGAAGTCAGTGGGCCGCTGCTTATc      <  1:1108660/139‑1 (MQ=255)
                                                                                                                                gATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTGCGGCGGCGGCAGAAGTCAGTGGGCCGCTGCTTATc      <  2:1473761/139‑1 (MQ=255)
                                                                                                                                 aTGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTGCGGCGGCGGCAGAAGTCAGTGGGCCGCTGCTTATCa     >  2:849962/1‑139 (MQ=255)
                                                                                                                                 aTGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTGCGGCGGCGGCAGAAGTCAGTGGGCCGCTGCTTATCa     <  2:1395282/139‑1 (MQ=255)
                                                                                                                                    cGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTGCGGCGGCGGCAGAAGTCAGTGGGCCGCTGCTTATCAGCt  <  2:386590/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                      
GACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTGCGGCGGCGGCAGAAGTCAGTGGGCCGCTGCTTATCAGCT  >  NC_000913/470507‑470777

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