Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,783,598 T→G 5.3% Q1254P (CAA→CCA)  pvdJ ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,783,5980TG 5.3% 233.9 / 4.1 76Q1254P (CAA→CCA) pvdJnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base T (34/38);  new base G (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 = 7.51e-01

GCATGCGCCGCCTGGCGCAGGTTGCATTCGCCTTGCGGGTTATCGCAAGGCAGCTCGCCCGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACC  >  NC_002947/4783451‑4783746
                                                                                                                                                   |                                                                                                                                                    
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                                                                                                                                     <  2:1163362/145‑1 (MQ=255)
                    gTTGCATTCGCCTTGCGGGTTATCGCAAGGCAGCTCGCCCGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGAc                                                                                                                                     >  1:1163362/1‑145 (MQ=255)
                    gTTGCATTCGCCTTGCGGGTTATCGCAAGGCAGCTCGCCCGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGa                                                                                                                                >  1:376143/1‑150 (MQ=255)
                         aTTCGCCTTGCGGGTTATCGCAAGGCAGCTCGCCCGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTg                                                                                                                                                     <  1:537532/124‑1 (MQ=255)
                         aTTCGCCTTGCGGGTTATCGCAAGGCAGCTCGCCCGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTg                                                                                                                                                     >  2:537532/1‑124 (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                                                                                                                                <  2:560117/128‑1 (MQ=255)
                                          tCGCAAGGCAGCTCGCCCGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGa                                                                                                                                >  1:560117/1‑128 (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                                                                                                  >  1:1359854/1‑146 (MQ=255)
                                                      tCGCCCGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGc                                                                                                  <  2:1359854/146‑1 (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:182431/150‑1 (MQ=255)
                                                                                 aGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAAc                                                                   <  1:1374542/150‑1 (MQ=255)
                                                                                  gCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACa                                                                  <  2:1168359/150‑1 (MQ=255)
                                                                                     aGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGc                                                                                                  >  2:405242/1‑115 (MQ=255)
                                                                                     aGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGc                                                                                                  <  1:405242/115‑1 (MQ=255)
                                                                                            aCGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCa                                                        >  1:1092520/1‑150 (MQ=255)
                                                                                              gTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAg                                                      >  2:426459/1‑150 (MQ=255)
                                                                                                 cAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGccc                                                                                                           >  2:1216456/1‑94 (MQ=255)
                                                                                                 cAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGccc                                                                                                           <  1:1216456/94‑1 (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                                                                                              >  1:1171118/1‑71 (MQ=255)
                                                                                                                                     gCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCctgc                                                                                              <  2:1171118/71‑1 (MQ=255)
                                                                                                                                        cctgcgcCAATGGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAgg                                  <  2:478882/125‑1 (MQ=255)
                                                                                                                                        cctgcgcCAATGGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAgg                                  >  1:478882/4‑128 (MQ=255)
                                                                                                                                        cctgcgcCAATGGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGAt                              <  2:1281265/129‑1 (MQ=255)
                                                                                                                                        cctgcgcCAATGGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGAt                              >  1:1281265/4‑132 (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)
                                                                                                                                              ccaaTTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCt      >  1:1421003/4‑150 (MQ=255)
                                                                                                                                                cATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGa    <  1:126347/150‑1 (MQ=255)
                                                                                                                                                  ttGGCGCAGGGGTGCCGACCTGGACGGCATGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACGGCTGTTGCTggcg  <  1:332611/150‑4 (MQ=255)
                                                                                                                                                   tGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGt                                 <  1:322624/118‑1 (MQ=255)
                                                                                                                                                   tGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGt                                 >  2:322624/1‑118 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCATGCGCCGCCTGGCGCAGGTTGCATTCGCCTTGCGGGTTATCGCAAGGCAGCTCGCCCGCGGGCTGCTGCGCGGCCTGGAGCCAGTAATCACGTTCAGCCTGCAGGTTTTCGCTGCGCGCATAGTCATCAAGCTGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACC  >  NC_002947/4783451‑4783746

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