Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,468,399 G→A G259S (GGT→AGT)  mgtA → Mg(2(+)) importing P‑type ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,468,3990GA100.0% 1068.1 / NA 256G259S (GGT→AGT) mgtAMg(2(+)) importing P‑type ATPase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (127/129);  total (127/129)

Only 100 of 256 total aligned reads displayed.

GCCAGCCGGAGCACAGCAATCCGCTGGAGTGCGACACGCTGTGTTTTATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAGGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGAATGCCTTTCAGCAAGGGATCAGCCGCGTCAGTATGCT  >  NC_000913/4468301‑4468496
                                                                                                  |                                                                                                 
gCCAGCCGGAGCACAGCAATCCGCTGGAGTGCGACACGCTGTGTTTTATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGt                                                                                                 >  1:5910909/1‑101 (MQ=255)
  cAGCCGGAGCACAGCAATCCGCTGGAGTGCGACACGCTGTGTTTTATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGc                                                                                               >  1:5910913/1‑101 (MQ=255)
     ccGGAGCACAGCAATCCGCTGGAGTGCGACACGCTGTGTTTTATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCaa                                                                                            <  2:5910526/101‑1 (MQ=255)
     ccGGAGCACAGCAATCCGCTGGAGTGCGACACGCTGTGTTTTATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCaa                                                                                            >  2:5910915/1‑101 (MQ=255)
      cGGAGCACAGCAATCCGCTGGAGTGCGACACGCTGTGTTTTATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAAt                                                                                           >  2:5910918/1‑101 (MQ=255)
       ggAGCACAGCAATCCGCTGGAGTGCGACACGCTGTGTTTTATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATa                                                                                          <  1:5910590/101‑1 (MQ=255)
       ggAGCACAGCAATCCGCTGGAGTGCGACACGCTGTGTTTTATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATa                                                                                          >  1:5910921/1‑101 (MQ=255)
         aGCACAGCAATCCGCTGGAGTGCGACACGCTGTGTTTTATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATAcc                                                                                        <  1:5910514/101‑1 (MQ=255)
         aGCACAGCAATCCGCTGGAGTGCGACACGCTGTGTTTTATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATAcc                                                                                        <  1:5910651/101‑1 (MQ=255)
           cacaGCAATCCGCTGGAGTGCGACACGCTGTGTTTTATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTg                                                                                      <  1:5910622/101‑1 (MQ=255)
             caGCAATCCGCTGGAGTGCGACACGCTGTGTTTTATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGt                                                                                    >  2:5910925/1‑101 (MQ=255)
              aGCAATCCGCTGGAGTGCGACACGCTGTGTTTTATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGtt                                                                                   <  1:5910450/101‑1 (MQ=255)
              aGCAATCCGCTGGAGTGCGACACGCTGTGTTTTATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGtt                                                                                   <  1:5910504/101‑1 (MQ=255)
                cAATCCGCTGGAGTGCGACACGCTGTGTTTTATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTg                                                                                 >  1:5910927/1‑101 (MQ=255)
                 aaTCCGCTGGAGTGCGACACGCTGTGTTTTATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTgg                                                                                <  2:5910563/101‑1 (MQ=255)
                   tCCGCTGGAGTGCGACACGCTGTGTTTTATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTc                                                                              <  2:5910609/101‑1 (MQ=255)
                    ccGCTGGAGTGCGACACGCTGTGTTTTATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCa                                                                             >  1:5910930/1‑101 (MQ=255)
                    ccGCTGGAGTGCGACACGCTGTGTTTTATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCa                                                                             >  2:5910931/1‑101 (MQ=255)
                     cGCTGGAGTGCGACACGCTGTGTTTTATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCaa                                                                            <  1:5910655/101‑1 (MQ=255)
                      gCTGGAGTGCGACACGCTGTGTTTTATGGGCACCACCGTGGTGAGCGGTACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAAc                                                                           <  1:5910628/101‑1 (MQ=255)
                      gCTGGAGTGCGACACGCTGTGTTTTATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAAc                                                                           <  1:5910636/101‑1 (MQ=255)
                       cTGGAGTGCGACACGCTGTGTTTTATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACt                                                                          <  1:5910600/101‑1 (MQ=255)
                        tGGAGTGCGACACGCTGTGTTTTATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTg                                                                         >  1:5910938/1‑101 (MQ=255)
                        tGGAGTGCGACACGCTGTGTTTTATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTg                                                                         >  2:5910935/1‑101 (MQ=255)
                         ggAGTGCGACACGCTGTGTTTTATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTgg                                                                        >  2:5910941/1‑101 (MQ=255)
                          gAGTGCGACACGCTGTGTTTTATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGc                                                                       <  2:5910644/101‑1 (MQ=255)
                          gAGTGCGACACGCTGTGTTTTATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGc                                                                       <  2:5910709/101‑1 (MQ=255)
                          gAGTGCGACACGCTGTGTTTTATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGc                                                                       <  1:5910634/101‑1 (MQ=255)
                           aGTGCGACACGCTGTGTTTTATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCg                                                                      <  1:5910492/101‑1 (MQ=255)
                            gTGCGACACGCTGTGTTTTATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCgg                                                                     >  2:5910943/1‑101 (MQ=255)
                             tGCGACACGCTGTGTTTTATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGt                                                                               >  1:5910945/1‑90 (MQ=255)
                               cGACACGCTGTGTTTTATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGc                                                                  <  1:5910678/101‑1 (MQ=255)
                               cGACACGCTGTGTTTTATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGc                                                                  <  2:5910541/101‑1 (MQ=255)
                                 acacGCTGTGTTTTATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCaa                                                                            >  2:5910948/1‑89 (MQ=255)
                                  cacGCTGTGTTTTATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCgtg                                                               >  2:5910950/1‑101 (MQ=255)
                                      cTGTGTTTTATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAg                                                           <  1:5910339/101‑1 (MQ=255)
                                        gtgtTTTATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTg                                                         >  1:5910952/1‑101 (MQ=255)
                                          gtTTTATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAg                                                       <  1:5910830/101‑1 (MQ=255)
                                           ttttATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGc                                                      >  1:5910954/1‑101 (MQ=255)
                                           ttttATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGc                                                      <  2:5910775/101‑1 (MQ=255)
                                           ttttATGGGCACCACCGTGGTGAGCGGCACGGCACAACCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGc                                                      <  1:5910697/101‑1 (MQ=255)
                                            tttATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGc                                                                       >  1:5910959/1‑83 (MQ=255)
                                             ttATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCg                                                                 >  1:5910961/1‑88 (MQ=255)
                                              tATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCgtg                                                               >  2:5910963/1‑89 (MQ=255)
                                              tATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAgg                                                   <  2:5910488/101‑1 (MQ=255)
                                              tATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAgg                                                   <  1:5910731/101‑1 (MQ=255)
                                               aTGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGa                                                  >  1:5910965/1‑101 (MQ=255)
                                                 gggCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGaaa                                                >  1:5910967/1‑101 (MQ=255)
                                                 gggCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGaaa                                                <  1:5910690/101‑1 (MQ=255)
                                                 gggCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGaaa                                                <  2:5910400/101‑1 (MQ=255)
                                                    caccacCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCg                                             >  2:5910970/1‑101 (MQ=255)
                                                     accaccGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGa                                            >  1:5910972/1‑101 (MQ=255)
                                                        accGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGc                                              <  2:5910548/96‑1 (MQ=255)
                                                        accGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGcc                                         >  1:5910976/1‑101 (MQ=255)
                                                         ccGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCg                                        <  1:5910605/101‑1 (MQ=255)
                                                         ccGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCg                                        <  2:5910561/101‑1 (MQ=255)
                                                          cGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGa                                       <  1:5910804/101‑1 (MQ=255)
                                                          cGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGa                                       >  1:5910979/1‑101 (MQ=255)
                                                           gtggtgAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGaa                                      <  1:5910577/101‑1 (MQ=255)
                                                           gtggtgAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGaa                                      >  2:5910984/1‑101 (MQ=255)
                                                           gtggtgAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGaa                                      >  1:5910982/1‑101 (MQ=255)
                                                            tggtgAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGc                                                                  >  2:5910988/1‑72 (MQ=255)
                                                            tggtgAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGAAt                                     >  1:5910985/1‑101 (MQ=255)
                                                             ggtgAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGAATg                                    <  2:5910617/101‑1 (MQ=255)
                                                              gtgAGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGAATGc                                   <  2:5910641/101‑1 (MQ=255)
                                                                 aGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGAATGCCtt                                >  1:5910991/1‑101 (MQ=255)
                                                                 aGCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGAATGCCtt                                <  1:5910421/101‑1 (MQ=255)
                                                                  gCGGCACGGCACAAGCAATGGTGATTGCTACAAGTGCCAAt                                                                                           <  2:5910602/41‑1 (MQ=255)
                                                                     gcacggcacAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGAATGCCTTTCAg                            >  2:5910994/1‑101 (MQ=255)
                                                                      cacggcacAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGaaa                                                >  2:5910999/1‑80 (MQ=255)
                                                                      cacggcacAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGAATGCCTTTCAGc                           >  1:5910996/1‑101 (MQ=255)
                                                                       acggcacAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGAATGCCTTTCAGCa                          <  2:5910445/101‑1 (MQ=255)
                                                                        cggcacAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGAATGCCTTTCAGCaa                         >  2:5911002/1‑101 (MQ=255)
                                                                          gcacAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGAATGCCTTTCAGCAAgg                       >  1:5911003/1‑101 (MQ=255)
                                                                          gcacAAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGAATGCCTTTCAGCAAgg                       <  1:5910464/101‑1 (MQ=255)
                                                                           cacaAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGAATGCCTTTCAGCAAggg                      >  2:5911006/1‑101 (MQ=255)
                                                                             caAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGAATGCCTTTCAGCAAGGGAt                    <  1:5910700/101‑1 (MQ=255)
                                                                             caAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGAATGCCTTTCAGCAAGGGAt                    <  1:5910555/101‑1 (MQ=255)
                                                                             caAGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGAATGCCTTTCAGCAAGGGAt                    <  2:5910666/101‑1 (MQ=255)
                                                                              aaGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGAATGCCTTTCAGCAAGGGATc                   <  2:5910799/101‑1 (MQ=255)
                                                                              aaGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGAATGCCTTTCAGCAAGGGATc                   <  1:5910531/101‑1 (MQ=255)
                                                                               aGCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGAATGCCTTTCAGCAAGGGATCa                  <  2:5910687/101‑1 (MQ=255)
                                                                                gCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGAATGCCTTTCAGCAAGGGATCAg                 >  2:5911009/1‑101 (MQ=255)
                                                                                gCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGAATGCCTTTCAGCAAGGGATCAg                 >  1:5911011/1‑101 (MQ=255)
                                                                                gCAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGAATGCCTTTCAGCAAGGGATCAg                 <  1:5910672/101‑1 (MQ=255)
                                                                                 cAATGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGAATGCCTTTCAGCAAGGGATCAGc                >  1:5911013/1‑101 (MQ=255)
                                                                                  aaTGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGAATGCCTTTCAGCAAg                        >  1:5911018/1‑92 (MQ=255)
                                                                                   aTGGTTATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGAATGCCTTTCAGCAAGGGATCAGCcg              >  1:5911020/1‑101 (MQ=255)
                                                                                    tGGTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGAATGCCTTTCAGCAAGGGATCAGCcgc             <  2:5910596/101‑1 (MQ=255)
                                                                                     ggTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGAATGCCTTTCAGCAAGGGATCAGCcgcg            >  1:5911023/1‑101 (MQ=255)
                                                                                     ggTGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGAATGCCTTTCAGCAAGGGATCAGCcgcg            <  1:5910747/101‑1 (MQ=255)
                                                                                       tGATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGc                                                      >  1:5911025/1‑57 (MQ=255)
                                                                                        gATTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGAATGCCTTTCAGCAAGGGATCAGCCGCGTCa         >  2:5911026/1‑101 (MQ=255)
                                                                                         aTTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGAATGCCTTTCAGCAAGGGATCAGCCGCGTCAg        <  2:5910575/101‑1 (MQ=255)
                                                                                         aTTGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGAATGCCTTTCAGCAAGGGATCAGCCGCGTCAg        >  2:5911028/1‑101 (MQ=255)
                                                                                          ttGCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGAATGCCTTTCAGCAAGGGATCAGCCGCGTCAGt       <  2:5910752/101‑1 (MQ=255)
                                                                                            gCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGAATGCCTTTCAGCAAGGGATCAGCCGCGTCAGTAt     >  1:5911032/1‑101 (MQ=255)
                                                                                            gCTACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGAATGCCTTTCAGCAAGGGATCAGCCGCGTCAGTAt     <  1:5910760/101‑1 (MQ=255)
                                                                                             cTACAAGTGCCAATACCCGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGAATGCCTTTCAGCAAGGGATCAGCCGCGTCAGTAtg    >  2:5911034/1‑101 (MQ=255)
                                                                                              tACAAGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGAATGCCTTTCAGCAAGGGATCAGCCGCGTCAGTAtgc   <  2:5910764/101‑1 (MQ=255)
                                                                                                  |                                                                                                 
GCCAGCCGGAGCACAGCAATCCGCTGGAGTGCGACACGCTGTGTTTTATGGGCACCACCGTGGTGAGCGGCACGGCACAAGCAATGGTGATTGCTACAGGTGCCAATACCTGGTTTGGTCAACTGGCGGGGCGTGTTAGTGAGCAGGAAAGCGAGCCGAATGCCTTTCAGCAAGGGATCAGCCGCGTCAGTATGCT  >  NC_000913/4468301‑4468496

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