Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 469,429 A→G K187E (AAG→GAG)  mdlA → putative multidrug ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913469,4290AG100.0% 302.7 / NA 85K187E (AAG→GAG) mdlAputative multidrug ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (46/39);  total (46/39)

CTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCG  >  NC_000913/469305‑469561
                                                                                                                            |                                                                                                                                    
cTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAggcg                                                                                                                        <  1:1311121/139‑1 (MQ=255)
   cgcTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGgcggcg                                                                                                                     >  2:1386671/1‑139 (MQ=255)
          ttGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCa                                                                                                              <  1:1386671/139‑1 (MQ=255)
                      tCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGAcc                                                                                                  >  1:439829/1‑139 (MQ=255)
                       cTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGgcggcg                                                                                                                     <  1:119757/119‑1 (MQ=255)
                       cTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGgcggcg                                                                                                                     >  2:119757/1‑119 (MQ=255)
                       cTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCg                                                                                                 >  2:1793809/1‑139 (MQ=255)
                        tACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGACACAGGCGGCGTTTTCCa                                                                                                              >  2:29231/1‑125 (MQ=255)
                        tACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGACACAGGCGGCGTTTTCCa                                                                                                              <  1:29231/125‑1 (MQ=255)
                         aCGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCa                                                                                               >  2:814527/1‑139 (MQ=255)
                                ttAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGa                                                                                        <  1:917273/139‑1 (MQ=255)
                                ttAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGa                                                                                        >  2:1406363/1‑139 (MQ=255)
                                 tAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGaa                                                                                       <  2:439829/139‑1 (MQ=255)
                                  aGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGaaa                                                                                      >  2:765229/1‑139 (MQ=255)
                                       gCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCa                                                                                                              <  2:271941/110‑1 (MQ=255)
                                       gCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCa                                                                                                              >  1:271941/1‑110 (MQ=255)
                                         aGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCAc                                                                               >  1:46886/1‑139 (MQ=255)
                                         aGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCAc                                                                               <  2:46886/139‑1 (MQ=255)
                                         aGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCAc                                                                               <  1:153068/139‑1 (MQ=255)
                                              aCCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCAc                                                                               >  1:478962/1‑134 (MQ=255)
                                              aCCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCAc                                                                               <  2:478962/134‑1 (MQ=255)
                                                 ttATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATcc                                                                       >  1:1450339/1‑139 (MQ=255)
                                                  tATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCg                                                                      >  1:177753/1‑139 (MQ=255)
                                                   aTTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAg                                                                                     <  1:1744738/123‑1 (MQ=255)
                                                   aTTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAg                                                                                     >  2:1744738/1‑123 (MQ=255)
                                                    ttttCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCa                                                                    <  1:765229/139‑1 (MQ=255)
                                                         ccTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGtttt                                                                                                                 >  2:1659160/1‑89 (MQ=255)
                                                         ccTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGtttt                                                                                                                 <  1:1659160/89‑1 (MQ=255)
                                                         ccTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGAc                                                                                                   <  1:668485/103‑1 (MQ=255)
                                                         ccTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGAc                                                                                                   >  2:668485/1‑103 (MQ=255)
                                                           tGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAgg                                                                                         <  1:1858394/111‑1 (MQ=255)
                                                           tGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAgg                                                                                         >  2:1858394/1‑111 (MQ=255)
                                                             ttGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTaa                                                                                                       >  2:1666715/1‑95 (MQ=255)
                                                             ttGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTaa                                                                                                       <  1:1666715/95‑1 (MQ=255)
                                                             ttGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCt                                                                                  <  1:715365/116‑1 (MQ=255)
                                                             ttGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCt                                                                                  >  2:715365/1‑116 (MQ=255)
                                                             ttGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAg                                                           >  2:1266582/1‑139 (MQ=255)
                                                               gCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTc                                                                                                                <  1:224543/84‑1 (MQ=255)
                                                               gCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTc                                                                                                                >  2:224543/1‑84 (MQ=255)
                                                                ccGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGAcc                                                                                                  >  1:176174/1‑97 (MQ=255)
                                                                ccGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGAcc                                                                                                  <  2:176174/97‑1 (MQ=255)
                                                                 cGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCtt                                                       <  1:1793809/139‑1 (MQ=255)
                                                                   aTGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAg                                                                                                                                    <  2:674115/60‑1 (MQ=255)
                                                                   aTGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAg                                                                                                                                    >  1:674115/1‑60 (MQ=255)
                                                                       cAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAggcg                                                                                                                        >  1:254884/1‑68 (MQ=255)
                                                                       cAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAggcg                                                                                                                        <  2:254884/68‑1 (MQ=255)
                                                                       cAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCt                                                 <  2:865160/139‑1 (MQ=255)
                                                                          tGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGgcggcg                                                                                                                     <  1:51732/68‑1 (MQ=255)
                                                                          tGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGgcggcg                                                                                                                     >  2:51732/1‑68 (MQ=255)
                                                                               gCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATc                                         >  2:33233/1‑139 (MQ=255)
                                                                                  aTCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGcc                                      >  2:340713/1‑139 (MQ=255)
                                                                                  aTCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGcc                                      >  1:704604/1‑139 (MQ=255)
                                                                                   tCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCa                                     >  1:892820/1‑139 (MQ=255)
                                                                                    cATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCAc                                                                                              >  2:1401472/1‑81 (MQ=255)
                                                                                    cATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCAc                                                                                              <  1:1401472/81‑1 (MQ=255)
                                                                                    cATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAg                                    <  1:340713/139‑1 (MQ=255)
                                                                                       gATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCg                                 <  2:1450339/139‑1 (MQ=255)
                                                                                       gATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCg                                 <  1:33233/139‑1 (MQ=255)
                                                                                        aTTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAgg                                                                                         >  1:1308748/1‑82 (MQ=255)
                                                                                        aTTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAgg                                                                                         <  2:1308748/82‑1 (MQ=255)
                                                                                          tAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTc                                                                                                                >  2:1788581/1‑57 (MQ=255)
                                                                                          tAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTc                                                                                                                <  1:1788581/57‑1 (MQ=255)
                                                                                          tAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCg                              >  2:605340/1‑139 (MQ=255)
                                                                                            aGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGtt                            >  1:1748252/1‑139 (MQ=255)
                                                                                            aGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGtt                            >  1:674466/1‑139 (MQ=255)
                                                                                              cgcAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCa                                                                                               >  1:1574355/1‑70 (MQ=255)
                                                                                              cgcAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCa                                                                                               <  2:1574355/70‑1 (MQ=255)
                                                                                              cgcAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTAt                          >  2:583779/1‑139 (MQ=255)
                                                                                                cAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCa                                                                                                              <  2:325218/53‑1 (MQ=255)
                                                                                                cAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCa                                                                                                              >  1:325218/1‑53 (MQ=255)
                                                                                                 aaCGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAg                                                                                                                           >  1:1092929/1‑39 (MQ=255)
                                                                                                 aaCGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAg                                                                                                                           <  2:1092929/39‑1 (MQ=255)
                                                                                                  aCGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTc                                                                                                                >  2:77472/1‑49 (MQ=255)
                                                                                                  aCGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTc                                                                                                                <  1:77472/49‑1 (MQ=255)
                                                                                                  aCGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGc                      >  2:584102/1‑139 (MQ=255)
                                                                                                   cGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGAc                                                                                                   >  2:1845411/1‑61 (MQ=255)
                                                                                                   cGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGAc                                                                                                   <  1:1845411/61‑1 (MQ=255)
                                                                                                   cGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTa                                                                          >  2:388592/1‑86 (MQ=255)
                                                                                                   cGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTa                                                                          <  1:388592/86‑1 (MQ=255)
                                                                                                       gATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCgg                 <  1:1266582/139‑1 (MQ=255)
                                                                                                        aTGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGcc                                                                                   <  1:1415394/72‑1 (MQ=255)
                                                                                                        aTGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGcc                                                                                   >  2:1415394/1‑72 (MQ=255)
                                                                                                              tGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCAc                                                                               <  2:908217/70‑1 (MQ=255)
                                                                                                              tGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCAc                                                                               >  1:908217/1‑70 (MQ=255)
                                                                                                                      cGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCg  >  1:1296514/1‑139 (MQ=255)
                                                                                                                            |                                                                                                                                    
CTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCG  >  NC_000913/469305‑469561

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