Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,827,055 T→G 8.5% S1379R (AGT→CGT) ‡ pvdL ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,827,0550TG 8.5% 214.5 / 10.2 71S1379R (AGT→CGT) ‡pvdLnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base T (30/35);  new base G (4/2);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.88e-01

TGCCAGTGGGTGACGTTGATGGCGGTTTCAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTGCAGCCGGTGACCA  >  NC_002947/4826907‑4827201
                                                                                                                                                    |                                                                                                                                                  
tGCCAGTGGGTGACGTTGATGGCGGTTTCAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACt                                                                                                                                                    >  6:357589/1‑149 (MQ=255)
  ccAGTGGGTGACGTTGATGGCGGTTTCAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACtt                                                                                                                                                   <  5:134629/148‑1 (MQ=255)
  ccAGTGGGTGACGTTGATGGCGGTTTCAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTg                                                                                                                                                  >  6:422950/1‑149 (MQ=255)
      tGGGTGACGTTGATGGCGGTTTCAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTg                                                                                                                                              <  8:10932/149‑1 (MQ=255)
         gTGACGTTGATGGCGGTTTCAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCgg                                                                                                                                           <  2:236614/149‑1 (MQ=255)
         gTGACGTTGATGGCGGTTTCAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCAAGTGCCTGcag                                                                                                                                           <  7:5184/149‑3 (MQ=255)
            aCGTTGATGGCGGTTTCAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTc                                                                                                                                        <  2:320590/149‑1 (MQ=255)
              gTTGATGGCGGTTTCAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCt                                                                                                                                      <  5:359193/149‑1 (MQ=255)
                tGATGGCGGTTTCAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTg                                                                                                                                     >  1:362798/1‑148 (MQ=255)
                 gATGGCGGTTTCAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGa                                                                                                                                   >  1:175781/1‑149 (MQ=255)
                   tGGCGGTTTCAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGAc                                                                                                                                  >  8:24163/1‑148 (MQ=255)
                      cGGTTTCAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGaa                                                                                                                               >  1:11090/1‑148 (MQ=255)
                            cAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTg                                                                                                                         >  4:183197/1‑148 (MQ=255)
                              gTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACAcc                                                                                                                           <  2:99899/144‑1 (MQ=255)
                              gTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACAcc                                                                                                                           >  1:99899/1‑144 (MQ=255)
                              gTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCaa                                                                                                                      <  3:210258/149‑1 (MQ=255)
                               tCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCaa                                                                                                                      >  2:123225/1‑148 (MQ=255)
                               tCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCaa                                                                                                                      >  6:224319/1‑148 (MQ=255)
                               tCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAAc                                                                                                                     <  4:306044/149‑1 (MQ=255)
                                 gggCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAg                                                                                                                   >  2:288756/1‑149 (MQ=255)
                                 gggCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAg                                                                                                                   >  6:157601/1‑149 (MQ=255)
                                  ggCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAgc                                                                                                                  <  7:36760/149‑1 (MQ=255)
                                     cATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAgcggc                                                                                                               <  2:228591/149‑1 (MQ=255)
                                     cATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAgcggc                                                                                                               <  8:238390/149‑1 (MQ=255)
                                       tACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGc                                                                                                                                             >  7:130847/1‑117 (MQ=255)
                                       tACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGc                                                                                                                                             <  8:130847/117‑1 (MQ=255)
                                          cGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTAc                                                                                                          <  5:224319/149‑1 (MQ=255)
                                           ggTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACa                                                                                                         <  1:123225/149‑1 (MQ=255)
                                             ttGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACaaa                                                                                                       <  5:357589/149‑1 (MQ=255)
                                                 gCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTg                                                                                                    >  3:360577/1‑148 (MQ=255)
                                                    gTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTTCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGc                                                                                                <  3:208374/149‑1 (MQ=255)
                                                     ttGCACCTGCGGCTGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGTCAACGCTTCGCCTCCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGATGAACACCTGCAACAGCGGCGGTACAAAGTGCAGc                                                                                                <  5:101599/148‑1 (MQ=255)
                                                        cACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTa                                                                                             <  5:157601/148‑1 (MQ=255)
                                                                 cagcaCCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCa                                                                                    >  5:143543/1‑148 (MQ=255)
                                                                   gcaCCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATaa                                                                                 <  6:325565/149‑1 (MQ=255)
                                                                   gcaCCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCACAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCAt                                                                                   >  6:277468/1‑147 (MQ=255)
                                                                    caCCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAg                                                                                >  5:159538/1‑149 (MQ=255)
                                                                     aCCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGc                                                                               <  8:161971/149‑1 (MQ=255)
                                                                       cTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCt                                                                             >  3:277084/1‑149 (MQ=255)
                                                                       cTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCt                                                                             <  3:199966/149‑1 (MQ=255)
                                                                          cagcaCGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCGTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGc                                                                           <  1:126233/148‑1 (MQ=255)
                                                                          cagcaCGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCTTAGTCACCCCATAAGCCTGCa                                                                          >  6:183708/1‑149 (MQ=255)
                                                                            gcaCGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGc                                                                                                                        >  6:87608/1‑101 (MQ=255)
                                                                            gcaCGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGc                                                                                                                        <  5:87608/101‑1 (MQ=255)
                                                                             caCGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCAc                                                                         >  2:374837/1‑147 (MQ=255)
                                                                              aCGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCACAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAg                                                                      <  5:277468/149‑1 (MQ=255)
                                                                                             gTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTg                                                       <  5:422950/149‑1 (MQ=255)
                                                                                                gCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCgg                                                    <  2:11090/149‑1 (MQ=255)
                                                                                                        cGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGAt                                            <  1:99591/149‑1 (MQ=255)
                                                                                                         gCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATg                                           <  5:123415/149‑1 (MQ=255)
                                                                                                         gCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATg                                           >  6:436298/1‑149 (MQ=255)
                                                                                                           ttCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGt                                          >  3:128966/1‑148 (MQ=255)
                                                                                                              gCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCAc                                      <  6:143543/149‑1 (MQ=255)
                                                                                                                   ccGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAggg                                 >  4:333856/1‑149 (MQ=255)
                                                                                                                        aaCACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGc                            <  3:183197/149‑1 (MQ=255)
                                                                                                                           aggagccgcAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAg                         >  1:300751/7‑149 (MQ=255)
                                                                                                                           acCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAg                         <  2:362798/149‑1 (MQ=255)
                                                                                                                           acCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAg                         >  6:414386/1‑149 (MQ=255)
                                                                                                                             ccAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCa                       >  7:184019/1‑149 (MQ=255)
                                                                                                                             ccAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCa                       >  1:148736/1‑149 (MQ=255)
                                                                                                                                     aggagccgCAGGCAAGTGCCTGCAGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTg               >  3:349824/8‑149 (MQ=255)
                                                                                                                                     aGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTg               <  4:128966/149‑1 (MQ=255)
                                                                                                                                       cTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTGCa             <  8:184019/149‑1 (MQ=255)
                                                                                                                                         gccgCAGGCAAGTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACACCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCAGGCGAGCACCAGCTTGCAGc           >  3:336827/4‑149 (MQ=255)
                                                                                                                                         aCTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTGCAGc           >  2:389973/1‑149 (MQ=255)
                                                                                                                                          ccgCAGGCAAGTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGc                                                                                                   >  3:164149/3‑60 (MQ=255)
                                                                                                                                          ccgCAGGCAAGTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGc                                                                                                   <  4:164149/58‑1 (MQ=255)
                                                                                                                                          ccgCAGGCAAGTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTGCAGcc          >  1:176216/3‑149 (MQ=255)
                                                                                                                                             cAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTGCAGCCGGt       <  7:24163/149‑1 (MQ=255)
                                                                                                                                             cAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTGCAGCCGGt       >  1:190035/1‑149 (MQ=255)
                                                                                                                                              aGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTGCAGCTGGt       >  8:157478/1‑148 (MQ=255)
                                                                                                                                               ggCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTGCAGCCGGTGa     <  1:288756/149‑1 (MQ=255)
                                                                                                                                                  aCTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTGCAGCCGGTGACCa  <  5:436298/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TGCCAGTGGGTGACGTTGATGGCGGTTTCAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTGCAGCCGGTGACCA  >  NC_002947/4826907‑4827201

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: ‑