Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,557,347 A→G 7.4% L231P (CTG→CCG)  ddpF ← D,D‑dipeptide ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,557,3470AG 7.4% 306.0 / 11.1 94L231P (CTG→CCG) ddpFD,D‑dipeptide ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (44/43);  new base G (3/4);  total (47/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

CAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAATTGAG  >  NC_000913/1557200‑1557495
                                                                                                                                                   |                                                                                                                                                    
cAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAg                                                                                                                                                     <  4:357027/149‑1 (MQ=255)
 aaCGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGc                                                                                                                                                    <  3:33263/149‑1 (MQ=255)
    gTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCTGACGGCAc                                                                                                                                                  <  4:362518/148‑1 (MQ=255)
            aaaaaaCAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTgcgc                                                                                                                                         >  4:1345823/1‑149 (MQ=255)
             aaaaaCAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGcgcg                                                                                                                                        <  3:1345823/149‑1 (MQ=255)
                 aCAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCTGACGGCACctgc                                                                                                                                              >  4:858689/1‑139 (MQ=255)
                 aCAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCTGACGGCACctgc                                                                                                                                              <  3:858689/139‑1 (MQ=255)
                           ggCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCtgctgc                                                                                                                                           <  2:951437/132‑1 (MQ=255)
                           ggCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCtgctgc                                                                                                                                           >  1:951437/1‑132 (MQ=255)
                             cAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTa                                                                                                                           >  3:1361398/1‑146 (MQ=255)
                             cAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTa                                                                                                                           <  4:1361398/146‑1 (MQ=255)
                               aCGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATc                                                                                                                      >  1:614920/1‑149 (MQ=255)
                                   gCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGcc                                                                                                                  >  1:282983/1‑149 (MQ=255)
                                             ggCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTa                                                                                                                           <  1:960829/130‑1 (MQ=255)
                                             ggCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATAGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCACTTCTa                                                                                                                           >  2:960829/1‑130 (MQ=255)
                                                  aTCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGccc                                                                                                                 <  1:944106/135‑1 (MQ=255)
                                                  aTCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGccc                                                                                                                 >  2:944106/1‑135 (MQ=255)
                                                  aTCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCg                                                                                                   >  1:478951/1‑149 (MQ=255)
                                                   tCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGc                                                                                                  >  1:783342/1‑149 (MQ=255)
                                                           aCGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGAt                                                                                          <  2:614920/149‑1 (MQ=255)
                                                                  gCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCt                                                                                      <  1:797511/146‑1 (MQ=255)
                                                                  gCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCt                                                                                      >  2:797511/1‑146 (MQ=255)
                                                                   cccATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAAACCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTc                                                                                                                             <  2:219242/106‑1 (MQ=255)
                                                                   cccATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAAACCAGCAATAATCGGGGGTATGGATGTTCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTc                                                                                                                             >  1:219242/1‑106 (MQ=255)
                                                                    ccATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCTGACGGCACctg                                                                                                                                               >  1:153216/1‑87 (MQ=255)
                                                                    ccATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCTGACGGCACctg                                                                                                                                               <  2:153216/87‑1 (MQ=255)
                                                                    ccATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATc                                                                                         <  1:51684/141‑1 (MQ=255)
                                                                    ccATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATc                                                                                         >  2:51684/1‑141 (MQ=255)
                                                                    ccATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATc                                                                                         <  3:871143/141‑1 (MQ=255)
                                                                    ccATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCatat                                                                                 >  2:1058436/1‑149 (MQ=255)
                                                                    ccATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAACACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATc                                                                                         >  4:871143/1‑141 (MQ=255)
                                                                     cATTCTTCGTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCAcc                                                                                                    >  4:122539/1‑129 (MQ=255)
                                                                     cATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCat                                                                                   <  2:1274247/146‑1 (MQ=255)
                                                                     cATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCat                                                                                   >  1:1274247/1‑146 (MQ=255)
                                                                     cATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGCCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCAcc                                                                                                    <  3:122539/129‑1 (MQ=255)
                                                                      aTTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCgaga                                                                                                             >  4:683057/1‑119 (MQ=255)
                                                                      aTTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCgaga                                                                                                             <  3:683057/119‑1 (MQ=255)
                                                                       ttcttcCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAAGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACAACACCGCCACCCGATCGCTCATATGAc                                                                              >  2:429955/1‑149 (MQ=255)
                                                                           tcctccAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTAt                                                                          >  2:1241815/1‑149 (MQ=255)
                                                                            cctccAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCAcc                                                                                                    <  3:839929/122‑1 (MQ=255)
                                                                            cctccAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCAcc                                                                                                    >  4:839929/1‑122 (MQ=255)
                                                                              tccAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCAc                                                                       >  1:636602/1‑149 (MQ=255)
                                                                                cAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACAt                                                                                                        <  3:405268/114‑1 (MQ=255)
                                                                                cAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACAt                                                                                                        >  4:405268/1‑114 (MQ=255)
                                                                                      tttATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACAt                                                               >  3:173275/1‑149 (MQ=255)
                                                                                       ttATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACAtt                                                              <  2:783342/149‑1 (MQ=255)
                                                                                        tATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATtg                                                             >  3:43531/1‑149 (MQ=255)
                                                                                         aTCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCga                                                                    >  4:1332469/1‑141 (MQ=255)
                                                                                         aTCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCga                                                                    <  3:1332469/141‑1 (MQ=255)
                                                                                          tCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGccc                                                                                                                 >  1:144879/1‑95 (MQ=255)
                                                                                          tCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGccc                                                                                                                 <  2:144879/95‑1 (MQ=255)
                                                                                           cAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATCTGACGTATCACCGAGACATagtgt                                                          <  1:1241815/149‑1 (MQ=255)
                                                                                           cAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATtgtgt                                                          <  1:1163976/149‑1 (MQ=255)
                                                                                           cAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATtgtgt                                                          >  2:1163976/1‑149 (MQ=255)
                                                                                                  ggggAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCa                                                   <  2:282983/149‑1 (MQ=255)
                                                                                                   gggagggaATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAg                                                  <  2:636602/149‑1 (MQ=255)
                                                                                                    ggagggaATCCAGCTATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGc                                                 >  3:315146/1‑149 (MQ=255)
                                                                                                    ggagggaATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCga                                                                    >  3:1051420/1‑130 (MQ=255)
                                                                                                    ggagggaATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCga                                                                    <  4:1051420/130‑1 (MQ=255)
                                                                                                     gagggaATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGca                                                <  3:598518/149‑1 (MQ=255)
                                                                                                        ggaATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTc                                                                                                                                      <  2:594224/60‑1 (MQ=255)
                                                                                                        ggaATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTc                                                                                                                                      >  1:594224/1‑60 (MQ=255)
                                                                                                         gaATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCAc                                                                                               >  2:1308369/1‑98 (MQ=255)
                                                                                                         gaATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCAc                                                                                               <  1:1308369/98‑1 (MQ=255)
                                                                                                             ccAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTc                                        >  1:316653/1‑149 (MQ=255)
                                                                                                              cAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATtgtgt                                                          <  3:128546/130‑1 (MQ=255)
                                                                                                              cAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGATCATATGACGTATCACCGAGACATtgtgt                                                          >  4:128546/1‑130 (MQ=255)
                                                                                                               aGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATa                                                                                                           <  1:475277/80‑1 (MQ=255)
                                                                                                               aGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATa                                                                                                           >  2:475277/1‑80 (MQ=255)
                                                                                                                gCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCa                                                                                                <  2:458420/90‑1 (MQ=255)
                                                                                                                gCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCa                                                                                                >  1:458420/1‑90 (MQ=255)
                                                                                                                gCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGc                                     >  3:595971/1‑149 (MQ=255)
                                                                                                                      atCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGa                               <  2:316653/149‑1 (MQ=255)
                                                                                                                           gtgtATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTc                          <  1:1058436/149‑1 (MQ=255)
                                                                                                                             gtATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCt                        <  4:315146/149‑1 (MQ=255)
                                                                                                                              tATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCa                                                                                    <  3:411297/88‑1 (MQ=255)
                                                                                                                              tATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCa                                                                                    >  4:411297/1‑88 (MQ=255)
                                                                                                                              tATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATtgt                                                            <  1:387151/112‑1 (MQ=255)
                                                                                                                              tATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATtgt                                                            >  2:387151/1‑112 (MQ=255)
                                                                                                                               aTGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGa                               >  3:604434/1‑140 (MQ=255)
                                                                                                                               aTGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGa                               <  4:604434/140‑1 (MQ=255)
                                                                                                                               aTGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGt                      <  1:429955/149‑1 (MQ=255)
                                                                                                                                 ggATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCga                                                                    <  1:552189/101‑1 (MQ=255)
                                                                                                                                 ggATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCga                                                                    >  2:552189/1‑101 (MQ=255)
                                                                                                                                   aTGTGCAGGTGCCTTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCg                  >  2:1293024/1‑149 (MQ=255)
                                                                                                                                       gCAGGTGCTGACGGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGAc                                                                 >  4:1392233/1‑98 (MQ=255)
                                                                                                                                       gCAGGTGCTGACGGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGAc                                                                 <  3:1392233/98‑1 (MQ=39)
                                                                                                                                       gCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTAc              <  2:478951/149‑1 (MQ=255)
                                                                                                                                           gTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTa                                                                           >  1:83631/1‑84 (MQ=255)
                                                                                                                                           gTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTa                                                                           <  2:83631/84‑1 (MQ=255)
                                                                                                                                           gTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGt          >  4:311266/1‑149 (MQ=255)
                                                                                                                                                gTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAAtt     <  4:173275/149‑1 (MQ=255)
                                                                                                                                                gTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAAtt     <  4:595971/149‑1 (MQ=255)
                                                                                                                                                gTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAAtt     >  1:324036/1‑149 (MQ=255)
                                                                                                                                                   aGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAATTGAg  <  3:311266/149‑1 (MQ=255)
                                                                                                                                                   aGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAATTGAg  <  2:324036/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAATTGAG  >  NC_000913/1557200‑1557495

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