Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,652,960 T→G 5.4% intergenic (‑61/+100) PP_4961 ← / ← PP_4962 lipoprotein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,652,9600TG 5.4% 194.1 / 2.9 74intergenic (‑61/+100)PP_4961/PP_4962lipoprotein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (36/31);  new base G (2/2);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.76e-01

GTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTA  >  NC_002947/5652813‑5653101
                                                                                                                                                    |                                                                                                                                             
gTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTa                                                                                                                                              <  2:393940/149‑1 (MQ=255)
   gTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGCGCATCCTGCCAGTGGCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGTTTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACa                                                                                                                                            <  2:282126/148‑1 (MQ=255)
   gTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACa                                                                                                                                            <  2:75092/148‑1 (MQ=255)
   gTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAg                                                                                                                                           >  4:95436/1‑149 (MQ=255)
      tGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGc                                                                                                                                        >  7:342102/1‑149 (MQ=255)
      tGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGCGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTg                                                                                                                                         >  6:297216/1‑148 (MQ=255)
       gCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTg                                                                                                                                         >  7:220010/1‑147 (MQ=255)
        ccGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCATGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCa                                                                                                                                       <  3:31644/148‑1 (MQ=255)
        ccGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCa                                                                                                                                       <  7:85539/148‑1 (MQ=255)
            ggCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGtctc                                                                                                                                  <  4:8116/149‑1 (MQ=255)
                ggTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCACCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCtgt                                                                                                                               <  1:278301/149‑1 (MQ=255)
                 gTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTgg                                                                                                                             <  7:70832/149‑1 (MQ=255)
                      gcCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGAGACTGCACTGTACCTGTGGGAGCg                                                                                                                        <  8:112879/149‑1 (MQ=25)
                               agcagTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGccc                                                                                                               >  8:39161/1‑149 (MQ=255)
                                gcagTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCcg                                                                                                              >  5:199905/1‑149 (MQ=255)
                                 cagTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGc                                                                                                                 <  3:73487/145‑1 (MQ=255)
                                 cagTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGc                                                                                                                 >  4:73487/1‑145 (MQ=255)
                                  agTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGt                                                                                                                                     >  4:95663/1‑124 (MQ=255)
                                  agTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGt                                                                                                                                     <  3:95663/124‑1 (MQ=255)
                                       gCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTgg                                                                                                                             <  7:50733/127‑1 (MQ=255)
                                       gCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTgg                                                                                                                             >  8:50733/1‑127 (MQ=255)
                                       gCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAggg                                                                                                       >  7:223909/1‑149 (MQ=255)
                                        cATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGc                                                                                                      >  8:142272/1‑149 (MQ=255)
                                            cTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGca                                                                                                  >  6:288572/1‑149 (MQ=255)
                                             tGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGca                                                                                                  >  5:307628/1‑148 (MQ=255)
                                                cAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCGGTGCGGCCCTTCGCGGGCTCGCCCGCTCCCACAGAGACTGCACTGTACCTGTGGGAGCGGGCGAGCCCGAGAAGGGCTGCACAGc                                                                                              <  6:4517/149‑1 (MQ=14)
                                                 aGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGGGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCg                                                                                             <  3:95436/149‑1 (MQ=255)
                                                     tCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGccc                                                                                         >  7:367341/1‑149 (MQ=255)
                                                      ccccAGGGCCGCCACGGCAAAACCCCGTTCCATTCCTTTTCCTCCATTGGGGCCGGCGGCGCCGCCCTTCCCGGCCGCCCCCGCCCCCCCCGGTCCAGTCCCGGCTCTGTGGGCGGGGGCGAGCCCGCGAAGGGCCGCACAGCGGcccc                                                                                        <  8:342102/147‑1 (MQ=255)
                                                       gCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCct                                                                                       <  6:307628/149‑1 (MQ=255)
                                                        ccAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCctt                                                                                      >  5:293744/1‑149 (MQ=255)
                                                         cAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACCGTGCAGGCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCcttc                                                                                     <  5:146477/149‑1 (MQ=255)
                                                              ccGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCaa                                                                                <  6:293744/149‑1 (MQ=255)
                                                                tacccGGCCAAACGGGGTTTCATTCCGTTCCCTCGTTTGGGGGCTGCGGGGCTGCCGTTCTCGGGCTGGCCCGGTCCCACTGTAACTTTCCCGTCGGTGTGGGAGGGGGCGTGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGa                                                                              <  2:8702/145‑1 (MQ=255)
                                                                          aaaGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTGGTCAAGACCACAAATc                                                                     >  4:258592/1‑148 (MQ=255)
                                                                          aaaGCGATTTAATACAGTTTCCTCGATTGGGGGCTGGTGTCCAGCCCTTCTCGGGCGCCCCCGCACCCACTGGTACCGGCCAGTCTCTGTGGTAGCGGGCTAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCa                                                                    <  7:142272/149‑1 (MQ=25)
                                                                              cgcTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGt                                                                <  5:288572/149‑1 (MQ=255)
                                                                               gcTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGtt                                                               <  6:199905/149‑1 (MQ=255)
                                                                                   tcattcaGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAgg                                                           >  3:306604/1‑149 (MQ=255)
                                                                                    cattcaGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAACGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGt                                                          >  6:379772/1‑149 (MQ=255)
                                                                                             tCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTAttt                                                 >  3:187654/1‑148 (MQ=255)
                                                                                              ccTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCGTCTTCAAGACGACAAATCAACGTTCAGGTTGCCTATTAc                                                >  5:228518/1‑149 (MQ=255)
                                                                                                tCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAccc                                              >  3:7547/1‑149 (MQ=255)
                                                                                                    tttGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGAt                                           >  2:1299/1‑148 (MQ=255)
                                                                                                    tttGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGGGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATc                                          <  8:220010/149‑1 (MQ=255)
                                                                                                       ggggCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGGCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAg                                       <  8:367341/149‑1 (MQ=255)
                                                                                                        gggCTGCTGTGCAGCCTTCCTCGGGCGCGCCCGCTCCCACAGGAACAGGCCAGTCTCTGTGGTTGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAg                                       <  7:49821/148‑1 (MQ=255)
                                                                                                        gggCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACCGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAgg                                      <  4:256457/149‑1 (MQ=255)
                                                                                                        gggCTGCTGTGCAGCCCGTGTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAg                                       <  5:297216/148‑1 (MQ=255)
                                                                                                        ggcctggtgTGCTGCCCTTGGGGGGCCGGCCGGCCCCCACCGTTCCTGTGCAGTCTCGGGGGCTGGGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAAGGTTCAGGTTGCCTATTACCCGATCAAgg                                      <  4:262536/146‑1 (MQ=255)
                                                                                                         gcctgcggGGCCGCCCTTCTCGGGCGCCCCCGCCCCCACCGCTACAGTGCCGCCTCTGGGGGCGGGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAggg                                     <  6:228518/147‑1 (MQ=255)
                                                                                                           ctgctgTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGc                                    >  3:41351/1‑148 (MQ=255)
                                                                                                               tgtgCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAg                                       <  3:348117/141‑1 (MQ=255)
                                                                                                               tgtgCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAg                                       >  4:348117/1‑141 (MQ=255)
                                                                                                                 tgctgcCCTTCTCGGGCTCGCCCGCTCCCACAGTTACTGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCgg                             <  5:322404/149‑1 (MQ=255)
                                                                                                                  gCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGgc                            >  8:123930/1‑149 (MQ=255)
                                                                                                                   cAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGgcg                           >  1:244203/1‑149 (MQ=255)
                                                                                                                    aGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGACGGGCCGCCCAGCGGCCCCGTCTTCAAGACCACAAAGCAACGTTCAGGTTGGCGATTTCCCGCTCAAGGGCAACCCGgcgc                          >  5:233400/1‑149 (MQ=38)
                                                                                                                    aGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGAGCAAGGGCAACCCGgcgc                          >  2:334216/1‑149 (MQ=255)
                                                                                                                         cTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGAAGTCTCTGTGGGAGCGGGCGAGCACGCGCAGGGCCGCAAAGCGGCCCCGTATTCAAGACCACAAAGCAACGTTCAGGTTGCCTATTTCGCGATCAAGGGCAACCCGGCGCCGTg                      >  4:337751/1‑149 (MQ=255)
                                                                                                                          ttCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTgc                     >  5:290206/1‑149 (MQ=255)
                                                                                                                                ggCGCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCTAGGGCAACCCGGCGCCTTGCGCtggt               <  1:11631/149‑1 (MQ=255)
                                                                                                                                ggCGCGCCCGCTCCCACAGGAACAGTGCAGTCTCTGGGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCtggt               <  8:223909/149‑1 (MQ=255)
                                                                                                                                 gCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACCGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTTTTACCCGATCAAGGGCAACCCGGCGCCTTGCGCtggtg              >  3:51105/1‑149 (MQ=255)
                                                                                                                                 gCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCtt                       <  3:225340/140‑1 (MQ=255)
                                                                                                                                 gCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCtt                       >  4:225340/1‑140 (MQ=255)
                                                                                                                                  cTCGCCCGCTCCCACAGAGACTGCACTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCtggtgg             >  4:13148/1‑149 (MQ=25)
                                                                                                                                     gCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCACGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCt          >  1:281917/1‑149 (MQ=255)
                                                                                                                                      cccGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCtt         >  5:245833/1‑149 (MQ=255)
                                                                                                                                      cccGCTCCCACAGAGACTGCACTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCtt         >  4:386470/1‑149 (MQ=25)
                                                                                                                                       ccGCTCCCTCTGCTACAGTGCTGGGTCTGTGGGTGGGGGCGAGCCCGCGAAGGGCCGCACAGCTGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTtc        <  4:51105/149‑1 (MQ=255)
                                                                                                                                             ccACCGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTa  <  4:41351/149‑1 (MQ=255)
                                                                                                                                             ccACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTa  >  4:246419/1‑149 (MQ=255)
                                                                                                                                                   gTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTtcg       >  4:311665/1‑138 (MQ=255)
                                                                                                                                                   ctACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTtcg       <  3:311665/137‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
GTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTA  >  NC_002947/5652813‑5653101

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