Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,087,488 C→G 6.7% intergenic (+137/‑289) PP_3597 → / → PP_3598 amino acid ABC transporter ATP‑binding protein/peptidase C26

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,087,4880CG 6.7% 220.5 / 2.8 70intergenic (+137/‑289)PP_3597/PP_3598amino acid ABC transporter ATP‑binding protein/peptidase C26
Reads supporting (aligned to +/- strand):  ref base C (33/32);  new base G (2/3);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCTGCACTGATTGAAAGCCTGTCACTGTGGGGGGTTTACCCGCGAAGAAGACACCGCCGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTCCGTCGTTTGACATATCGAAC  >  NC_002947/4087341‑4087633
                                                                                                                                                   |                                                                                                                                                 
tCCTGCACTGATTGAAAGCCTGTCACTGTGGGGGGTTTACCCGCGAAGAAGACACCGCCGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCt                                                                                                                                                  <  3:35232/149‑1 (MQ=255)
     cACTGATTGAAAGCCTGTCACTGTGGGGGGTTTACCCGCGAAGAAGACACCGCCGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGCGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAAt                                                                                                                                             >  6:175091/1‑149 (MQ=255)
     cACTGATTGAAAGCCTGTCACTGTGGGGGGTTTACCCGCGAAGAAGACACCGCCGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAAt                                                                                                                                             >  3:13972/1‑149 (MQ=255)
       cTGATTGAAAGCCTGTCACTGTGGGGGGTTTACCCGCGAAGAAGACACCGCCGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGa                                                                                                                                           <  3:175791/149‑1 (MQ=255)
           ttGAAAGCCTGTCACTGTGGGGGGTTTACCCGCGAAGAAGACACCGCCGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCAGCGGCGCAATGattag                                                                                                                                       <  5:35624/149‑5 (MQ=25)
           ttGAAAGCCTGTCACTGTGGGGGGTTTACCCGCGAAGAAGACACCGCCGTGGATGGCACGGGCTTCGCCCGAGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCa                                                                                                                                       <  1:119875/149‑1 (MQ=255)
              aaaGCCTGTCACTGTGGGGGGTTTACCCGCGAAGAAGACACCGCCGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCAGCGGCGCAATGattaggta                                                                                                                                    <  3:20907/149‑8 (MQ=25)
                   cTGTCACTGTGGGGGGTTTACCCGCGAAGAAGACACCGCCGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATc                                                                                                                               <  4:13972/149‑1 (MQ=255)
                   cTGTCACTGTGGGGGGTTTACCCGCGAAGAAGACACCGCCGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATc                                                                                                                               <  4:126778/149‑1 (MQ=255)
                        aCTGTGGGGGGTTTACCCGCGAAGAAGACACCGCCGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTaa                                                                                                                          <  3:119932/149‑1 (MQ=255)
                          tgtgGGGGGTTTACCCGCGAAGAAGACACCGCCGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCACCTGCATCGTTaaa                                                                                                                         >  7:328469/1‑148 (MQ=255)
                            tgGGGGGTTTACCCGCGAAGAAGACACCGCCGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAatat                                                                                                                      >  5:126426/1‑149 (MQ=255)
                                gggTTTACCCGCGAAGAAGACACCGCCGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATattatt                                                                                                                  <  7:143335/149‑1 (MQ=255)
                                  gTTTACCCGCGAAGAAGACACCGCCGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGa                                                                                                                >  8:335319/1‑149 (MQ=255)
                                      aCCCGCGAAGAAGACACCGCCGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCg                                                                                                             >  1:138558/1‑148 (MQ=255)
                                        ccGCGAAGAAGACACCGCCGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCCCTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGCCCCAACTGCATCGTTAAATATTATTGACCGCAg                                                                                                          >  5:293892/1‑149 (MQ=255)
                                        ccGCGAAGAAGACACCGCCGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCg                                                                                                             >  3:3064/1‑146 (MQ=255)
                                              agaagaCACCGCCGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTg                                                                                                    <  8:20786/149‑1 (MQ=255)
                                               gaagaCACCGCCGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGc                                                                                                   >  5:323910/1‑149 (MQ=255)
                                                aagaCACCGCCGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGAc                                                                                                        <  1:179843/143‑1 (MQ=255)
                                                aagaCACCGCCGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGAc                                                                                                        >  2:179843/1‑143 (MQ=255)
                                                aagaCACCGCCGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGc                                                                                                   >  3:298878/1‑148 (MQ=255)
                                                  gaCACCGCCGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCa                                                                                                 >  8:294104/1‑148 (MQ=255)
                                                   acacCGCCGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCTAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCaaa                                                                                               <  5:126316/149‑1 (MQ=255)
                                                   acacCGCCGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCaaa                                                                                               <  7:110157/149‑1 (MQ=255)
                                                               cgaacACGGGCGAAGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCca                                                                                   >  8:168079/5‑149 (MQ=255)
                                                                tGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCcac                                                                                  <  6:126426/149‑1 (MQ=255)
                                                                  gCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAg                                                                                <  6:323910/149‑1 (MQ=255)
                                                                   cACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAgg                                                                               <  7:294104/149‑1 (MQ=255)
                                                                    aCGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCAcgccgc                                                                                     >  1:251429/1‑142 (MQ=255)
                                                                    aCGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCAcgccgc                                                                                     <  2:251429/142‑1 (MQ=255)
                                                                     cGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACCCCGCCACAGGcc                                                                             <  8:328469/149‑1 (MQ=255)
                                                                        gCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGGTGCAATGACCCAACTGCATCGTTTAATATTATTGACCGCAGACCCTGCCAAACCCACGCTGCCACAGGCCCg                                                                           >  3:15793/1‑148 (MQ=255)
                                                                              cccGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAATTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGt                                                                    <  7:335319/149‑1 (MQ=255)
                                                                              cccGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGt                                                                    >  2:374542/1‑149 (MQ=255)
                                                                               ccGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGtt                                                                   <  2:138558/149‑1 (MQ=255)
                                                                               ccGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCAATGCGCCGCTGAAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACACATCCAAACCCACGCCGCCACAGGCCCGCATCGGtt                                                                   >  1:165430/1‑149 (MQ=255)
                                                                                        cgGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGccccc                                                          <  4:3064/149‑1 (MQ=255)
                                                                                             aaGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGaa                                                      >  8:112868/1‑148 (MQ=255)
                                                                                               gCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAAc                                                   <  7:168079/149‑1 (MQ=255)
                                                                                                  cACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGc                                                <  7:112868/149‑1 (MQ=255)
                                                                                                     tCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTa                                             >  2:11850/1‑149 (MQ=255)
                                                                                                           aGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGCCCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATccccc                                       >  3:36605/1‑149 (MQ=255)
                                                                                                                 cGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTa                                 <  5:175091/149‑1 (MQ=255)
                                                                                                                   aGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTaaa                               <  4:298878/149‑1 (MQ=255)
                                                                                                                   aGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTaaa                               <  1:374542/149‑1 (MQ=255)
                                                                                                                          ttCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCg                                                                      <  2:237234/103‑1 (MQ=255)
                                                                                                                          ttCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCg                                                                      >  1:237234/1‑103 (MQ=255)
                                                                                                                          ttCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGc                                                >  5:12114/1‑125 (MQ=255)
                                                                                                                          ttCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGc                                                <  6:12114/125‑1 (MQ=255)
                                                                                                                          ttCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCt                        >  4:77144/1‑149 (MQ=255)
                                                                                                                          ttCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTACCAAACCCACGCCGCCACAGGc                                                                              >  1:272571/1‑95 (MQ=255)
                                                                                                                          ttCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTACCAAACACACGCCGCCACAGGc                                                                              <  2:272571/95‑1 (MQ=255)
                                                                                                                                  aCCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTc                                                                 >  1:89960/1‑100 (MQ=255)
                                                                                                                                  aCCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTc                                                                 <  2:89960/100‑1 (MQ=255)
                                                                                                                                   ccTAATCATTGCGCCGCGGCAATGACCCCACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTCcgtcgt                <  6:293892/148‑1 (MQ=255)
                                                                                                                                    cTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTTTTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGAGTTAAAGCCCGCTTCCGTCGatt              >  6:149307/1‑149 (MQ=255)
                                                                                                                                    cTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCTCGCGTTAAAGCCCGCTTCCGTCGtt               >  3:3902/1‑148 (MQ=255)
                                                                                                                                    cTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCTCGCGTTAAAGCCCGCTTCCGTCGtt               <  4:3902/148‑1 (MQ=255)
                                                                                                                                    cTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTCCGTCGttt              >  4:128070/1‑149 (MQ=255)
                                                                                                                                     tAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCTCAGACCATGCCAAACAAACGCCGCAACAGGCCCGAATCGGTTTCAGCCCCAAGAAAACTGCGTATCCACCGAGTTAAAGACAGCTTCAGTCGttt              >  5:69212/1‑148 (MQ=255)
                                                                                                                                      ggtCATTGCAGCGGCGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAATCCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTCcgtcgt                <  5:31265/143‑1 (MQ=255)
                                                                                                                                      ggtCATTGCAGCGGCGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTCcgtcgt                >  6:31265/3‑145 (MQ=255)
                                                                                                                                      aaTCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTcc                     <  1:171890/140‑1 (MQ=255)
                                                                                                                                      aaTCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTCCGTCGTTTGa            >  4:93471/1‑149 (MQ=255)
                                                                                                                                      aaTCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGACCGCTTcc                     >  2:171890/1‑140 (MQ=255)
                                                                                                                                        tCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATccc                                         >  5:268980/1‑118 (MQ=255)
                                                                                                                                        tCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATccc                                         <  6:268980/118‑1 (MQ=255)
                                                                                                                                          aTTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTCCGTCGTTTGACata        >  7:261131/1‑149 (MQ=255)
                                                                                                                                          aTTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTCCGTCGTTTGACata        >  8:194536/1‑149 (MQ=255)
                                                                                                                                                ccgcTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGATTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTCCGTCGTTTGACATATCGAAc  >  2:90157/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TCCTGCACTGATTGAAAGCCTGTCACTGTGGGGGGTTTACCCGCGAAGAAGACACCGCCGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTCCGTCGTTTGACATATCGAAC  >  NC_002947/4087341‑4087633

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