Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 427,133 A→G 6.6% Y222Y (TAT→TAC)  PP_0350 ← ferrichrome‑iron receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947427,1330AG 6.6% 174.3 / 2.5 60Y222Y (TAT→TAC) PP_0350ferrichrome‑iron receptor
Reads supporting (aligned to +/- strand):  ref base A (28/28);  new base G (2/2);  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 = 5.03e-01

CTCGATTTCGTGAGCGCTTTCGGCTTGCGGACGGCGGCTGACCATGTCCAGCAGGCCGCCAGGTGGAGTCTGGCCGTAAATCGAAGAGGCCGGGCCGCGCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCCGCGCACGCGCATCCAGTCATACCGCGTATCGCTGCCGAAGCTGGCAGAGACGATGCCCGGCATATAGCGCACTGCATCATCCAGGCTGTGTACAC  >  NC_002947/426986‑427281
                                                                                                                                                   |                                                                                                                                                    
cTCGATTTCGTGAGCGCTTTCGGCTTGCGGACGGCGGCTGACCATGTCCAGCAGGCCGCCAGGTGGAGTCTGGCCGTAAATCGAAGAGGCCGGGCCGCGCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCa                                                                                                                                                      >  6:253698/1‑148 (MQ=255)
                cTTTCGGCTGGCGGACGTCGGCTGACCTTGTCCAGCAGGCCTCCAGGTGGAGTCTGGCAGTAAATCGCAGAGGCCGGGCCGCGCAGCAGGGCGAGGCGGTCGAGGGTCCTGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAg                                                                                                                                     <  5:260827/149‑1 (MQ=255)
                 tttCGGCTTGCGGACGGCGGCTGACCATGTCCAGCAGGCCGCCAGGTGGAGTCTGGCCGTAAATCGAAGAGGCCGGGCCGCGCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAgg                                                                                                                                    >  7:100062/1‑149 (MQ=255)
                   tCGGCTTGCGGACGGCGGCTGACCATGTCCAGCAGGCCGCCAGGTGGAGTCTGGCCGTAAATCGAAGAGGCCGGGCCGCGCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAggg                                                                                                                                   <  1:139781/148‑1 (MQ=255)
                      gCTTGCGGACGGCGGCTGACCATGTCCAGCAGGCCGCCAGGTGGAGTCTGGCCGTAAATCGAAGAGGCCGGGCCGCGCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATaca                                                                                                                                                  >  4:39876/1‑130 (MQ=255)
                      gCTTGCGGACGGCGGCTGACCATGTCCAGCAGGCCGCCAGGTGGAGTCTGGCCGTAAATCGAAGAGGCCGGGCCGCGCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATaca                                                                                                                                                  <  3:39876/130‑1 (MQ=255)
                       cTTGCGGACGGCGGCTGACCATGTCCAGCAGGCCGCCAGGTGGAGTCTGGCCGTAAATCGAAGAGGCCGGGCCGCGCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAgg                                                                                                                              >  1:235808/1‑149 (MQ=255)
                       cTTGCGGACGGCGGCTGACCATGTCCAGCAGGCCGCCAGGTGGAGTCTGGCCGTAAATCGAAGAGGCCGGGCCGCGCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAgg                                                                                                                              <  6:122120/149‑1 (MQ=255)
                            ggACGGCGGCTGACCATGTCCAGCAGGCCGCCAGGTGGAGTCTGGCCGTAAATCGAAGAGGCCGGGCCGCGCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATc                                                                                                                         <  1:268133/149‑1 (MQ=255)
                              aCGGCGGCTGACCATGTCCAGCAGGCCGCCAGGTGGAGTCTGGCCGTAAATCGAAGAGGCCGGGCCGCGCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGa                                                                                                                       <  5:253698/149‑1 (MQ=255)
                                       gACCATGTCCAGCAGGCCGCCAGGTGGAGTCTGGCCGTAAATCGAAGAGGCCGGGCCGCGCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCGTACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCg                                                                                                              >  5:80902/1‑149 (MQ=255)
                                       gACCATGTCCAGCAGGCCGCCAGGTGGAGTCTGGCCGTAAATCGAAGAGGCCGGGCCGCGCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCGTACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCg                                                                                                              >  3:143090/1‑149 (MQ=255)
                                       gACCATGTCCAGCAGGCCGCCAGGTGGAGTCTGGCCGTAAATCGAAGAGGCCGGGCCGCGCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCg                                                                                                              >  2:228549/1‑149 (MQ=255)
                                        aCCATGTCCAGCAGGCCGCCAGGTGGAGTCTGGCCGTAAATCGAAGAGGCCGGGCCGCGCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGt                                                                                                             <  8:140021/149‑1 (MQ=255)
                                         ccATGTCCAGCAGGCCGCCAGGTGGAGTCTGGCCGTAAATCGAAGAGGCCGGGCCGCGCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGtt                                                                                                            >  3:42361/1‑149 (MQ=255)
                                         ccATGTCCAGCAGGCCGCCAGGTGGAGTCTGGCCGTAAATCGAAGAGGCCGGGCCGCGCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGtt                                                                                                            <  8:100062/149‑1 (MQ=255)
                                          cATGTCCAGCAGGCCGCCAGGTGGAGTCTGGCCGTAAATCGAAGAGGCCGGGCCGCGCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTg                                                                                                           <  2:84506/149‑1 (MQ=255)
                                          cATGTCCAGCAGGCCGCCAGGTGGAGTCTGGCCGTAAATCGAAGAGGCCGGGCCGCGCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTg                                                                                                           >  7:242199/1‑149 (MQ=255)
                                              tCCAGCAGGCCGCCAGGTGGAGTCTGGCCGTAAATCGAAGAGGCCGGGCCGCGCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTc                                                                                                       >  5:183388/1‑149 (MQ=255)
                                                cagcagGCCGCCAGGTGGAGTCTGGCCGTAAATCGAAGAGGCCGGGCCGCGCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCg                                                                                                      <  4:42361/148‑1 (MQ=255)
                                                         gccAGGTGGAGTCTGGCCGTAACTCGAAGAGGCCGGGCCGCGCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCcgcgca                                                                                            <  4:86439/149‑1 (MQ=255)
                                                         gccAGGTGGAGTCTGGCCGTAAATCGAAGAGGCCGGGCCGCGCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCcgcgca                                                                                            <  4:151267/149‑1 (MQ=255)
                                                                 gAGTCTGGCCGTAAATCGAAGAGGCCGGGCCGCGCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCGTACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCCGCGCACGCGCATc                                                                                    <  4:143090/149‑1 (MQ=255)
                                                                 gAGTCTGGCCGTAAATCGAAGAGGCCGGGCCGCGCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCGTACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCCGCGCACGCGCATc                                                                                    <  6:80902/149‑1 (MQ=255)
                                                                   gTCTGGCCGTAAATCGAAGAGGCCGGGCCGCGCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCcgcgcacgcgc                                                                                       <  6:46274/144‑1 (MQ=255)
                                                                   gTCTGGCCGTAAATCGAAGAGGCCGGGCCGCGCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCcgcgcacgcgc                                                                                       >  5:46274/1‑144 (MQ=255)
                                                                   gTCTGGCCGTAAATCGAAGAGGCCGGGCCGCGCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCCGCGCACGCGCATCCa                                                                                  <  1:51529/149‑1 (MQ=255)
                                                                      tGGCCGTAAATCGAAGAGGCCGGGCCGCGCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCCGCGCACGCGCATCCAGt                                                                                >  7:41953/1‑148 (MQ=255)
                                                                      tGGCCGTAAATCGAAGAGGCCGGGCCGCGCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCCGCGCACGCGCATCCAGt                                                                                >  7:26035/1‑148 (MQ=255)
                                                                         ccGTAAATCGAAGAGGCCGGGCCGCGCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCCGCGCACGCGCATCCAGTCATa                                                                            <  6:233388/149‑1 (MQ=255)
                                                                             aaaTCGAAGAGGCCGGGCCGCGCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCCGCGCACGCGCATCCAGTCATACcgc                                                                        >  3:174036/1‑149 (MQ=255)
                                                                              aaTCGAAGAGGCCGGGCCGCGCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCCGCGCACGCGCATCCAGTCATACcgcg                                                                       >  6:279746/1‑149 (MQ=255)
                                                                                    agagGCCGGGCCGCGCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCCGCGCACGCGCATCCAGTCATACCGCGTATCGc                                                                 <  8:41953/149‑1 (MQ=255)
                                                                                     gagGCCGGGCCGCGCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCCGCGCACGCGCATCCAGTCATACCGCGTATCGCt                                                                >  6:32799/1‑149 (MQ=255)
                                                                                     gagGCCGGGCCGCGCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCCGCGCACGCGCATCCAGTCATACCGCGTATCGCt                                                                <  8:242199/149‑1 (MQ=255)
                                                                                       ggccgggccgCGCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCCGCGCACGCGCATCCAGTCATACCGCGTATCGCTGc                                                              <  1:228549/149‑1 (MQ=255)
                                                                                        gccgggccgCGCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCCGCGCACGCGCATCCAGTCATACCGCGTATCGCTGcc                                                             >  1:7219/1‑149 (MQ=255)
                                                                                          cgggccgCGCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCCGCGCACGCGCATCCAGTCATACCGCGTATCGCTGCCg                                                            >  4:223739/1‑148 (MQ=255)
                                                                                               cgcgCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCCGCGCACGCGCATCCAGTCATACCGCGTATCGCTGCCGAAGCTg                                                      >  1:64752/1‑149 (MQ=255)
                                                                                               cgcgCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCCGCGCACGCGCATCCAGTCATACCGCGTATCGCTGCCGAAGCTg                                                      <  6:183388/149‑1 (MQ=255)
                                                                                                gcgcAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCCGCGCACGCGCATCCAGTCATACCGCGTATCGCTGCCGAAGCTgg                                                     <  2:7219/149‑1 (MQ=255)
                                                                                                gcgcAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCCGCGCACGCGCATCCAGGCATACCGCGTATCGCTGCCGAAGCTgg                                                     <  2:235808/149‑1 (MQ=255)
                                                                                                  gcagcaGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCCGCGCACGCGCATCCAGTCATACCGCGTATCGCTGCCGAAGCTgg                                                     >  8:295745/1‑147 (MQ=255)
                                                                                                    agcagGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCCGCGCACGCGCATCCAGTCATACCGCGTATCGCTGCCGAAGCTGGCa                                                   >  1:207253/1‑147 (MQ=255)
                                                                                                                  gTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCCGCGCACGCGCATCCAGTCATACCGCGTATCGCTGCCGAAGCTGGCAGAGACGATGCCCGGCa                                   <  8:26035/149‑1 (MQ=255)
                                                                                                                         ttCCAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCCGCGCACGCGCATCCAGTCATACCGCGTATCGCTGCCGAAGCTGGCAGAGACGATGCCCGGCATATAgcg                            >  5:252296/1‑149 (MQ=255)
                                                                                                                           ccAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCTGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCCACGCACGCGCATCCAGTCATACCGCGTATCGCTGCCGAAGCTGGCAGAGACGATGCCCGGCATATAgcgc                           >  5:26361/1‑148 (MQ=255)
                                                                                                                           ccAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCTTTGGATCGAAGCCGCGCACGCGCATCCAGTCATACCGCGTATCGCTGCCGAAGCTGGCAGAGACGATGCCCGGCATATAgcgc                           >  1:154964/1‑148 (MQ=255)
                                                                                                                           ccAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCCGCGCACGCGCATCCAGTCATACCGCGTATCGCTGCCGAAGCTGGCAGAGACGATGCCCGGCATATAgcgc                           <  3:223739/148‑1 (MQ=255)
                                                                                                                             aGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCCGCGCACGCGCATCCAGTCATACCGCGTATCGCTGCCGAAGCTGGCAGAGACGATGCCCGGCATATAGCGCACt                        >  1:180244/1‑149 (MQ=255)
                                                                                                                              ggTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCCGCGCACGCGCATCCAGTCATACCGCGTATCGCTGCCGAAGCTGGCAGAGACGATGCCCGGCATATAGCGCACTg                       >  6:55957/1‑149 (MQ=255)
                                                                                                                                 ttCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCCGCGCACGCGCATCCAGTCATACCGCGTATCGCTGCCGAAGCTGGCAGAGACGATGCCCGGCATATAGCGCACTGca                     >  7:178597/1‑148 (MQ=255)
                                                                                                                                    ggCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCCGCGCACGCGCATCCAGTCATACCGCGTATCGCTGCCGAAGCTGGCAGAGACGATGCCCGGCATATAGCGCACTGcatcat                 <  2:64752/149‑1 (MQ=255)
                                                                                                                                    ggCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCCGCGCACGCGCATCCAGTCATACCGCGTATCGCTGCCGAAGCTGGCAGAGACGATGCCCGGCATATAGCGCACTGcatcat                 >  6:117096/1‑149 (MQ=255)
                                                                                                                                       tttCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCcgc                                                                                               <  1:258230/68‑1 (MQ=255)
                                                                                                                                       tttCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCcgc                                                                                               >  2:258230/1‑68 (MQ=255)
                                                                                                                                       tttCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCCGCGCACGCGCATCCAGTCATACCGCGTATCGCTGCCGAAGCTGGCAGAGACGATGCCCGGCata                                 <  7:256386/130‑1 (MQ=255)
                                                                                                                                       tttCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCCGCGCACGCGCATCCAGTCATACCGCGTATCGCTGCCGAAGCTGGCAGAGACGATGCCCGGCata                                 >  8:256386/1‑130 (MQ=255)
                                                                                                                                          cGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCCGCGCACGCGCATCCAGTCATACCGCGTATCGCTGCCGAAGCTGGCAGAGACGATGCCCGGCATATAGCGCACTGCATCATCCAgg            <  4:174036/148‑1 (MQ=255)
                                                                                                                                            ggTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCCGCGCACGCGCATCCAGTCATACCGCGTATCGCTGCCGAAGCTGGCAGAGACGATGCCCGGCATATAGCGCACTGCATCATCCAGGCtg         <  5:32799/149‑1 (MQ=255)
                                                                                                                                            ggTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCCGCGCACGCGCATCCAGTCATACCGCGTATCGCTGCCGAAGCTGGCAGAGACGATGCCCGGCATATAGCGCACTGCATCATCCAGGCtg         >  4:59582/1‑149 (MQ=255)
                                                                                                                                                   aTACACGCCGCGGGGCAGGGGTAGGCCATAGAGGAACTGCGTTGGCTCGAAGCCGCGCACGCGCATCCAGTCATACCGCGTATCGCTGCCGAAGCTGGCAGAGACGATGCCCGGCATATAGCGCACTGCATCATCCAGGCTGTGTacac  >  2:263240/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CTCGATTTCGTGAGCGCTTTCGGCTTGCGGACGGCGGCTGACCATGTCCAGCAGGCCGCCAGGTGGAGTCTGGCCGTAAATCGAAGAGGCCGGGCCGCGCAGCAGGGCGAGGCGGTCGAGGTTCCAGGTTTCGGCTTTCGGGTTGGCATACACGCCGCGGGGCAGGGGTAGGCCATCGAGGAACTGCGTTGGCTCGAAGCCGCGCACGCGCATCCAGTCATACCGCGTATCGCTGCCGAAGCTGGCAGAGACGATGCCCGGCATATAGCGCACTGCATCATCCAGGCTGTGTACAC  >  NC_002947/426986‑427281

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