Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,788,412 T→A 8.0% M1825L (ATG→TTG)  pvdI ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,788,4120TA 8.0% 178.5 / 4.8 60M1825L (ATG→TTG) pvdInon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base T (28/27);  new base A (2/3);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

GTGGCAGCGCCTGCGCCGTGGCACCTAGGTAAGGCGCGCAGTCGTGGGCGGTAACCGTGACGGTCGCCTCAGTCGGGCCGTAGGTATTGAGCAGCGTTACATGGCCAAGGCCCGCCGCCTGCCAGGCGAGCAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGAGGGCCGTGCTCGGCAAAGTCCTTGGCCAGCATGAACCAGTAGGCGGTGGTCACATCTACCACGCTGATGTCATTGCTGATCAGTTCGCGGTAGAACCGCTCA  >  NC_002947/4788265‑4788560
                                                                                                                                                    |                                                                                                                                                    
gTGGCAGCGCCTGCGCCGTGGCACCTAGGTAAGGCGCGCAGTCGGGGGCGGTAACCGTGACGGTCGCCTCAGTCGGGCCGTAGGTATTGAGCAGCGTTACATGGCCAAGGCCCGTCGCCTGCCAGGCGAGCAAGCCCTCCGGCGGCATgg                                                                                                                                                    <  1:932756/150‑1 (MQ=255)
    cAGCGCCTGCGCCGTGGCACCTAGGTAAGGCGCGCAGTCGTGGGCGGTAACCGTGACGGTCGCCTCAGTCGGGCCGTAGGTATTGAGCAGCGTTACATGGCCAAGGCCCGCCGCCTGCCAGGCGAGCAAGCCCTCCGGCGGCATGGCTTc                                                                                                                                                <  2:1063052/150‑1 (MQ=255)
     aGCGCCTGCGCCGTGGCACCTAGGTAAGGCGCGCAGTCGTGGGCGGTAACCGTGACGGTCGCCTCAGTCGGGCCGTAGGTATTGAGCAGCGTTACATGGCCAAGGCCCGCCGCCTGCCAGGCGAGCAAGCCCTCCGGCGGCATGGCTTcg                                                                                                                                               >  2:359732/1‑150 (MQ=255)
                    gCACCTAGGTAAGGCGCGCAGTCGTGGGCGGTAACCGTGACGGTCGCCTCAGTCGGGCCGTAGGTATTGAGCAGCGTTACATGGCCAAGGCCCGCCGCCTGCCAGGCGAGCAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGt                                                                                                                                     <  2:169781/145‑1 (MQ=255)
                    gCACCTAGGTAAGGCGCGCAGTCGTGGGCGGTAACCGTGACGGTCGCCTCAGTCGGGCCGTAGGTATTGAGCAGCGTTACATGGCCAAGGCCCGCCGCCTGCCAGGCGAGCAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGt                                                                                                                                     >  1:169781/1‑145 (MQ=255)
                       ccTAGGTAAGGCGCGCAGTCGTGGGCGGTAACCGTGACGGTCGCCTCAGTCGGGCCGTAGGTATTGAGCAGCGTTACATGGCCAAGGCCCGCCGCCTGCCAGGCGAGCAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAACTgg                                                                                                                             <  2:917163/150‑1 (MQ=255)
                        cTAGGTAAGGCGCGCAGTCGTGGGCGGTAACCGTGACGGCCGCTTCACTCTTGCCGTAGGTATTGAGCAGAATAACGTGGCCAAGGCCCGCCGCCTGCCAGGCGAGCAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAACTGGt                                                                                                                            <  1:74388/150‑1 (MQ=255)
                                 gcgcgcAGTCGTGGGCGGTAATCGTGACGGTCGCCTCAGTCGGGCCGTAGGTATTGAGCAGCGTTACATGGCCAAGGCCCGCCGCCTGCCAGGCGAGCAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAACTGGTGCAggcggc                                                                                                                   >  1:1108618/1‑150 (MQ=255)
                                 gcgcgcAGTCGTGGGCGGTAACCGTGACGGTCGCCTCAGTCGGGCCGTAGGTATTGTGCAGCGTTACATGGCCAAGGCCCGCCGCCTGCCAGGCGAGCAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAAc                                                                                                                                <  1:619933/137‑1 (MQ=255)
                                 gcgcgcAGTCGTGGGCGGTAACCGTGACGGTCGCCTCAGTCGGGCCGTAGGTATTGAGCAGCGTTACATGGCCAAGGCCCGCCGCCTGCCAGGCGAGCAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAAc                                                                                                                                >  2:619933/1‑137 (MQ=255)
                                 gcgcgcAGTCGTGGGCGGTAACCGTGACGGTCGCCTCAGTCGGGCCGTAGGTATTGAGCAGCGTTACATGGCCAAGGCCCGCCGCCTGCCAGGCGAGCAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAACTGGTGCAtgcggc                                                                                                                   >  1:296504/1‑150 (MQ=255)
                                 gcgcgcAGTCGTGGGCGGTAACCGTGACGGTCGCCTCAGTCGGGCCGTAGGTATTGAGCAGCGTTACATGGCCAAGGCCCGCCGCCTGCCAGGCGAGCAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAACTGGTGCAggcggc                                                                                                                   >  2:1075050/1‑150 (MQ=255)
                                  cgcgcAGTCGTGGGCGGTAACCGTGACGGTCGCCTCAGTCGGGCCGTAGGTATTGAGCAGCGTTACATGGCCAAGGCCCGCCGCCTGCCAGGCGAGCAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAACTGGTGCATGCGGcc                                                                                                                  <  2:296504/150‑1 (MQ=255)
                                  cgcgcAGTCGTGGGCGGTAACCGTGACGGTCGCCTCAGTCGGGCCGTAGGTATTGAGCAGCGTTACATGGCCAAGGCCCGCCGCCTGCCAGGCGAGCAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAACTGGTGCAGGCGGcc                                                                                                                  <  1:1075050/150‑1 (MQ=255)
                                      cAGTCGTGGGCGGTAACCGTGACGGTCGCCTCAGTCGGGCCGTAGGTATTGAGCAGCGTTACATGGCCAAGGCCCGCCGCCTGCCAGGCGAGCAAGCCCTCCGGCGGCCACTCTCCTCCGCCGGCGTGAAACTGGTGCAGGCGGCCGTAg                                                                                                              >  2:1112409/1‑150 (MQ=255)
                                               gCGGTAACCGTGACGGTCGCCTCAGTCGGGCCGTAGGTATTGAGCAGCGTTACATGGCCAAGGCCCGCCGCCTGCCAGGCGAGCAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGAggg                                                                                                     <  1:1032887/150‑1 (MQ=255)
                                                 ggTAACCGTGACGGTCGCCTCAGTCGGGCCGTAGGTATTGAGCAGCGTTACATGGCCAAGGCCCGCCGCCTGCCAGGCGAGCAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGAGGGcc                                                                                                   <  1:1063693/150‑1 (MQ=255)
                                                     aCCGTGACGGTCGCCTCAGTCGGGCCGTAGGTATTGAGCAGCGTTACATGGCCAAGGCCCGCCGCCTGCCAGGCGAGCAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGAGGGCCGTGc                                                                                               >  2:300882/1‑150 (MQ=255)
                                                        gTGACGGTCGCCTCAGTCGGGCCGTAGGTATTGAGCAGCGTTACATGGCCAAGGCCCGCCGCCTGCCAGGCGAGCAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGAGGGCCGTGCTCg                                                                                            <  2:756988/150‑1 (MQ=255)
                                                             ggTCGCCTCAGTCGGGCCGTAGGTATTGAGCAGCGTTACATGGCCAAGGCCCGCCGCCTGCCAGGCGAGCAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGAGGGCCGTGCTCGGCaaa                                                                                       <  1:1069931/150‑1 (MQ=255)
                                                              gTCGCCTCAGTCGGGCCGTAGGTATTGAGCAGCGTTACATGGCCAAGGCCCGCCGCCTGCCAGGCGAGCAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGt                                                                                                                                     <  2:795782/103‑1 (MQ=255)
                                                              gTCGCCTCAGTCGGGCCGTAGGTATTGAGCAGCGTTACATGGCCAAGGCCCGCCGCCTGCCAGGCGAGCAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGt                                                                                                                                     >  1:795782/1‑103 (MQ=255)
                                                              gTCGCCTCAGTCGGGCCGTAGGTATTGAGCAGCGTTACATGGCCAAGGCCCGCCGCCTGCCAGGCGAGCAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGAGGGCCGTGCTCGGCAAAg                                                                                      <  2:296125/150‑1 (MQ=255)
                                                                               gTAGGTATTGAGCAGCGTCACATCAACAAACAATGCCCCCTTGCAGGCCAGTAAGGTCTCCGTCGGCATTGCTTTTCTGCCGTCGTTAAACTGGTGAAGGCGGCCGTAGGCGCGAGGGCCGTGCTCGGCAAAGTCCTTGGCCAGCATGaa                                                                     >  2:11905/1‑150 (MQ=255)
                                                                                tAGGTATTGAGCAGCGTTACATGGCCAAGGCCCGCCGCCTGCCAGGCGAGCAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGAGGGCCGTGCTCGGCAAAGTCCTTGGCCAGCATGAAc                                                                    <  2:19436/150‑1 (MQ=255)
                                                                                  ggTATTGAGCAGCGTTACATGGCCAAGGCCCGCCGCCTGCCAGGCGAGCAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCTTGAAACTGGTGCAGGCGGCCGTAGTCGCGAGGGCCGTGCTCGGCAAAGTCCTTGGCCAGCATGAACCa                                                                  <  2:637147/150‑1 (MQ=255)
                                                                                   gTATTGAGCAGCGTTACATGGCCAAGGCCCGCCGCCTGCCAGGCGAGCAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGAGGGTCGTGCTCGGCAAAGTCCTTGGCCAGCATGAACCAg                                                                 <  1:1112409/150‑1 (MQ=255)
                                                                                      ttGAGCAGCGTTACATGGCCAAGGCCCGCCGCCTGCCAGGCGAGCAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGAGGGCCGTGCTCGGCAAAGTCCTTGGCCAGCATGAACCAGTAg                                                              <  1:1118158/150‑1 (MQ=255)
                                                                                          gcagcGTTACATGGCCAAGGCCCGCCGCCTGCCAGGCGAGCAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGAGGGCCGTGCTCGGCAAAGTCCTTGGCCAGCATGAACCAGTAGGCgg                                                          >  1:719744/1‑150 (MQ=255)
                                                                                            agcGTTACATGGCCAAGGCCCGCCGCCTGCCAGGCGAGCAAGCCCTCCGGCGGCATGGCTTTGCCGCCCGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGATGGCCGTGCTCGGCAAAGTCCTTGGCCGGCATGAACCAGGAGGCggtg                                                        >  2:254771/1‑150 (MQ=255)
                                                                                            agcGTTACATGGCCAAGGCCCGCCGCCTGCCAGGCGAGCAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGAGGGCCGTGCTCGGCAAAGTCCTTGGCCAGCATGAACCAGTAGGCggtg                                                        >  2:543906/1‑150 (MQ=255)
                                                                                                        ccAAGGCCCGCCGCCTGCCAGGCGAGCAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGAGGGCCGTGCTCGGCAAAGTCCTTGGCCAGCATGAACCAGTAGGCGGTGGTCACATCTAcc                                            >  2:802537/1‑150 (MQ=255)
                                                                                                        ccAAGGCCCGCCGCCTGCCAGGCGAGCAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGAGGGCCGTGCTCGGCAAAGTCCTTGGCCAGCATGAACCAGTAGGCGGTGGTCACATCTAcc                                            >  2:388620/1‑150 (MQ=255)
                                                                                                         cAAGGCCCGCCGCCTGCCAGGCGAGCAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGAGGGCCGTGCTCGGCAAAGTCCTTGGCCAGCATGAACCAGTAGGCggtggt                                                      >  2:1021552/1‑139 (MQ=255)
                                                                                                         cAAGGCCCGCCGCCTGCCAGGCGAGCAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGAGGGCCGTGCTCGGCAAAGTCCTTGGCCAGCATGAACCAGTAGGCggtggt                                                      <  1:1021552/139‑1 (MQ=255)
                                                                                                           aGGCCCGCCGCCTGCCAGGCGAGCAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGAGGGCCGTGCTCGGCAAAGTCCTTGGCCAGCATGAACCAGTAGGCGGTGGTCACATCTACCACg                                         >  2:645106/1‑150 (MQ=255)
                                                                                                            ggCCCGCCGCCTGCCAGGCGAGCAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGAGGGCCGTGGTCGGCAAAGTCCTTGGCCAGCATGAACCAGTAGGCGGTGGTCACATCTACCACGc                                        >  1:986831/1‑150 (MQ=255)
                                                                                                                cgccgcCTGCCAGGCGAGCAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGAGGGCCGTGCTCGGCAAAGTCCTTGGCCAGCATGAACCAGTAGGCGGTGGTCACATCTACCACGCTGAt                                    <  1:543906/150‑1 (MQ=255)
                                                                                                                 gccgccTGCCAGGCGAGCAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGAgg                                                                                                      <  2:335742/83‑1 (MQ=255)
                                                                                                                 gccgccTGCCAGGCGAGCAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGAgg                                                                                                      >  1:335742/1‑83 (MQ=255)
                                                                                                                       ggCCAGGCGAGCAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGAGGGCCGTGCTCGGCAAAGTCCTTGGCCAGCATGAACCAGTAGGCGGTGGTCACATCTACCACGCTGATGTCATTg                             <  2:717381/149‑1 (MQ=255)
                                                                                                                        gcctgcaccagtttcacgCCGGCGGCGAAGCCATGCCGCCGGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGAGGGCCGTGCTCGGCAAAGTCCTTGGCCAGCATGAACCAGTAGGCGGTGGTCACATCTACCACGCTGATGTCATTGc                            <  1:707114‑M1/132‑1 (MQ=255)
                                                                                                                           tgcaccagtttcacgCCGGCGGCGAAGCCATGCCGCCGGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGAGGGCCGTGCTCGGCAAAGTCCTTGGCCAGCATGAACCAGTAGGCGGTGGTCACATCTACCACGCTGATGTCATTGCTGa                         <  2:222386‑M1/135‑1 (MQ=255)
                                                                                                                           aGGCGAGCAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGAGGGCCGTGCTCGGCAAAGTCCTTGGCCAGCATGAACCAGTAGGCGGTGGTCACATCTACCACGCTGATGTCATTGCTGa                         >  1:680624/1‑150 (MQ=255)
                                                                                                                           aGGCGAGCAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGAGGGCCGTGCTCGGCAAAGTCCTTGGCCAGCATGAACCAGTAGGCGGTGGTCACATCTACCACGCTGATGTCATTGCTGa                         <  1:388620/150‑1 (MQ=255)
                                                                                                                            ggCGAGCAAGCCCTCCGGCGGCATTGCTTCGCCGCCGGCGTGAAACTCGTGCAGGCGGCCGTAGTCGCGAGGGCCGGGCCCGGCAAAGTCCTTGGCCAGCATGAACCAGTAGGCGGTGGTCACATCTACCACGCTGATGTCATTGCTGac                        <  2:1097086/150‑2 (MQ=255)
                                                                                                                            ggCGAGCAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGAGGGCCGTGCTCGGCAAAGTCCTTGGCCAGCATGAACCAGTAGGCGGTGGTCACATCTACCACGCTGATGTCATTGCTGAt                        >  2:1006367/1‑150 (MQ=255)
                                                                                                                               ccagtttcacgCCGGCGGCGAAGCCATGCCGCCGGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGAGGGCCGTGCTCGGc                                                                                          >  1:797485‑M1/12‑81 (MQ=255)
                                                                                                                               ccagtttcacgCCGGCGGCGAAGCCATGCCGCCGGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGAGGGCCGTGCTCGGc                                                                                          <  2:797485‑M1/70‑1 (MQ=255)
                                                                                                                                  cAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGAGGGCCGTGCTCGGCAAAGTCCTTGGCCAGCATGAACCAGTAGGCGGTGGTCACATCTACCACGCTGATGTCATTGCTGATCAGTTc                  >  2:840352/1‑150 (MQ=255)
                                                                                                                                    aGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGAGGGCCGTGCTCGGCAAAGTCCTTGGCCAGCATGAACCAGTAGGCGGTGGTCACATCTACCACGCTGATGTCATTGCTGATCAGTTcgc                >  1:1079324/1‑150 (MQ=255)
                                                                                                                                     gCCCTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAACTGGTGCGGGCGCGAGAGGTATGGCGGGTCTGGGTCTGCAAAGTCCTTGGCCAGCATGAACCAGTAGGCGGTGGTCACATCTACCACGCTGATGTCATTGCTGATCAGTTcgcg               <  1:407550/150‑1 (MQ=255)
                                                                                                                                     gCCCTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGAGGGCCGTGCTCGGCAAAGTCCTTGGCCAGCATGAACCAGTAGGCGGTGGTCACATCTACCACGCTGATGTCATTGCTGATCAGTTcgcg               <  1:840352/150‑1 (MQ=255)
                                                                                                                                        cTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGAGGGCCGTGCTCGGCAAAGTCCTTGGCCAGCATGAACCAGTAGGCGGTGGTCACATCTACCACGCTGATGTCAtt                              <  1:489470/132‑1 (MQ=255)
                                                                                                                                        cTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGAGGGCCGTGCTCGGCAAAGTCCTTGGCCAGCATGAACCAGTAGGCGGTGGTCACATCTACCACGCTGATGTCAtt                              >  2:489470/1‑132 (MQ=255)
                                                                                                                                        cTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGAGGGCCGTGCTCGGCAAAGTCCTTGGCCAGCATGAACCAGTAGGCGGTGGTCACATCTACCACGCTGATGTCATTGCTGATCAGTTCGCGGTa            >  2:296347/1‑150 (MQ=255)
                                                                                                                                        cTCCGGCGGCATGGCTTCGCCGCCAGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGAGGGCCGTGCTCGGCAAAGTCCTTGGCCAGCATGAACCAGTAGGCGGTGGTCACATCTACCACGCTGATGTCATTGCTGATCAGTTCGCGGTa            >  2:296435/1‑150 (MQ=255)
                                                                                                                                               GGCATGGCTTCGCCGCCGGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGAGGGCCGTGCTCGGCAAAGTCCTTGGCCAGCATGAACCAGTAGGCGGTGGTCACATCTACCACGCTGATGTCATTGCTGATCAGTTCGCGGTAGAACCGc     >  2:987948‑M1/1‑150 (MQ=255)
                                                                                                                                                   aTGGCTTCGCCGCCGGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGAGGGCCGTGCTCGGCAAAGTCCTTGGCCAGCATGAACCAGTAGGCGGTGGTCACATCTACCACGCTGATGTCATTGCTGATCAGTTCGCGGTAGAACCGCTCa  >  2:981380/1‑150 (MQ=255)
                                                                                                                                                   aTGGCTTCGCCGCCGGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGAGGGCCGTGCTCGGCAAAGTCCTTGGCCAGCATGAACCAGTAGGCGGTGGTCACATCTACCACACTGATGTCATTGCTGATCAGTTCGCGGTAGAACCGCTCa  >  1:1059226/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTGGCAGCGCCTGCGCCGTGGCACCTAGGTAAGGCGCGCAGTCGTGGGCGGTAACCGTGACGGTCGCCTCAGTCGGGCCGTAGGTATTGAGCAGCGTTACATGGCCAAGGCCCGCCGCCTGCCAGGCGAGCAAGCCCTCCGGCGGCATGGCTTCGCCGCCGGCGTGAAACTGGTGCAGGCGGCCGTAGTCGCGAGGGCCGTGCTCGGCAAAGTCCTTGGCCAGCATGAACCAGTAGGCGGTGGTCACATCTACCACGCTGATGTCATTGCTGATCAGTTCGCGGTAGAACCGCTCA  >  NC_002947/4788265‑4788560

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