Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,982,954 G→T 5.8% intergenic (+113/‑33) PP_3510 → / → ilvE hypothetical protein/branched‑chain‑amino acid aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,982,9540GT 5.8% 210.3 / 3.1 69intergenic (+113/‑33)PP_3510/ilvEhypothetical protein/branched‑chain‑amino acid aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (35/30);  new base T (2/2);  total (37/32)
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.97e-01

ACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGACTACATCAAGACCGACAAACGCTACCTGTCCGTATGGCGCAATGGCGAGTGGGACAAGGGCACCCTGACCGAA  >  NC_002947/3982813‑3983100
                                                                                                                                             |                                                                                                                                                  
aCTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCgg                                                                                                                                             <  5:200681/149‑1 (MQ=255)
  tGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCggg                                                                                                                                            >  8:227748/1‑148 (MQ=255)
      gAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACtt                                                                                                                                       <  1:185608/149‑1 (MQ=255)
      gAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACtt                                                                                                                                       >  8:30215/1‑149 (MQ=255)
                 cAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATg                                                                                                                            <  4:325698/149‑1 (MQ=255)
                  agcggCATTGACCCTGCGGGTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCAGATCCGCCAAGGATGTATTAAGGTGGCATCGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGa                                                                                                                           <  6:22453/146‑1 (MQ=255)
                   ggggCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCATCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAg                                                                                                                          <  5:320556/149‑1 (MQ=255)
                    gggCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCCCTTGCAAGACAATGAgg                                                                                                                         >  8:203699/1‑149 (MQ=255)
                    gggCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAgg                                                                                                                         >  7:288111/1‑149 (MQ=255)
                    gggCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAAATCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAgg                                                                                                                         <  2:278550/149‑1 (MQ=255)
                      gCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGgt                                                                                                                       >  6:135768/1‑149 (MQ=255)
                                tGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGc                                                                                                              >  7:268280/1‑148 (MQ=255)
                                         aCGATAGTCTCGCTTCGCACCCGATCAAATCGTGACTGGCCGTTCCCGATCCGCCAAGGATGGAGTTTGGTGGCATGGTGCGCTTATAACAGCAGCCCGCGCGCTCGGGCACTGGCAAGCCAGTGAGGGTGTCATGAGCACCggagg                                                                                                      >  8:144010/1‑143 (MQ=255)
                                         aCGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCa                                                                                                    >  4:341640/1‑149 (MQ=255)
                                            aTAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAg                                                                                                               >  2:34553/1‑135 (MQ=255)
                                            aTAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAg                                                                                                               <  1:34553/135‑1 (MQ=255)
                                                tATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTg                                                                                             <  5:135768/149‑1 (MQ=255)
                                                 aTCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTg                                                                                             <  5:261725/148‑1 (MQ=255)
                                                  tCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTgg                                                                                            <  4:364581/148‑1 (MQ=255)
                                                           aCCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTgg                                                                                  <  8:268280/149‑1 (MQ=255)
                                                           aCCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTgg                                                                                  >  3:26423/1‑149 (MQ=255)
                                                                 tCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTc                                                                             >  1:76417/1‑148 (MQ=255)
                                                                   aaTTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAAtt                                                                                              >  6:13255/1‑129 (MQ=255)
                                                                   aaTTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAAtt                                                                                              <  5:13255/129‑1 (MQ=255)
                                                                   aaTTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGa                                                                           <  5:123338/148‑1 (MQ=255)
                                                                    aTTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTgg                                                                                            <  1:314073/130‑1 (MQ=255)
                                                                    aTTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTgg                                                                                            >  2:314073/1‑130 (MQ=255)
                                                                     ttCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGACTa                                                                        <  5:313369/149‑1 (MQ=255)
                                                                     ttCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGACTa                                                                        >  2:185600/1‑149 (MQ=255)
                                                                       cGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGACTAc                                                                       <  7:30215/148‑1 (MQ=255)
                                                                        gTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACATGCTGGGTTTCGACTACAt                                                                     >  6:48758/1‑149 (MQ=255)
                                                                         tGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGTTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGACTACAt                                                                     >  6:364823/1‑148 (MQ=255)
                                                                           aCTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGCGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGCCTACATc                                                                    >  2:307009/1‑147 (MQ=255)
                                                                                      ccGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGACTACATCAAGACCGACAAAc                                                       <  5:48758/149‑1 (MQ=255)
                                                                                        gATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGACTACATCAAGACCGACAAACGc                                                     >  4:145839/1‑149 (MQ=255)
                                                                                        gATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGACTACATCAAGACCGACAAACGc                                                     >  1:25879/1‑149 (MQ=255)
                                                                                                 aGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGACTACATCAAGACCGACAAACGCTACCTGTcc                                            <  2:25879/149‑1 (MQ=255)
                                                                                                  ggATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGACTACATCAAGACCGACAAACGCTACCTGTCCg                                           <  5:364823/149‑1 (MQ=255)
                                                                                                       tATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGACTACATCAAGACCGACAAACGCTACCTGTCCGTATg                                       <  7:227748/148‑1 (MQ=255)
                                                                                                         ttAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGACTACATCAAGACCGACAAACGCTACCTGTCCGTATGgcg                                    >  1:208604/1‑149 (MQ=255)
                                                                                                         ttAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGACTACATCAAGACCGACAAACGCTACCTGTCCGTATGgcg                                    >  4:143594/1‑149 (MQ=255)
                                                                                                            aGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGACTACATCAAGACCGACAAACGCTACCTGTCCATATGGCGCaa                                 <  3:341640/149‑1 (MQ=255)
                                                                                                             ggTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGACTACATCAAGACCGACAAACGCTACCTGTCCGTATGGCGCAAt                                >  4:29537/1‑149 (MQ=255)
                                                                                                                   aTGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGACTACATCAAGACCGACAAACGCTAc                                                  <  3:237338/125‑1 (MQ=255)
                                                                                                                   aTGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGACTACATCAAGACCGACAAACGCTAc                                                  >  4:237338/1‑125 (MQ=255)
                                                                                                                    tGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGACTACATCAAGACCGACAAACGCTACCTGTCCGTATGGCGCAATGGCGAGt                         >  1:337995/1‑149 (MQ=255)
                                                                                                                        gcgcATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGACTACATCAAGACCGACAAACGCTAc                                                  >  5:210193/1‑120 (MQ=255)
                                                                                                                        gcgcATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGACTACATCAAGACCGACAAACGCTAc                                                  <  6:210193/120‑1 (MQ=255)
                                                                                                                        gcgcATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGAAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGACTACATCAAGAACGACAAACGCTCCATGTCCGTATGGCGCAATGGCGAGTGGGa                     >  2:331760/1‑149 (MQ=255)
                                                                                                                         cgcATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGACTACATCAAGACCGACAAACGCTACCTGTCCGTATGGCGCAATGGCGAGTGGGAc                    <  1:185600/149‑1 (MQ=255)
                                                                                                                          gcATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGACTACATCAAGACCGACAAACGCTACCTGTCCGTATGgc                                     >  8:74403/1‑131 (MQ=255)
                                                                                                                          gcATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGACTACATCAAGACCGACAAACGCTACCTGTCCGTATGgc                                     <  7:74403/131‑1 (MQ=255)
                                                                                                                            atatAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGACTACATCAAGACCGACAAACGCTACCt                                                >  2:116619/1‑118 (MQ=255)
                                                                                                                            atatAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGACTACATCAAGACCGACAAACGCTACCTGTCCGTATGGCGCAATGGCGAGTGGGACAAg                 <  8:288111/149‑1 (MQ=255)
                                                                                                                            atatAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGACTACATCAAGACCGACAAACGCTACCTGTCCGTATGGCGCAATGGCGAGTGGGACAAg                 <  3:143594/149‑1 (MQ=255)
                                                                                                                                aaCAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGACTACATCAAGACCGACAAACGCTACCTGTCCGTATGGCGCAATGGCGAGTGGGACAAGGGCa             >  5:360838/1‑149 (MQ=255)
                                                                                                                                aaCAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGACTACATCAAGACCGACAAACGCTACCTGTCCGTATGGCGCAATGGCGAGTGGGACAAGGGCa             >  6:182082/1‑149 (MQ=255)
                                                                                                                                   agcagcCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGACTACATCAAGACCGACAAACGCTACCTGTCCGTATGGCGCAATGGCGAGTGGGACAAGGGCAccc          >  5:365096/1‑149 (MQ=255)
                                                                                                                                   agcagcCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGACTACATCAAGACCGACAAACGCTACCTGTCCGTATGGCGCAATGGCGAGTGGGACAAGGGCAccc          >  1:196545/1‑149 (MQ=255)
                                                                                                                                    gcagcCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGATTCGACTACATCAAGACCGACAAACGCTACCTGTCCGTATGGCGCAATGGCGAGTGGGACAAGGGCACCCt         >  7:269671/1‑149 (MQ=255)
                                                                                                                                    agtgcCCGATCGCTCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGACTACATCAAGACCGACAAACGCTACCTGTCCGTATGGCGCAATGGCGAGTgg                       <  1:245210/132‑1 (MQ=255)
                                                                                                                                    agtgcCCGATCGCTCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGACTACATCAAGACCGACAAACGCTACCTGTCCGTATGGCGCAATGGCGAGTgg                       >  2:245210/4‑135 (MQ=255)
                                                                                                                                     gtgcCCGATCGCTCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGACTACATCAAGACCGACAAACGCTACCTGTCCGTATg                                       >  1:254805/3‑118 (MQ=255)
                                                                                                                                     gtgcCCGATCGCTCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGACTACATCAAGACCGACAAACGCTACCTGTCCGTATg                                       <  2:254805/116‑1 (MQ=255)
                                                                                                                                       gcCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGACTACATCAAGACCGACAAACGCTACCTGTCCGTATGGCGCAATGGCGAGTGGGACAAGGGCACCCTGAc      >  2:273540/1‑149 (MQ=255)
                                                                                                                                       gcCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGACTACATCAAGACCGACAAACGCTACCTGTCCGTATGGCGCAATGGCGAGTGGGACAAGGGCACCCTGAc      >  7:240552/1‑149 (MQ=255)
                                                                                                                                           gAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTACTTGGGACAAGCTGGGGTTCGACTACATCAAGCCCGACAAACGCTACCTGTCCGTATGGCGCAATGGCGAGTGGGACAAGGGCACCCTGACCGaa  <  8:98330/149‑1 (MQ=255)
                                                                                                                                           gAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGACTACATCAAGACCGACAAACGCTACCTGTCCGTATGGCGCAATGGCGAGTGGGACAAGGGCACCCTGACCGaa  <  3:145839/149‑1 (MQ=255)
                                                                                                                                           gAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGACTACATCAAGACCGACAAACGCTACCTGTCCGTATGGCGCAATGGCGAGTGGGACAAGGGCACCCTGACCGaa  <  4:26423/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
ACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCATGAGCAACGAAAGCATTAATTGGGACAAGCTGGGTTTCGACTACATCAAGACCGACAAACGCTACCTGTCCGTATGGCGCAATGGCGAGTGGGACAAGGGCACCCTGACCGAA  >  NC_002947/3982813‑3983100

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: