Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,772,509 G→T 8.2% P2253Q (CCA→CAA)  pvdD ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,772,5090GT 8.2% 220.4 / 8.1 73P2253Q (CCA→CAA) pvdDnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base G (33/34);  new base T (3/3);  total (36/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 = 5.83e-01

GCGCTGATCGTCCAGGCGCACCAGCAGCATGCGCTGCAGCGGCGGCCTGGCCAGGTCGAAGCCGCGTTCACGCTCTTCGTTGGCCAGTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGGCGCGCCACCACCGTTTCCCAGGCTTGGCGCAGGCGCTCGGCAACCAGGCCCTGTACCGGTACGCTCAACTGGTT  >  NC_002947/4772361‑4772657
                                                                                                                                                    |                                                                                                                                                    
gcgcTGATCGTCCAGGCGCACCAGCAGCATGCGCTGCAGCGGCGGCCTGGCCAGGTCGAAGCCGCGTTCACGCTCTTCGTTGGCCAGTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTg                                                                                                                                                      <  2:34468/149‑1 (MQ=255)
  gcTGATCGTCCAGGCGCACCAGCAGCATGCGCTGCAGCGGCGGCCTGGCCAGGTCGAAGCCGCGTTCACGCTCTTCGTTGGCCAGTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTgg                                                                                                                                                     >  3:159574/1‑148 (MQ=255)
    tGATCGTCCAGGCGCACCAGCAGCATGCGCTGCAGCGGCGGCCTGGCCAGGTCGAAGCCGCGTTCACGCTCTTCGTTGGCCAGTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTTGcca                                                                                                                                                  >  8:229655/1‑147 (MQ=255)
         gTCCAGGCGCACCAGCAGCATGCGCTGCAGCGGCGGCCTGGCCAGGTCGAAGCCGCGTTCACGCTCTTCGTTGGCCAGTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGcc                                                                                                                                             <  5:201677/149‑1 (MQ=255)
          tCCAGGCGCACCAGCAGCATGCGCTGCAGCGGCGGCCTGGCCAGGTCGAAGCCGCGTTCACGCTCTTCGTTGGCCAGTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTg                                                                                                                                                      >  1:132162/1‑139 (MQ=255)
          tCCAGGCGCACCAGCAGCATGCGCTGCAGCGGCGGCCTGGCCAGGTCGAAGCCGCGTTCACGCTCTTCGTTGGCCAGTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTg                                                                                                                                                      <  2:132162/139‑1 (MQ=255)
            cAGGCGCACCAGCAGCATGCGCTGCAGCGGCGGCCTGGCCAGGTCGAAGCCGCGTTCACGCTCTTCGTTGGCCAGTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCa                                                                                                                                          >  8:33488/1‑149 (MQ=255)
            cAGGCGCACCAGCAGCATGCGCTGCAGCGGCGGCCCGGCCAGGTCGAAGCCGCGTTCACGCTCTTCGTTGGCCAGTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTTGCCACTGccatt                                                                                                                                          <  1:139818/149‑4 (MQ=255)
             aGGCGCACCAGCAGCATGCGCTGCAGCGGCGGCCTGGCCAGGTCGAAGCCGCGTTCACGCTCTTCGTTGGCCAGTGCCTCCATCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTTGCCACTGccattg                                                                                                                                         <  5:3562/149‑5 (MQ=255)
              ggCGCACCAGCAGCATGCGCTGCAGCGGCGGCCTGGCCAGGTCGAAGCCGCGTTCACGCTCTTCGTTGGCCAGTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTTGCCACTGccattgc                                                                                                                                        <  7:229655/149‑6 (MQ=255)
                      agcagcATGCGCTGCAGCGGCGGCCTGGCCAGGTCGAAGCCGCGTTCACGCTCTTCGTTGGCCAGTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTg                                                                                                                                <  7:30587/149‑1 (MQ=255)
                          gcaTGCGCTGCAGCGGCGGCCTGGCCAGGTCGAAGCCGCGTTCACGCTCTTCGTTGGCCAGTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAAt                                                                                                                            >  5:77173/1‑149 (MQ=255)
                               cgcTGCAGCGGCGGCCTGGCCAGGTCGAAGCCGCGTTCACGCTCTTCGTTGGCCAGTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGc                                                                                                                                              <  6:278607/126‑1 (MQ=255)
                               cgcTGCAGCGGCGGCCTGGCCAGGTCGAAGCCGCGTTCACGCTCTTCGTTGGCCAGTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGc                                                                                                                                              >  5:278607/1‑126 (MQ=255)
                                gcTGCAGCGGCGGCCTGGCCAGGTCGAAGCCGCGTTCACGCTCTTCGTTGGCCAGTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCa                                                                                                                      >  6:184015/1‑149 (MQ=255)
                                    cAGCGGCGGCCTGGCCAGGTCGAAGCCGCGTTCACGCTCTTCGTTGGCCAGTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCa                                                                                                                  <  3:75050/149‑1 (MQ=255)
                                       cggcggCCTGGCCAGGTCGAAGCCGCGTTCACGCTCTTCGTTGGCCAGTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATTGCAGTGGCAACTTCCCCATAACGTACTGCAATGGCTCATACAgg                                                                                                                >  6:48813/1‑148 (MQ=255)
                                              cTGGCCAGGTCGAAGCCGCGTTCACGCTCTTCGTTGGCCAGTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGcc                                                                                                         >  8:20525/1‑148 (MQ=255)
                                                  ccAGGTCGAAGCCGCGTTCACGCTCTTCGTTGGCCAGTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAgg                                                                                                    <  7:33488/149‑1 (MQ=255)
                                                  ccAGGTCGAAGCCGCGTTCACGCTCTTCGTTGGCCAGTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAAATGCTCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAg                                                                                                     >  3:52000/1‑148 (MQ=255)
                                                   cAGGTCGAAGCCGCGTTCACGCTCTTCGTTGGCCAGTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAgg                                                                                                    >  5:288973/1‑148 (MQ=255)
                                                   cAGGTCGAAGCCGCGTTCACGCTCTTCGTTGGCCAGTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGa                                                                                                   <  2:94580/149‑1 (MQ=255)
                                                      gTCGAAGCCGCGTTCACGCTCTTCGTTGGCCAGTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGc                                                                                                          <  6:247384/139‑1 (MQ=255)
                                                      gTCGAAGCCGCGTTCACGCTCTTCGTTGGCCAGTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGCCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGc                                                                                                          >  5:247384/1‑139 (MQ=255)
                                                          aaGCCGCGTTCACGCTCTTCGTTGGCCAGTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGt                                                                                            <  4:159574/149‑1 (MQ=255)
                                                             ccGCGTTCACGCTCTTCGTTGGCCAGTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCa                                                                                                                                       >  4:197906/1‑103 (MQ=255)
                                                             ccGCGTTCACGCTCTTCGTTGGCCAGTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCa                                                                                                                                       <  3:197906/103‑1 (MQ=255)
                                                              cgcgTTCACGCTCTTCGTTGGCCAGTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTgcg                                                                                         <  7:115417/148‑1 (MQ=255)
                                                                 gTTCACGCTCTTCGTTGGCCAGTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAAt                                                                                     <  7:20525/149‑1 (MQ=255)
                                                                 gTTCACGATCTTCTTTGGCCAGTGCCTCCAGCGGCTGCGCGCAATCTCCACGACCGCGCCAGTCCAGCTGACGCAATGGCAGTGGCAACTGCCCCATCACGTAGTTCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAAt                                                                                     <  3:45098/149‑1 (MQ=255)
                                                                       gCTCTTCGTTGGCCAGTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATAcgtcg                                                                               >  7:199807/1‑149 (MQ=255)
                                                                       gCTCTTCGTTGGCCAGTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATAcgtcg                                                                               >  5:299592/1‑149 (MQ=255)
                                                                             cGTTGGCCAGTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTAACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAATCCGTTGCGCAATACGTCGTGgcgc                                                                         >  4:53163/1‑149 (MQ=255)
                                                                              gTTGGCCAGTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGgcgcg                                                                        >  4:42956/1‑149 (MQ=255)
                                                                              gTTGGCCAGTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGgcgc                                                                         <  6:288973/148‑1 (MQ=255)
                                                                              gTTGGCCAGTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGgcgc                                                                         <  5:184015/148‑1 (MQ=255)
                                                                              gTTGGCCAGTGCCTCCAGCGCCTGCGCGAAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGTCAACTGCCCCATCAAGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGgcgc                                                                         <  5:48813/148‑1 (MQ=255)
                                                                               ttGGCCAGTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGgcgcgc                                                                       >  5:201957/1‑149 (MQ=255)
                                                                                  gCCAGTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGGCGCGccac                                                                    >  7:37403/1‑149 (MQ=255)
                                                                                   ccAGTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGGCGCGccacc                                                                   >  7:133236/1‑149 (MQ=255)
                                                                                   ccAGTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGGCGCGccac                                                                    >  2:103079/1‑148 (MQ=255)
                                                                                      gTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGGCGCGCCaccacc                                                                <  6:299592/149‑1 (MQ=255)
                                                                                          cTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATg                                                                                                                           <  6:71642/86‑1 (MQ=255)
                                                                                           tCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGGCGCGCCACCACCGTTTc                                                           <  6:163942/149‑1 (MQ=255)
                                                                                             cAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGGCGCGCCACCACCGTTTcc                                                          >  4:102647/1‑148 (MQ=255)
                                                                                              aGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGGCGCGCCACCACCGTTTcc                                                          >  3:239914/1‑147 (MQ=255)
                                                                                                     gcgcgcAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGGCGCGCCACCACCGTTTCCCAGGCTTgg                                                 >  4:277349/1‑149 (MQ=255)
                                                                                                          cAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGGCGCGCCACCACCGTTTCCCAGGCTTGGCGCAg                                            <  8:37403/149‑1 (MQ=255)
                                                                                                          cAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCGCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGGCGCGCCACCACCGTTTCCCAGGCTTGGCGCAg                                            <  5:216508/149‑1 (MQ=255)
                                                                                                          cAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAAATGCTCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGGCGCGCCACCACCGTTTCCCAGGCTTGGCGCAg                                            <  4:52000/149‑1 (MQ=255)
                                                                                                           aaTCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGGCGCGCCACCACCGTTTCCCAGGCTTGGCGCAgg                                           <  3:42956/149‑1 (MQ=255)
                                                                                                              cGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGGCGCGCCACCACCGTTTc                                                           <  1:129124/130‑1 (MQ=255)
                                                                                                              cGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGGCGCGCCACCACCGTTTc                                                           >  2:129124/1‑130 (MQ=255)
                                                                                                               gCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGGCGCGCCACCACCGTTTCCCAGGCTTGGCGCAGGCGCt                                       <  6:201957/149‑1 (MQ=255)
                                                                                                                ccACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGGCGCGCCACCACCGTTTCCCAGGCTTGGCGCAGGCGCTc                                      >  2:118394/1‑149 (MQ=255)
                                                                                                                 cACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGGCGCGCCACCACCGTTTCCCAGGCTTGGCGCAGGCGCTCg                                     <  3:277349/149‑1 (MQ=255)
                                                                                                                   cGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGGCGCGCCACCACCGTTTCCCAGGCTTGGCGCAGGCGCTCGGc                                   <  3:102647/149‑1 (MQ=255)
                                                                                                                      ccGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGGCGCGCCACCACCGTTTCCCAGGCTTGGCGCAGGCGCTCGGCAAc                                <  1:103079/149‑1 (MQ=255)
                                                                                                                       cgcgCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGGCGCGCCACCACCGtt                                                             <  1:271070/119‑1 (MQ=255)
                                                                                                                       cgcgCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGGCGCGCCACCACCGtt                                                             >  2:271070/1‑119 (MQ=255)
                                                                                                                       cgcgCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGGCGCGCCACCACCGTTTCCCAGGCTTGGCGCAGGCGCTCg                                     <  4:219358/143‑1 (MQ=255)
                                                                                                                       cgcgCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGGCGCGCCACCACCGTTTCCCAGGCTTGGCGCAGGCGCTCg                                     >  3:219358/1‑143 (MQ=255)
                                                                                                                          gcCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGGCGCGCCACCACCGTTTCCCAGGCTTGGCGCAGGCGCTCGGCAACCAgg                            >  6:1795/1‑149 (MQ=255)
                                                                                                                           ccAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGGCGCGCCACCACCGTTTCCCAGGCTTGGCGCAGGCGCTCGGCAACCAGGc                           >  5:28414/1‑149 (MQ=255)
                                                                                                                                    cTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGAGCAATACGTCGTGGCGCGCCACCACCGTTTCCCATGCTTGGCGCAGGCGCTCGGCAACCAGGCCCTGTACCg                  <  3:53163/149‑1 (MQ=255)
                                                                                                                                    cTCACGCAATGGCAGTGGCAAATGCCCCATCACGTACTTCAATGGCGCATCCAGGTCTTTCCAAAGGAAGCCGGTGCGCAATACGTCGTGGCGCGCCACCACCGTTTCCCAGGCTTGGCGCAGGCGCTCGGCAACCAGGCCCTGTACCg                  <  8:134068/149‑1 (MQ=255)
                                                                                                                                       aCGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGGCGCGCCACCACCGTTTCCCAGGCTTGGCGCAGGCGCTCGGCAACCAGGCCCTGTACCgg                 <  5:107594/147‑1 (MQ=255)
                                                                                                                                       aCGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGGCGCGCCACCACCGTTTCCCAGGCTTGGCGCAGGCGCTCGGCAACCAGGCCCTGTACCGGTa               >  1:27922/1‑149 (MQ=255)
                                                                                                                                       aCGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGGCGCGCCACCACCGTTTCCCAGGCTTGGCCCAGGCGCTCGGCAACCAGGCCCTGTACCGGTa               >  8:248061/1‑149 (MQ=255)
                                                                                                                                         gatggggCAGTTGCCACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGGCGCGCCACCACCGTTTCCCAGGCTTGGCGCAGGCGCTCGGCAACCAGGCCCTGTACCGGTACg             >  5:63645/6‑149 (MQ=255)
                                                                                                                                         gCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGGCGCGCCACCACCGTTTCCCAGGCTTGGCGCAGGCGCTCGGCAACCAGGCCCTGTACCGGTACg             >  3:230489/1‑149 (MQ=255)
                                                                                                                                           aaTGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGGCGCGCCACCACCGTTTCCCAGGCTTGGCGCAGGCGCTCGGCAACCAGGCCCTGTACCGGTACGCt           >  7:60579/1‑149 (MQ=255)
                                                                                                                                           aaTGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGGCGCGCCACCACCGTTTCCCAGGCTTGGCGCAGGCGCTCGGCAACCAGGCCCTGTACCGGTACGCt           <  8:60579/149‑1 (MQ=255)
                                                                                                                                             gggCAGTTGCCACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGGCGCGCCACCACCGTTTCCCAGGCTTGGCGCAGGCGCTCGGCAACCAGGCCCTGTACCGGTACGCTc          >  4:152859/2‑148 (MQ=255)
                                                                                                                                              ggCAGTGGCAACTGCCCAATCACGTACTGCAATGTCTCATCCAGGTCTTGCAAAAGGAATCCGGTGCGCAATACGTCGTGGCGAGCCACCAACGTATCACAGGCTTGTCGCAGAAGCTCTGAAACCCGGACCTTAACCGGTACgaaaa         >  4:32629/1‑144 (MQ=255)
                                                                                                                                               gCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGGCGCGCCACCACCGTTTCCCAGGCTTGGCGCAGGCGCTCGGCAACCAGGCCCTGTACCGGTACGCTCAAc       <  8:199807/149‑1 (MQ=255)
                                                                                                                                                   tGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAG‑CTGTGCGCAATACGTCGTGGCGCGCCACCACCGTTTCCCAGGCTTGGCGCAGGCGCTCGGCAACCAGGCCCTGTACCGGTACGCTCAACTGGtt  >  2:22282/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCGCTGATCGTCCAGGCGCACCAGCAGCATGCGCTGCAGCGGCGGCCTGGCCAGGTCGAAGCCGCGTTCACGCTCTTCGTTGGCCAGTGCCTCCAGCGCCTGCGCGCAATCGCCACGACCGCGCCAGTCCAGCTCACGCAATGGCAGTGGCAACTGCCCCATCACGTACTGCAATGGCTCATCCAGGTCTTGCCAAAGGAAGCCGGTGCGCAATACGTCGTGGCGCGCCACCACCGTTTCCCAGGCTTGGCGCAGGCGCTCGGCAACCAGGCCCTGTACCGGTACGCTCAACTGGTT  >  NC_002947/4772361‑4772657

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