Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,789,173 G→C 9.2% A1556G (GCC→GGC)  pvdI ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,789,1730GC 9.2% 185.4 / 8.4 65A1556G (GCC→GGC) pvdInon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base G (28/31);  new base C (3/3);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.85e-01

CGGCGCCGAGCGCTTGCAGGCGATGGGCCAGGCGGTTGGCGCGCCGGTTGAGTTCGGCATAGCTCCAGGCCTGGCCTTCATGGATGACAGCCACGGCATCCGGTGCGCGCAGCACTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCGGCCGTGGCCAGGGCCAGGCGATCAACGCAGACGTGCGGCTGGGCACACAGCGCTTGCAGCAGGTTCAGCCAGTGCGTG  >  NC_002947/4789033‑4789316
                                                                                                                                            |                                                                                                                                               
cGGCGCCGAGCGCTTGCAGGCGATGGGCCAGGCGGTTGGCGCGCCGGTTGAGTTCGGCATAGCTCCAGGCCTGGCCTTCATGGATGACAGCCACGGCATCCGGTGCGCGCAGCACTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCg                                                                                                                                         <  3:256618/149‑1 (MQ=255)
  gcgcCGAGCGCTTGCAGGCGATGGGCCAGGCGGTTGGCGCGCCGGTTGAGTTCGGCATAGCTCCAGGCCTGGCCTTCTTGGATGACAGCCACGGCATCCGGTGCGCGCAGCACTTGCGCTTCGAACTGCCGTGGCCCGCCAGttcg                                                                                                                                          >  1:120697/1‑143 (MQ=255)
  gcgcCGAGCGCTTGCAGGCGATGGGCCAGGCGGTTGGCGCGCCGGTTGAGTTCGGCATAGCTCCAGGCCTGGCCTTCATGGATGACAGCCACGGCATCCGGTGCGCGCAGCACTTGCGCTTCGAACTGCCGTGGCCCGCCAGttcg                                                                                                                                          <  2:120697/146‑4 (MQ=255)
      cGAGCGCTTGCAGGCGATGGGCCAGGCGGTTGGCGCGCCGGTTGAGTTCGGCATAGCTCCAGGCCTGGCCTTCATGGATGACAGCCACGGCATCCGGTGCGCGCAGCACTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCg                                                                                                                                         >  8:183198/1‑143 (MQ=255)
      cGAGCGCTTGCAGGCGATGGGCCAGGCGGTTGGCGCGCCGGTTGAGTTCGGCATAGCTCCAGGCCTGGCCTTCATGGATGACAGCCACGGCATCCGGTGCGCGCAGCACTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCg                                                                                                                                         <  7:183198/143‑1 (MQ=255)
      cGAGCGCTTGCAGGCGATGGGCCAGGCGGTTGGCGCGCCGGTTGAGTTCGGCATAGCTCCAGGCCTGGCCTTCATGGATGACAGCCACGGCATCCGGTGCGCGCAGCACTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAAt                                                                                                                                   >  8:43359/1‑149 (MQ=255)
             ttGCAGGCGATGGGCCAGGCGGTTGGCGCGCCGGTTGAGTTCGGCATAGCTCCAGGCCTGGCCTTCATGGATGACAGCCACGGCATCCGGTGCGCGCAGCACTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTg                                                                                                                            >  6:99211/1‑149 (MQ=255)
                cAGGCGATGGGCCAGGCGGTTGGCGCGCCGGTTGAGTTCGGCATAGCTCCAGGCCTGGCCTTCATGGATGACAGCCACGGCATCCGGTGCGCGCAGCACTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCa                                                                                                                          <  3:95044/148‑1 (MQ=255)
                   gCGATGGGCCAGGCGGTTGGCGCGCCGGTTGAGTTCGGCATAGCTCCAGGCCTGGCCTTCATGGATGACAGCCACGGCGTCCGGTGCGCGCAGCACTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGcc                                                                                                                      >  4:22450/1‑149 (MQ=255)
                     gATGGGCCAGGCGGTTGGCGCGCCGGTTGAGTTCGGCATAGCTCCAGGCCTGGCCTTCATGGATGACAGCCACGGCATCCGGTGCGCGCAGCACTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGcccc                                                                                                                    >  4:78853/1‑149 (MQ=255)
                      aTGGGCCAGGCGGTTGGCGCGCCGGTTGAGTTCGGCATAGCTCCAGGCCTGGCCTTCATGGATGACAGCCACGGCATCCGGTGCGCGCAGCACTTGCGCTTCGAACTGCCGTGGCCCGCCAGtt                                                                                                                                            >  6:237289/1‑123 (MQ=255)
                      aTGGGCCAGGCGGTTGGCGCGCCGGTTGAGTTCGGCATAGCTCCAGGCCTGGCCTTCATGGATGACAGCCACGGCATCCGGTGCGCGCAGCACTTGCGCTTCGAACTGCCGTGGCCCGCCAGtt                                                                                                                                            >  6:242698/1‑123 (MQ=255)
                      aTGGGCCAGGCGGTTGGCGCGCCGGTTGAGTTCGGCATAGCTCCAGGCCTGGCCTTCATGGATGACAGCCACGGCATCCGGTGCGCGCAGCACTTGCGCTTCGAACTGCCGTGGCCCGCCAGtt                                                                                                                                            <  5:237289/124‑2 (MQ=255)
                      aTGGGCCAGGCGGTTGGCGCGCCGGTTGAGTTCGGCATAGCTCCAGGCCTGGCCTTCATGGATGACAGCCACGGCATCCGGTGCGCGCAGCACTTGCGCTTCGAACTGCCGTGGCCCGCCAGtt                                                                                                                                            <  5:242698/124‑2 (MQ=255)
                           ccAGGCGGTTGGCGCGCCGGTTGAGAGCGGCATAGCTCCAGGCCTGGCCTTCATGGATGACAGCCACGGCATCCGGTGCGCGCAGCACTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCAc                                                                                                              <  8:250564/149‑1 (MQ=255)
                               gCGGTTGGCGCGCCGGTTGAGTTCGGCATAGCTCCAGGCCTGGCCTTCATGGATGACAGCCACGGCATCCGGTGCGCGCAGCACTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGt                                                                                                          <  1:36372/149‑1 (MQ=255)
                                      gcgcgcCGGTTGAGTTCGGCATAGCTCCAGGCCTGGCCTTCATGGATGACAGCCACGGCATCCGGTGCGCGCAGCACTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTg                                                                                                                            >  5:369186/1‑124 (MQ=255)
                                      gcgcgcCGGTTGAGTTCGGCATAGCTCCAGGCCTGGCCTTCATGGATGACAGCCACGGCATCCGGTGCGCGCAGCACTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTg                                                                                                                            <  6:369186/124‑1 (MQ=255)
                                       cgcgcCGGTTGAGTTCGGCATAGCTCCAGGCCTGGCCTTCATGGATGACAGCCACGGCATCCGGTGCGCGCAGCACTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCg                                                                                                  <  3:78853/149‑1 (MQ=255)
                                        gcgcCGGTTGAGTTCGGCATAGCTCCAGGCCTGGCCTTCATGGATGACAGCCACGGCATCCGGTGCGCGCAGCACTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGa                                                                                                 >  5:351485/1‑149 (MQ=255)
                                         cgcCGGTTGAGTTCGGCATAGCTCCAGGCCTGGCCTTCATGGATGACAGCCACGGCATCCGGTGCGCGCAGCACTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGAc                                                                                                <  3:301752/149‑1 (MQ=255)
                                          gcCGGTTGAGTTCGGCATAGCTCCGGGCCTGGCCTTCATGGATGACAGCCACGGCATCCGGTGCGCGCAGCACTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGAcc                                                                                               <  4:344866/149‑1 (MQ=255)
                                                       ggCATAGCTCCAGGCCTGGCCTTCATGGATGACAGCCACGGCATCCGGTGCGCGCAGCACTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCtcgt                                                                                  >  6:200224/1‑149 (MQ=255)
                                                       ggCATAGCTCCAGGCCTGGCCTTCATGGATGACAGCCACGGCATCCGGTGCGCGCAGCACTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCtcgt                                                                                  >  2:144884/1‑149 (MQ=255)
                                                                   ggcctggcctTCATGGATGACAGCCACGGCATCCGGTGCGCGCAGCACTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCGGCCGTggc                                                                        >  2:283481/1‑147 (MQ=255)
                                                                     cctggcctTCATGGATGACAGCCACGGCATCCGGTGCGCGCAGCACTTGCGATTCGAACTGTCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCACTGGAACCAGTCTCTCTCGACCAGACGGCTCTCGTCGGCCGTggccagg                                                                    >  6:334349/1‑149 (MQ=255)
                                                                            ttCATGGATGACAGCCACGGCATCCGGTGCGCGCAGCACTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCGGCCGTGGCCAGGGCCAGGc                                                             >  6:237504/1‑149 (MQ=255)
                                                                               acggATGACAGCCACGGCATCCGGTGCGCGCAGCACTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCGGCCGTGGCCAGGGCCAGGCGAt                                                          <  5:334349/147‑1 (MQ=255)
                                                                                tGGATGACAGCCACGGCATCCGGTGCGCGCAGCACTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCGGCCGTGGCCAGGGCCAGGCGATc                                                         <  5:200224/149‑1 (MQ=255)
                                                                                  gATGACAGCCACGGCATCCGGTGCGCGCAGCACTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCGGCCGTGGCCAGGGCCAGGCGATCaa                                                       >  7:94188/1‑149 (MQ=255)
                                                                                  gATGACAGCCACGGCATCCGGTGCGCGCAGCACTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCGGCCGTGGCCAGGGCCAGGCGATCaa                                                       >  6:159868/1‑149 (MQ=255)
                                                                                     gACAGCCACGGCATCCGGTGCGCGGAGCACTTGAGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCGGCCGTgg                                                                         <  6:198599/128‑1 (MQ=255)
                                                                                     gACAGCCACGGCATCCGGTGCGCGCAGCACTAGCCCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGACACTGGCACCAGTCTCACACGACCAGACGGCGCTCGTCGGCCGTgg                                                                         >  5:198599/1‑128 (MQ=255)
                                                                                       cAGCCACAGCATCCGGTGCGCGCAGCACTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCgg                                                                               >  2:254993/1‑120 (MQ=255)
                                                                                       cAGCCACAGCATCCGGTGCGCGCAGCACTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCgg                                                                               <  1:254993/120‑1 (MQ=255)
                                                                                        aGCCACGGCATCCGGTGCGCGCAGCACTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCGGCCGTGGCCAGGGCCAGGCGATCAACGCAGa                                                 >  3:172382/1‑149 (MQ=255)
                                                                                        aGCCACGGCATCCGGTGCGCGCAGCACTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCGGCCGTGGCCAGGGCCAGGCGATCAACGCAGa                                                 >  4:318496/1‑149 (MQ=255)
                                                                                         gCCACGGCATCCGGTGCGCGCAGCACTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCGGCCGTGGCCAGGGCCAGGCGATCAACGCAGAc                                                <  8:94188/149‑1 (MQ=255)
                                                                                         gCCACGGCATCCGGTGCGCGCAGCACTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCGGCCGTGGCCAGGGCCAGGCGATCAACGCAGAc                                                >  7:177574/1‑149 (MQ=255)
                                                                                               gCATCCGGTGCGCGCAGCACTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCGGCCGTggccaggg                                                                   >  6:178115/1‑124 (MQ=255)
                                                                                               gCATCCGGTGCGCGCAGCACTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCGGCCGTggccaggg                                                                   <  5:178115/124‑1 (MQ=255)
                                                                                                cATCCGGTGCGCGCAGCACTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCGGCCGTGGCCAGGGCCAGGCGATCAACGCAGACGTGCGGc                                         <  5:237504/149‑1 (MQ=255)
                                                                                                    cGGTGCGCGCAGCACTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCGGCCGTGGCCAGGGCCAGGCGATCAACGCAGACGTGCGGCTggg                                     <  6:179708/149‑1 (MQ=255)
                                                                                                    cGGTGCGCGCAGCACTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCGGCCGTGGCCAGGGCCAGGCGATCAACGCAGACGTGCGGCTggg                                     <  1:153711/149‑1 (MQ=255)
                                                                                                     ggTGCGCGCAGCACTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCGGCCGTGGCCAGGGCCAGGCGATCAACGCAGACGTGCGGCTGGGc                                    <  5:159868/149‑1 (MQ=255)
                                                                                                      gTGCGCGCAGCACTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCGGCCGTGGCCAGGGCCAGGCGATCAACGCAGACGTGCGGCTGGGca                                   <  1:283481/149‑1 (MQ=255)
                                                                                                            gcagcaCTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCGGCCGTGGCCAGGGCCAGGCGATCAACGCAGACGTGCGGCTgg                                      >  5:265588/1‑140 (MQ=255)
                                                                                                            gcagcaCTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCGGCCGTGGCCAGGGCCAGGCGATCAACGCAGACGTGCGGCTgg                                      <  6:265588/140‑1 (MQ=255)
                                                                                                                caCTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCGGCCGTGGCCAGGGCCAGGCGATCAACGCAGACGTGCGGCTGGGCACACAGCGCtt                         <  5:279628/149‑1 (MQ=255)
                                                                                                                 aCTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCGGCCGTGGCCAGGGCCAGGCGATCAACGCAGACGTGCGGCTGGGCACACAGCGCTTg                        <  7:43359/149‑1 (MQ=255)
                                                                                                                   ttGCGCTTAGAACTGGCGTGACAGTGCAGTATGTCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCGGCCGTGGCCAGGGCCAGGCGATCAACGCAGACGTGCGGCTGGGCACACAGCGCTTgca                      <  3:61213/149‑1 (MQ=255)
                                                                                                                     gcgcTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCGGCCGTGGCCAGGGCCAGGCGATCAACGCAGACGTGCGGCTGGGCACACAGCGCTTgca                      <  8:177574/147‑1 (MQ=255)
                                                                                                                      cgcTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCGGCCGTGGCCAGGGCCAGGCGATCAACGCAGACGTGCGGCTGGGCACACAGCGCTTgcagc                    <  1:144884/148‑1 (MQ=255)
                                                                                                                      cgcTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCGGCCGTGGCCAGGGCCAGGCGATCAACGCAGACGTGCGGCTGGGCACACAGCGCTTgcagc                    <  3:22450/148‑1 (MQ=255)
                                                                                                                        cTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCGGCCGTGGCCAGGGCCAGGCGATCAACGCAGACGTGCGGCTGGGCACACAGCGCTTGCAGCAgg                 >  7:226193/1‑149 (MQ=255)
                                                                                                                            gAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCGGCCGTGGCCAGGGCCAGGCGATCAACGCAGACGTGCGGCTGGGCACACAGCGCTTGCAGCAGGTTCa             <  3:318496/149‑1 (MQ=255)
                                                                                                                            gAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCGGCCGTGGCCAGGGCCAGGCGACCAACGCAGACGTGCGGCTGGGCACACAGCGCTTGCAGCAGGtt               >  1:362391/1‑147 (MQ=255)
                                                                                                                             aaCTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCGGCCGTGGCCAGGGCCAGGCGATCAACGCAg                                                  <  7:233340/111‑1 (MQ=255)
                                                                                                                             aaCTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCGGCCGTGGCCAGGGCCAGGCGATCAACGCAg                                                  >  8:233340/1‑111 (MQ=255)
                                                                                                                             aaCTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCGGCCGTGGCCAGGGCCAGGCGATCAACGCAGACGTGCGGCTGGGCACACAGCGCTTGCAGCAGGTTCAg            >  1:175500/1‑149 (MQ=255)
                                                                                                                              aCTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCGGCCGTGGCCAGGGCCAGGCGATCAACGCAGACGTGCGGCTGGGCACACAGCGCTTGCAGCAGGTTCAGc           <  5:99211/149‑1 (MQ=255)
                                                                                                                                  gCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCGGCCGTGGCCAGGGCCAGGCGATCAACGCAGACGTGCGGCTGGGCACACAGCGCTTGCAGCAGGTTCAGCCAGt       <  6:351485/149‑1 (MQ=255)
                                                                                                                                    gggCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCGGCCGTGGCCAGGGCCAGGCGATCAACGCAGACGTGCGGCTGGGCACACAGCGCTTGCAGCAGGTTCAGCCAGTGc     >  3:289056/1‑149 (MQ=255)
                                                                                                                                      gCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCGGCCGTGGCCAGGGCCAGGCGATCAACGCAGACGTGCGGCTGGGCACACAGCGCTTGCAGCAGGTTCAGCCAGTGCGt   >  5:225617/1‑149 (MQ=255)
                                                                                                                                       ccACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCGGCCGTGGCCAGGGCCAGGCGATCAACGCAGACGTGCGGCTGGGCACACAGCGCTTGCAGCAGGTTCAGCCAGTGCGTg  >  6:41792/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                               
CGGCGCCGAGCGCTTGCAGGCGATGGGCCAGGCGGTTGGCGCGCCGGTTGAGTTCGGCATAGCTCCAGGCCTGGCCTTCATGGATGACAGCCACGGCATCCGGTGCGCGCAGCACTTGCGCTTCGAACTGGCGGGCCACGGCAGTATCGCTCAATTGCTCTGCAGGCCCCTGGCACCAGTCTCTCTCGACCAGACGGCGCTCGTCGGCCGTGGCCAGGGCCAGGCGATCAACGCAGACGTGCGGCTGGGCACACAGCGCTTGCAGCAGGTTCAGCCAGTGCGTG  >  NC_002947/4789033‑4789316

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