Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 395,476 G→T 75.1% S271Y (TCT→TAT)  ampH ← D‑alanyl‑D‑alanine‑ carboxypeptidase/endopeptidase; penicillin‑binding protein; weak beta‑lactamase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913395,4760GT75.1% 158.9 / 63.6 96S271Y (TCT→TAT) ampHD‑alanyl‑D‑alanine‑ carboxypeptidase/endopeptidase; penicillin‑binding protein; weak beta‑lactamase
Reads supporting (aligned to +/- strand):  ref base G (12/12);  new base T (38/34);  total (50/46)
Fisher's exact test for biased strand distribution p-value = 8.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

CTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCATCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCAGATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCACGGACTGGCACCGCGCTCAGCAACCATTAAACGGCGGCACTGATCCGGTGAGGGGGT  >  NC_000913/395331‑395621
                                                                                                                                                 |                                                                                                                                                 
ctctTTCGGGGCCATATATACCCAACCTAAACCAAGCGCATCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCGGATAAAAATCATATg                                                                                                                                                <  3:22675/149‑1 (MQ=255)
 tctTTCGGGGCCATATATACCCAACCTAAACCAAGCGCATCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGa                                                                                                                                               >  8:336316/1‑149 (MQ=255)
 tctTTCGGGGCCATATATACCCAACCTAAACCAAGCGCATCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCAGATGa                                                                                                                                               >  2:258846/1‑149 (MQ=255)
     tCGGGGCCATATATACCCAACCTAAACCAAGCGCATCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGAc                                                                                                                                              <  6:353776/146‑1 (MQ=255)
             atatatACCCAACCTAAACCACGCGCATCCGCTTTGCCGGCGACATCCGTGCCCATCACTTTAGTAAACTGCGCGCGCTGATAAATCAACATCTGCATGCGATCTGACTCTTTACTGCGCTGATAACAATCATATTACAAATActgctg                                                                                                                                   >  8:96591/1‑149 (MQ=255)
                tataCCCAACCTAAACCAAGCGCATCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCAGATGACAAATACTGCTGCAt                                                                                                                                >  4:168971/1‑149 (MQ=255)
                  taCCCAACCTAAACCAAGCGCATCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATc                                                                                                                               <  1:361446/148‑1 (MQ=255)
                    cacAACCTAAACCAAGCGCATCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAg                                                                                                                            <  7:396359/147‑1 (MQ=255)
                     ccAACCTAAACCAAGCGCATCCGCTTTGCCGGGGACATCCATGCCAATCCCTTTCGTAAACTGCGCGCTCTGATAAATCAACGTCTGCATGTGATCTGCCTGGTTACTGCGCTGATAAAAATCAGATGACAAATACTGCTGCATCCAgc                                                                                                                           >  3:27473/1‑149 (MQ=255)
                     ccAACCTAAACCAAGCGCATCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTCCTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAgc                                                                                                                           <  1:127330/149‑1 (MQ=255)
                     ccAACCTAAACCAAGCGCATCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAgc                                                                                                                           >  7:421839/1‑149 (MQ=255)
                       aaCCTAAACCAAGCGCATCCGCTTTGCCGGGGACATCCGTGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAgcg                                                                                                                          >  2:185223/1‑148 (MQ=255)
                       aaCCTAAACCAAGCGCATCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAgcg                                                                                                                          >  8:357577/1‑148 (MQ=255)
                         ccTAAACCAAGCGCATCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGcat                                                                                                                       <  5:363407/149‑1 (MQ=255)
                         ccTAAACCAAGCGCATCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCAGa                                                                                                                                                  <  8:223113/122‑1 (MQ=255)
                         ccTAAACCAAGCGCATCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCAGa                                                                                                                                                  >  7:223113/1‑122 (MQ=255)
                         ccTAAACCAAGCGCATCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGAGAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGcat                                                                                                                       >  4:92607/1‑149 (MQ=255)
                           tAAACCAAGCGCATCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGcatca                                                                                                                     >  4:138584/1‑149 (MQ=255)
                             aaCCAAGCGCATCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCata                                                                                                                   <  3:138584/149‑1 (MQ=255)
                             aaCCAAGCGCATCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGAAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCAGATGACAAATACTGCTGCATCCAGCGcatca                                                                                                                     <  7:68493/147‑1 (MQ=255)
                              aCCAAGCGCATCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCAGATGACAAATACTGCTGCATCCAGCGCATCatat                                                                                                                  <  1:258846/149‑1 (MQ=255)
                               ccAAGCGCATCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCTCTGATAAAAATCATATGACAAATACTGCTGCATCAAGCGAATCATATc                                                                                                                 >  6:341042/1‑149 (MQ=255)
                                cAAGCGCATCCGCTTTGCCGGGGACATACATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCg                                                                                                                >  1:240794/1‑149 (MQ=255)
                                 aaGCGCATCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGc                                                                                                               >  3:180975/1‑149 (MQ=255)
                                  aGCGCATCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGc                                                                                                               >  4:91799/1‑148 (MQ=255)
                                   gcgcATCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCAGATGACAAATACTGCTGCATCCAGCGCATCATATCGCCa                                                                                                             <  8:75718/149‑1 (MQ=255)
                                       aTCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGGTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCCTCCAGCGCATCATATCGCCAGGCg                                                                                                         <  6:133647/149‑1 (MQ=255)
                                          cGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTg                                                                                                       <  4:180975/148‑1 (MQ=255)
                                          cGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTg                                                                                                       <  8:421839/148‑1 (MQ=255)
                                           gCTTTGCCGGTGACATCCCTGCAAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTg                                                                                                       <  4:279324/147‑1 (MQ=255)
                                           gCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTg                                                                                                       >  3:279324/1‑147 (MQ=255)
                                            cTTTGCCTTGGACATACATGCCAATCACTTTCGTAAACTGAGCGAGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGaa                                                                                                    <  5:341042/149‑1 (MQ=255)
                                            cTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGaa                                                                                                    >  5:341224/1‑149 (MQ=255)
                                               tGCCGGGGACATCCGTGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATa                                                                                                  <  1:185223/148‑1 (MQ=255)
                                               tGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATa                                                                                                  <  7:357577/148‑1 (MQ=255)
                                                gCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATGAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAg                                                                                                 <  2:429439/148‑1 (MQ=255)
                                                 ccGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGAc                                                                                               <  3:64185/149‑1 (MQ=255)
                                                          tCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGaa                                                                                                    >  6:162723/1‑135 (MQ=255)
                                                           ccATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTg                                                                                     >  3:95315/1‑149 (MQ=255)
                                                             aTGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCAGATGACAAATACTGCTGCATCCAGCGCATCATATCGCAAGGCGTGGAATAGACGCCGCCGCTGc                                                                                    >  1:428569/1‑148 (MQ=255)
                                                              tGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCAGATGACAAATACTGCTGCATCCAGCGcatca                                                                                                                     >  4:155620/1‑114 (MQ=255)
                                                              tGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCAGATGACAAATACTGCTGCATCCAGCGcatca                                                                                                                     <  3:155620/114‑1 (MQ=255)
                                                              tGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCAGATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGcc                                                                                   >  2:149919/1‑148 (MQ=255)
                                                                 cAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCCGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAAtt                                                                               >  1:316232/1‑149 (MQ=255)
                                                                        tttCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCa                                                                        >  1:152915/1‑149 (MQ=255)
                                                                             tAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCgtg                                                                   >  3:357574/1‑149 (MQ=255)
                                                                               aaCTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCgtgt                                                                  <  7:336316/148‑1 (MQ=255)
                                                                                  tGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCAGATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTAtt                                                              >  3:157531/1‑149 (MQ=255)
                                                                                    cgcgcgcTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTATTg                                                             <  5:382367/148‑1 (MQ=255)
                                                                                    cgcgcgcTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCAGATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTAtt                                                              <  8:162324/147‑1 (MQ=255)
                                                                                     gcgcgcTGATAAATCAACGGCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCAGATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCa                                                           <  3:168971/149‑1 (MQ=255)
                                                                                       gcgcTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCAc                                                          >  3:397675/1‑148 (MQ=255)
                                                                                         gcTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATGAGATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCACGGa                                                       <  1:149919/149‑1 (MQ=255)
                                                                                             aTAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCACGGACTg                                                    <  4:384221/148‑1 (MQ=255)
                                                                                                aaTCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCgtgt                                                                  <  8:168910/131‑1 (MQ=255)
                                                                                                aaTCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCgtgt                                                                  >  7:168910/1‑131 (MQ=255)
                                                                                                aaTCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCACGGACTGGCAc                                                <  4:95315/149‑1 (MQ=255)
                                                                                                  tCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGca                                                                                                                        <  6:411398/75‑1 (MQ=255)
                                                                                                  tCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGca                                                                                                                        >  5:411398/1‑75 (MQ=255)
                                                                                                  tCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGA‑‑‑CGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCACGGACTGGCACCgcgc                                           <  3:92607/149‑1 (MQ=255)
                                                                                                          tGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCAGATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCACGGACTGGCACCGCGCTCAg                                       >  3:5910/1‑148 (MQ=255)
                                                                                                            cATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCGTCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCACGGACTGGCACCGCGCTCAGCaa                                    >  8:71112/1‑149 (MQ=255)
                                                                                                            cATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCAGATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCACGGACTGGCACCGCGCTCAGCaa                                    >  5:20687/1‑149 (MQ=255)
                                                                                                              tGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCACGGACTGGCACCGCGCTCAGCAAc                                   >  8:23334/1‑148 (MQ=255)
                                                                                                              tGCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCACGGACTGGCACCGCGCTCAGCAAc                                   >  2:294115/1‑148 (MQ=255)
                                                                                                               gCGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCACGGACTGGCACCGCGCTCAGCAACCa                                 >  1:19093/1‑149 (MQ=255)
                                                                                                                cGATCTGCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGGGTTATTGCACGGACTGGCACCGCGCTCAGCAACCAt                                <  1:20878/149‑1 (MQ=255)
                                                                                                                cGATCTGCCTGGTTACTGCGCTGATAAAAATCAGATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCACGGACTGGCACCGCGCTCAGCAACCAt                                <  6:20687/149‑1 (MQ=255)
                                                                                                                  aTCTGCCTGGTTACTGCGCTGATAAAAATCAGATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCAGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCACGGACTGGCACCGCGCTCAGCAACCAtt                               <  4:27473/148‑1 (MQ=255)
                                                                                                                      gCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCACGGAc                                                      <  1:179996/121‑1 (MQ=255)
                                                                                                                      gCCTGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCACGGAc                                                      >  2:179996/1‑121 (MQ=255)
                                                                                                                         tGGTTACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCACGGACTGGCACCGCGCTCAGCAACCATTAAAcggcg                       >  2:113119/1‑149 (MQ=255)
                                                                                                                           gTTACTGCGCTGATAAAAATCAGATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCACGGACTGGCACCGCGCTCAGCAACCATTAAAcggcgg                      <  4:157531/148‑1 (MQ=255)
                                                                                                                            ttACTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCTTCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCACGGACTGGCACCGCGCTCAGCAACCATTAAACGGCGGCa                    <  7:71112/149‑1 (MQ=255)
                                                                                                                              aCTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCACGGACTGGCACCGCGCTCAGCAAc                                   <  6:128929/132‑1 (MQ=255)
                                                                                                                              aCTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCACGGACTGGCACCGCGCTCAGCAAc                                   >  5:128929/1‑132 (MQ=255)
                                                                                                                              aCTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCACGGACTGGCACCGCGCTCAGCAACCATTAAACGGCGGCAc                   <  2:152915/148‑1 (MQ=255)
                                                                                                                              aCTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCACGGACTGGCACCGCGCTCAGCAACCATTAAACGGCGGCAc                   <  6:341224/148‑1 (MQ=255)
                                                                                                                              aCTGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCACGGACTGGCACCGCGCTCAGCAACCATTAAACGGCGGCACt                  >  5:372767/1‑149 (MQ=255)
                                                                                                                                tGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCACGGACTGGCACCGCGCTCAGCAACCATTAAACGGCGGCACTGa                >  3:331965/1‑149 (MQ=255)
                                                                                                                                tGCGCTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCACGGACTGGCACCGCGCTCAGCAACCATTAAACGGCGGCACTGa                >  4:231869/1‑149 (MQ=255)
                                                                                                                                 gcgcTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCACGGACTGGCACCGCGCTCAGCAACCATTAAACGGCGGCACTGAt               >  4:116068/1‑149 (MQ=255)
                                                                                                                                   gcTGATAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCACGGACTGGCACCGCGCTCAGCAACCATTAAACGGCGGCACTGATcc             >  3:169338/1‑149 (MQ=255)
                                                                                                                                       aTAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCACGGACTGGCACCGCGCTCAGCAACCATTAAACGGCGGCACTGATCCGGt          <  7:194883/148‑1 (MQ=255)
                                                                                                                                       aTAAAAATCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCACGGACTGGCACCGCGCTCAGCAACCATTAAACGGCGGCACTGATCCGGt          <  5:97936/148‑1 (MQ=255)
                                                                                                                                       aTAAAAATCAGATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCACGGACTGGCACCGCGCTCAGCAACCATTAAACGGCGGCACTGAt               >  6:206346/1‑143 (MQ=255)
                                                                                                                                       aTAAAAATCAGATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCACGGACTGGCACCGCGCTCAGCAACCATTAAACGGCGGCACTGAt               <  5:206346/143‑1 (MQ=255)
                                                                                                                                         aaaaaTCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCgtg                                                                   <  4:193457/89‑1 (MQ=255)
                                                                                                                                         aaaaaTCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCgtg                                                                   >  3:193457/1‑89 (MQ=255)
                                                                                                                                         aaaaaTCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCACGGACTGGCACCGCGCTCAGCAACCATTAAACGGCGGCACTGATCCGGTGAg       >  2:403671/1‑149 (MQ=255)
                                                                                                                                         aaaaaTCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCACGGACTGGCACCGCGCTCAGCAACCATTAAACGGCGGCACTGATCCGGTGAg       >  8:98237/1‑149 (MQ=255)
                                                                                                                                           aaaTCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCACGGACTGGCACCGCGCTCAGCAACCATTAAACGGCGGCACTGATCCGGTGAggg     >  2:39572/1‑149 (MQ=255)
                                                                                                                                            aaTCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCACGGACTGGCACCGCGCTCAGCAACCATTAAACGGCGGCACTGATCCGGTGAgggg    <  4:357574/149‑1 (MQ=255)
                                                                                                                                            aaTCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCACGGACTGGCACCGCGCTCAGCAACCATTAAACGGCGGCACTGATCCGGTGAgggg    <  3:91799/149‑1 (MQ=255)
                                                                                                                                              tCATATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCACGGACTGGCACCGCGCTCAGCAACCATTAAACGGCGGCACTGATCCGGTGAGGGGgt  >  3:94973/1‑149 (MQ=255)
                                                                                                                                              tCAGATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCACGGACTGGCACCGCGCTCAGCAACCATTAAACGGCGGCACTGATCCGGTGAGGGGgt  >  3:374772/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCATCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTGCGCTGATAAAAATCAGATGACAAATACTGCTGCATCCAGCGCATCATATCGCCAGGCGTGGAATAGACGCCGCCGCTGCCAATTGCCGCCAGCGTGTTATTGCACGGACTGGCACCGCGCTCAGCAACCATTAAACGGCGGCACTGATCCGGTGAGGGGGT  >  NC_000913/395331‑395621

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