Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,783,595 C→A 5.3% W1255L (TGG→TTG)  pvdJ ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,783,5950CA 5.3% 246.4 / 2.3 76W1255L (TGG→TTG) pvdJnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base C (34/38);  new base A (2/2);  total (36/40)
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.55e-01

GCATGCGCCGCCTGGCGCAGGTTGCATTCGCCTTGCGGGTTATCGCAAGGCAGCTCGCCCGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGA  >  NC_002947/4783451‑4783744
                                                                                                                                                |                                                                                                                                                     
gCATGCGCCGCCTGGCGCAGGTTGCATTCGCCTTGCGGGTTATCGCAAGGCAGCTCGCCCGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTgg                                                                                                                                                  <  1:665954/150‑1 (MQ=255)
  aTGCGCCGCCTGGCGCAGGTTGCATTCGCCTTGCGGGTTATCGCAAGGCAGTTCGCCCGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGgcg                                                                                                                                                >  1:551643/1‑150 (MQ=255)
    gcgcCGCCTGGCGCAGGTTGCATTCGCCTTGCGGGTTATCGCAAGGCAGCTCGCCCGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCa                                                                                                                                              >  2:889869/1‑150 (MQ=255)
     cgccgcCTGGCGCAGGTTGCATTCGCCTTGCGGGTTATCGCAAGGCAGCTCGCCCGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAg                                                                                                                                             <  2:443906/150‑1 (MQ=255)
      gccgccTGGCGCAGGTTGCATTCGCCTTGCGGGTTATCGCAAGGCAGCTCGCCCGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAgg                                                                                                                                            <  2:468978/150‑1 (MQ=255)
              gcgcAGGTTGCATTCGCCTTGCGGGTTATCGCAAGGCAGCTCGCCCGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGa                                                                                                                                    >  1:250009/1‑150 (MQ=255)
                 cAGGTTGCATTCGCCTTGCGGGTTATCGCAAGGCAGCTCGCCCGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCt                                                                                                                                 <  1:889869/150‑1 (MQ=255)
                    gTTGCATTCGCCTTGCGGGTTATCGCAAGGCAGCTCGCCCGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGAc                                                                                                                                   >  1:1163362/1‑145 (MQ=255)
                    gTTGCATTCGCCTTGCGGGTTATCGCAAGGCAGCTCGCCCGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGAc                                                                                                                                   <  2:1163362/145‑1 (MQ=255)
                    gTTGCATTCGCCTTGCGGGTTATCGCAAGGCAGCTCGCCCGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGa                                                                                                                              >  1:376143/1‑150 (MQ=255)
                         aTTCGCCTTGCGGGTTATCGCAAGGCAGCTCGCCCGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTg                                                                                                                                                   >  2:537532/1‑124 (MQ=255)
                         aTTCGCCTTGCGGGTTATCGCAAGGCAGCTCGCCCGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTg                                                                                                                                                   <  1:537532/124‑1 (MQ=255)
                         aTTCGCCTTGCGGGTTATCGCAAGGCAGCTCGCCCGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCa                                                                                                                         >  2:712134/1‑150 (MQ=255)
                         aTTCGCCTTGCGGGTTATCGCAAGGCAGCTCGCCCGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCa                                                                                                                         <  1:303059/150‑1 (MQ=255)
                         aTTCGCCTTGCGGGTTATCGCAAGGCAGCTCGCCCGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTGCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCa                                                                                                                         >  1:1168359/1‑150 (MQ=255)
                              ccTCGCGGGTCATGGCAAGCCAGCCCTCCCGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTc                                                                                                                    <  2:661780/150‑1 (MQ=255)
                                     ggTTATCGCAAGGCAGCTCGCCCGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCTCGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCt                                                                                                             >  2:126347/1‑150 (MQ=255)
                                         aTCGCAAGGCAGCTCGCCCGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGccc                                                                                                         <  1:488691/150‑1 (MQ=255)
                                          tCGCAAGGCAGCTCGCCCGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGa                                                                                                                              >  1:560117/1‑128 (MQ=255)
                                          tCGCAAGGCAGCTCGCCCGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGa                                                                                                                              <  2:560117/128‑1 (MQ=255)
                                             cAAGGCAGCTCGCCCGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAgg                                                                                                     <  1:756876/150‑1 (MQ=255)
                                                ggCAGCTCGCCCGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTa                                                                                                  >  1:510587/1‑150 (MQ=255)
                                                 gCAGCTCGCCCGCGGGCTGCTGCGCGGCCTGGCGCCAGTAATCCCGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGCGTGCCGACCTGGCCGGCAGGCTCACCGCCTGCCCTAGGCTAg                                                                                                 >  2:174626/1‑150 (MQ=255)
                                                 gCAGCTCGCCCGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAg                                                                                                 >  1:63408/1‑150 (MQ=255)
                                                      tCGCCCGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGc                                                                                                <  2:1359854/146‑1 (MQ=255)
                                                      tCGCCCGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGc                                                                                                >  1:1359854/1‑146 (MQ=255)
                                                      tCGCCCGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCctgc                                                                                            <  2:1350091/150‑1 (MQ=255)
                                                       cGCCCGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCctgct                                                                                           <  2:5875/150‑1 (MQ=255)
                                                        gCCCGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCctgctg                                                                                          <  1:154727/150‑1 (MQ=255)
                                                          ccGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTa                                                                                                  >  2:1252474/1‑140 (MQ=255)
                                                          ccGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTa                                                                                                  <  1:1252474/140‑1 (MQ=255)
                                                          ccGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGt                                                                                        >  2:1316321/1‑150 (MQ=255)
                                                            gcgGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAg                                                                                      <  2:63408/150‑1 (MQ=255)
                                                                  tgctgcGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTc                                                                                <  2:125898/150‑1 (MQ=255)
                                                                     tgcGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCGTTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAgg                                                                                     >  1:1062400/1‑142 (MQ=255)
                                                                     tgcGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAgg                                                                                     <  2:1062400/142‑1 (MQ=255)
                                                                              tGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTCCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTACCCTACCCTATCCTGCTGGTAGGCCGTCTGCAGGTCTTCg                                                                    >  2:472800/1‑150 (MQ=255)
                                                                                 aGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAAc                                                                 <  1:1374542/150‑1 (MQ=255)
                                                                                 aGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAAc                                                                 <  1:182431/150‑1 (MQ=255)
                                                                                  gCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACa                                                                <  2:1168359/150‑1 (MQ=255)
                                                                                     aGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGc                                                                                                <  1:405242/115‑1 (MQ=255)
                                                                                     aGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGc                                                                                                >  2:405242/1‑115 (MQ=255)
                                                                                            aCGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCa                                                      >  1:1092520/1‑150 (MQ=255)
                                                                                              gTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCAt                                                                                                                                                     <  1:588035/53‑1 (MQ=255)
                                                                                              gTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCAt                                                                                                                                                     >  2:588035/1‑53 (MQ=255)
                                                                                              gTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAg                                                    >  2:426459/1‑150 (MQ=255)
                                                                                                 cAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGccc                                                                                                         <  1:1216456/94‑1 (MQ=255)
                                                                                                 cAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGccc                                                                                                         >  2:1216456/1‑94 (MQ=255)
                                                                                                  aGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGacac                                                <  2:510587/150‑1 (MQ=255)
                                                                                                               tCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATccaccacc                                   <  2:551643/150‑1 (MQ=255)
                                                                                                                 gCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAg                                 >  2:1165914/1‑150 (MQ=255)
                                                                                                                  cTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCGCCACCAgg                                <  2:376143/150‑1 (MQ=255)
                                                                                                                     cgcgcATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGt                               >  1:546880/1‑148 (MQ=255)
                                                                                                                     cgcgcATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGt                               <  2:546880/148‑1 (MQ=255)
                                                                                                                        gcATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGc                          <  2:417359/150‑1 (MQ=255)
                                                                                                                        gcATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGc                          >  1:1064807/1‑150 (MQ=255)
                                                                                                                         cATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCa                         >  1:890762/1‑150 (MQ=255)
                                                                                                                          aTAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCg                                                                                  >  2:213587/1‑92 (MQ=255)
                                                                                                                          aTAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCg                                                                                  <  1:213587/92‑1 (MQ=255)
                                                                                                                          aTAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCAc                        <  2:250009/150‑1 (MQ=255)
                                                                                                                             gTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGc                     <  2:794615/150‑1 (MQ=255)
                                                                                                                               catcaAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCa                   <  2:1092520/150‑1 (MQ=255)
                                                                                                                                atcaAGCTGTCGCGCCCATTGGCGCAGGGGTGCAGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGATGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAg                  >  1:473535/1‑150 (MQ=255)
                                                                                                                                 tcaAGCTGTCACGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGAt                                                            <  2:975949/107‑1 (MQ=255)
                                                                                                                                 tcaAGCTGTCACGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGAt                                                            >  1:975949/1‑107 (MQ=255)
                                                                                                                                    aGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGAGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAAc              >  1:456997/1‑150 (MQ=255)
                                                                                                                                     gCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCctgc                                                                                            <  2:1171118/71‑1 (MQ=255)
                                                                                                                                     gCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCctgc                                                                                            >  1:1171118/1‑71 (MQ=255)
                                                                                                                                        cctgcgcCAATGGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAgg                                >  1:478882/4‑128 (MQ=255)
                                                                                                                                        cctgcgcCAATGGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAgg                                <  2:478882/125‑1 (MQ=255)
                                                                                                                                        cctgcgcCAATGGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGAt                            >  1:1281265/4‑132 (MQ=255)
                                                                                                                                        cctgcgcCAATGGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGAt                            <  2:1281265/129‑1 (MQ=255)
                                                                                                                                         tCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCGGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCg         >  1:455612/1‑150 (MQ=255)
                                                                                                                                         tCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGccc         <  2:802938/150‑2 (MQ=255)
                                                                                                                                            cgcCCATTGGCGCAGGTGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTACGCAGTTTCCAGGTCTTCAAACGCTATTCGCGAAGACACCTCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTg      <  1:174626/150‑1 (MQ=255)
                                                                                                                                             gcCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGc     >  1:1183411/1‑150 (MQ=255)
                                                                                                                                                cATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGa  <  1:126347/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
GCATGCGCCGCCTGGCGCAGGTTGCATTCGCCTTGCGGGTTATCGCAAGGCAGCTCGCCCGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGA  >  NC_002947/4783451‑4783744

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