Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,523,857 G→T 5.6% S329S (TCG→TCT mrcA → penicillin‑binding protein 1a, murein transglycosylase and transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,523,8570GT 5.6% 239.1 / 3.5 72S329S (TCG→TCTmrcApenicillin‑binding protein 1a, murein transglycosylase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base G (33/35);  new base T (2/2);  total (35/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.21e-01

CTATGAAGACGGTTATCGCATTTACACCACCATCACCCGCAAAGTGCAGCAGGCCGCGCAGCAGGCGGTACGTAATAACGTGCTGGACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGC  >  NC_000913/3523710‑3524004
                                                                                                                                                   |                                                                                                                                                   
cTATGAAGACGGTTATCGCATTTACACCACCATCACCCGCAAAGTGCAGCAGGCCGCGCAGCAGGCGGTACGTAATAACGTGCTGGACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCgg                                                                                                                                                    >  6:58701/1‑149 (MQ=255)
cTATGAAGACGGTTATCGCATTTACACCACCAACCCCCGCAAAGTGCATCAGGCCGCGCAGCAGGCGGTACGTAATAACGAGCAGGACTACGAAATCCGCCACGGCTATCTCGGCCCCGCAAATGAGCTCTAGAAAGTGGCCCAATCgg                                                                                                                                                    >  5:201723/1‑149 (MQ=255)
     aaGACGGTTATCGCATTTACACCACCATCACCCGCAAAGTGCAGCAGGCCGCGCAGCAGGCGGTACGTAATAACGTGCTGGACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTgg                                                                                                                                               <  2:318167/149‑1 (MQ=255)
     aaGACGGTTATCGCATTTACACCACCATCACCCGCAAAGTGCAGCAGGCCGCGCAGCAGGCGGTACGTAATAACGTGCTGGACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCCAGTCGGCGTgg                                                                                                                                               >  5:107949/1‑149 (MQ=255)
               tCGCATTTACACCACCATCACCCGCAAAGTGCAGCAGGCCGCGCAGCAGGCGGTACGTAATAACGTGCTGGACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATaacaac                                                                                                                                      >  1:392196/1‑148 (MQ=255)
               tCGCATTTACACCACCATCACCCGCAAAGTGCAGCAGGCCGCGCAGCAGGCGGTACGTAATAACGTGCTGGACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATaacaac                                                                                                                                      <  6:215464/148‑1 (MQ=255)
                 gCATTTACACCACCATCACCCGCAAAGTGCAGCAGGCCGCGCAGCAGGCGGTACGTAATAACGTGCTGGACTAGGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAg                                                                                                                                   <  2:392196/149‑1 (MQ=255)
                  cATTTACACCACCATCACCCGCAAAGTGCAGCAGGCCGCTCAGCAGGCGGTACGTAATAACGTGCTGGACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATatcaa                                                                                                                                       <  3:88058/144‑1 (MQ=255)
                  cATTTACACCACCATCACCCGCAAAGTGAAGCAGGCCGCGCAGCAGGCGGTACGTAATAACGTGCTGGACTACGACATGCGCCACGGCTATCGCGGACAGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATaacaa                                                                                                                                       >  4:88058/1‑144 (MQ=255)
                        caccacCATCACCCGCAAAGTGCAGCAGGCCGCGCAGCAGGCGGTACGTAATAACGTGCTGGACTACGACATGCGCCAGGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTAcc                                                                                                                             <  5:58701/148‑1 (MQ=255)
                        caccacCATCACCCGCAAAGTGCAGCAGGCCGCGCAGCAGGCGGTACGTAATAACGTGCTGGACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTAcc                                                                                                                             <  7:114296/148‑1 (MQ=255)
                         accaccATCACCCGCAAAGTGCAGCAGGCCGCGCAGCAGGCGGTACGTAATAACGTGCTGGACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGa                                                                                                                           <  8:249030/149‑1 (MQ=255)
                           caccaTCACCCGCAAAGTGCAGCAGGCCGCGCAGCAGGCGGTCCGTAATAACGTGCTGGACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATa                                                                                                                         >  2:181825/1‑149 (MQ=255)
                           caccaTCACCCGCAAAGTGCAGCAGGCCGCGCAGCAGGCGGTACGTAATAACGTGCTGGACTACGACATGCTCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTAc                                                                                                                              >  6:173723/1‑144 (MQ=255)
                           caccaTCACCCGCAAAGTGCAGCAGGCCGCGCAGCAGGCGGTACGTAATAACGTGCTGGACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTAc                                                                                                                              <  5:173723/144‑1 (MQ=255)
                           caccaTCACCCGCAAAGTGCAGCAGGCCGCGCAGCAGGCGGTACGTAATAACGTGCTGGACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGa                                                                                                                           <  7:62959/147‑1 (MQ=255)
                           caccaTCACCCGCAAAGTGCAGCAGGCCGCGCAGCAGGCGGTACGTAATAACGTGCTGGACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGa                                                                                                                           >  8:62959/1‑147 (MQ=255)
                              caTCACCCGCAAAGTGCAGCCGGCCGCGCAGCAGGCGGTACGTAATAACGTGCTGGACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACCACCAGATTACCGATACGc                                                                                                                      >  7:351974/1‑149 (MQ=255)
                                tCACCCGCAAAGTGCAGCAGGCCGCGCAGCAGGCGGTACGTAATAACGTGCTGGACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCt                                                                                                                     >  7:190327/1‑148 (MQ=255)
                                  aCCCGCACACTGCCGCAGGCCGCGCAGCAGGCGGTACGTAATAACGTGCTGGACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGaa                                                                                                                  <  8:190327/149‑1 (MQ=255)
                                  aCCAGCAAAGTGCATGAGGCCGCGCAGTAGGCGGTACGTAATAACGTGATTGACTTCTAAATGCGCCAGGGCTATCGCGGCCCGGCAATTGTGTTGTGGAATGTGGGCGTGTCGGCGTGGGATTACTTGATGATTACCGATACGCTGa                                                                                                                   <  6:107949/148‑1 (MQ=255)
                                   cccGCAAAGTGCAGCAGGCCGCGCAGCAGGCGGTACGTAATAACGTGCTGGACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAg                                                                                                                 <  6:189070/149‑1 (MQ=255)
                                   cccGCAAAGTGCAGCAGGCCGCGCAGCAGGCGGTACGTAATAACGTGCTGGACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAg                                                                                                                 >  5:189070/1‑149 (MQ=255)
                                      gCAAAGTGCAGCAGGCCGCGCAGCAGGCGGTACGTAATAACGTGCTGGACTACGACATGCGCCCCGGCTATCGCCGCCCGGCAAATGGGCGGTGGAAAGTGGGCGAGGCGGGGTGGGTGAAGAACAAGATTACGGGTAGGCTGCAGgcg                                                                                                              >  2:392250/1‑149 (MQ=255)
                                       cAAAGTGCAGCAGGCCGCGCAGCAGGCGGTACGTAATAACGTGCTTGACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTTGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGgcgc                                                                                                             <  3:279896/149‑1 (MQ=255)
                                        aaaGTGCAGCAGGCCGCGCAGCAGGCGGTACGTAATAACGTGCTGGACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCt                                                                                                            <  8:386054/149‑1 (MQ=255)
                                          aGTGCAGCAGGCCGCGCAGCAGGCGGTACGTAATAACGTGCTGGACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGc                                                                                                          <  8:190501/149‑1 (MQ=255)
                                              tagcagGCCTCGCAGCAGGCTGTTCGTACTAACGTGCTGGACTACGACATTCGCCACGGCTATCGCGGCACTGCAAATGGGCTGTGGAAAGTTGGCGAGGCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAAc                                                                                                      <  8:380899/148‑1 (MQ=255)
                                              cagcagGCCGCGCAGCAGGCGGTACGTAATAACGTGCTGGACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAAc                                                                                                      <  8:135752/149‑1 (MQ=255)
                                               agcagGCCGCGCAGCATGCGGTACGTAATAACGTGCTGGACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACACGATTACCGATACGCTGa                                                                                                                   <  5:116430/135‑1 (MQ=255)
                                               agcagGCCGCGCAGCAGGCGGTACGTAATAACGTGCTGGACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATaacaac                                                                                                                                      <  1:251193/116‑1 (MQ=255)
                                               agcagGCCGCGCAGCAGGCGGTACGTAATAACGTGCTGGACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATaacaac                                                                                                                                      >  2:251193/1‑116 (MQ=255)
                                               agcagGCCGCGCAGCAGGCGGTACGTAATAACGTGCTGGACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGa                                                                                                                   >  6:116430/1‑135 (MQ=255)
                                                      cgcgCAGCAGGCGGTACGTAATAACGTGCTGGACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGCGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTc                                                                                              >  5:374211/1‑149 (MQ=255)
                                                                 cGGTACGTAATAACGTGCTGGACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGTAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGctgcct                                                                                   <  5:94553/149‑1 (MQ=255)
                                                                          ataaCGTGCTGGACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGACGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGt                                                                           <  1:181825/148‑1 (MQ=255)
                                                                               gTGCTGGACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAg                                                                     >  3:16098/1‑149 (MQ=255)
                                                                                 gCTGGACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAgcg                                                                   >  5:81592/1‑149 (MQ=255)
                                                                                 gCTGGACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAgcg                                                                   >  2:270164/1‑149 (MQ=255)
                                                                                     gACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCaa                                                               >  7:400571/1‑149 (MQ=255)
                                                                                     gACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCa                                                                >  4:164001/1‑148 (MQ=255)
                                                                                     gACTACGACATGCCCCACGGCTATCGAGGCCATGAAAAAGTGAAGAGGAAAGTTGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTCCCAACCTATCGTCCGCTGCTTCCTGCCGCAGACACAAgctaaa                                                                >  6:49653/1‑144 (MQ=255)
                                                                                        tACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCTGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCg                                                                                            <  7:329072/117‑1 (MQ=255)
                                                                                        tACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCTGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCg                                                                                            <  7:323374/117‑1 (MQ=255)
                                                                                        tACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCTGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCg                                                                                            >  8:329072/1‑117 (MQ=255)
                                                                                        tACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCTGCGAGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCg                                                                                            >  8:323374/1‑117 (MQ=255)
                                                                                          cGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGAt                                                                                                                                            <  4:198911/67‑1 (MQ=255)
                                                                                          cGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGAt                                                                                                                                            >  3:198911/1‑67 (MQ=255)
                                                                                           gACATGCGCCCCGGCTATCGCGGCCCGGCAAATGGGCTGGGGAACGGGGGCGAGTCGGCGGGGGAAAACAACAAGATAACCGATACGCGGAAGGCGCTGCCAACCTATGGTCCGCTGCGGCCTGCCGCAGTCACCAGCGCCAAGCCTCa                                                         <  6:374211/149‑1 (MQ=255)
                                                                                              aTGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGc                                                       >  8:55765/1‑148 (MQ=255)
                                                                                              aTGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGACTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGc                                                       >  2:363673/1‑148 (MQ=255)
                                                                                                gcgcCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAg                                                    >  1:183513/1‑149 (MQ=255)
                                                                                                  gcCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCa                                                                                                        <  6:167536/95‑1 (MQ=255)
                                                                                                  gcCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCa                                                                                                        >  5:167536/1‑95 (MQ=255)
                                                                                                   ccACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCCCCAGCGCCAATCCTCAGCAAGCGa                                                 >  8:318087/1‑149 (MQ=255)
                                                                                                      cGGCTATCGCGGCCCGGCAAATGTGCTGGGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACgg                                              <  3:298735/149‑1 (MQ=255)
                                                                                                           aTCGCGGCCCGGCAAAGGTGCTGTGGAAAGTGGGCGAGGCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCGGCCGCCCTATGGTCCGCTGCTGCCTGCCGCAGGCACCAGCGCCAATCCTCCGCAAGCGCCGGCGagg                                         >  6:100140/1‑147 (MQ=255)
                                                                                                              gcgGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTg                                      <  1:270164/149‑1 (MQ=255)
                                                                                                              gcgGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTg                                      >  7:128610/1‑149 (MQ=255)
                                                                                                              gcgGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTg                                      >  5:40908/1‑149 (MQ=255)
                                                                                                                 gCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCg                                   >  4:125000/1‑149 (MQ=255)
                                                                                                                    cGGCAAATGTGCTGTGGAAAGTGGGCGAGGCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGAc                                <  4:380204/149‑1 (MQ=255)
                                                                                                                      gCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACgg                              <  6:40908/149‑1 (MQ=255)
                                                                                                                      gCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGACTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACgg                              <  1:363673/149‑1 (MQ=255)
                                                                                                                        aaaTGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGAc                                <  5:294654/145‑1 (MQ=255)
                                                                                                                        aaaTGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGAc                                >  6:294654/1‑145 (MQ=255)
                                                                                                                        aaaTGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGt                            <  3:164001/149‑1 (MQ=255)
                                                                                                                        aaaTGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGt                            >  2:310483/1‑149 (MQ=255)
                                                                                                                                  tgGAAAGTGCGCGAGTCGGCGTGGGATAACAACACGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTc                    <  2:218027/147‑1 (MQ=255)
                                                                                                                                        gTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAg            <  1:310483/149‑1 (MQ=255)
                                                                                                                                           ggCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTAt         >  7:132372/1‑149 (MQ=255)
                                                                                                                                            gCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATg        >  4:215963/1‑149 (MQ=255)
                                                                                                                                              gAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGa      <  3:215963/149‑1 (MQ=255)
                                                                                                                                                  cGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGc  <  8:400571/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CTATGAAGACGGTTATCGCATTTACACCACCATCACCCGCAAAGTGCAGCAGGCCGCGCAGCAGGCGGTACGTAATAACGTGCTGGACTACGACATGCGCCACGGCTATCGCGGCCCGGCAAATGTGCTGTGGAAAGTGGGCGAGTCGGCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGC  >  NC_000913/3523710‑3524004

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