Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,415,443 C→T 6.0% A112T (GCC→ACC)  gabP‑V ← gamma‑aminobutyrate permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,415,4430CT 6.0% 199.8 / 2.1 66A112T (GCC→ACC) gabP‑Vgamma‑aminobutyrate permease
Reads supporting (aligned to +/- strand):  ref base C (32/30);  new base T (2/2);  total (34/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 = 7.19e-01

ACCTTGAGGATGGCAAACCAGAACTCGAACTCGCCGTACTTGGCCACGCTGAACAGGTTGGTACCGGCCAGCAGCAACACCGAAGCCAGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGCCGATGGTATAGCCGGCCCAGCGGCCAATGGCCTGGTCGGCATAGGTGGAAAACGATCCGGTGTCGGGGTGCGCCACGGCCATTTCGCCGAGCATGCG  >  NC_002947/5415298‑5415590
                                                                                                                                                 |                                                                                                                                                   
aCCTTGAGGATGGCAAACCAGAACTCGAACTCGCCGTACTTGGCCACGCTGAACAGGTTGGTACCGGCCAGCAGCAACACCGAAGCCAGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCa                                                                                                                                                   <  4:161993/148‑1 (MQ=255)
      aGGATGGCAAACCAGAACTCGAACTCGCCGTACTTGGCCACGCTGAACAGGTTGGTACCGGCCAGCAGCAACACCGAAGCCAGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCtt                                                                                                                                             <  5:121784/148‑1 (MQ=255)
        gATGGCAAACCAGAACTCGAACTCGCCGTACTTGGCCACGCTGAACAGGTTGGTACCGGCCAGCAGCAACACCGAAGCCAGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCg                                                                                                                                           >  1:181416/1‑148 (MQ=255)
        gATGGCAAACCAGAACTCGAACTCGCCGTACTTGGCCACGCTGAACAGGTTGGTACCGGCCAGCAGCAACACCGAAGCCAGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTACAGCACATGCCCGGCGGCCAATGCTTCg                                                                                                                                           >  7:126970/1‑148 (MQ=255)
              aaaCCAGAACTCGAACTCGCCGTACTTGGCCACGCTGAACAGGTTGGTACCGGCCAGCAGCAACACCGAAGCCAGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGaa                                                                                                                                    >  1:114096/1‑149 (MQ=255)
              aaaCCAGAACTCGAACTCGCCGTACTTGGCCACGCTGAACAGGTTGGTACCGGCCAGCAGCAACACCGAAGCCAGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGaa                                                                                                                                    >  4:162152/1‑149 (MQ=255)
                aCCAGAACTCGAACTCGCCGTACTTGGACACGCTGAACAGGTTGCTACCGGCCAGCAGCAACACCGAAGCCAGGGCGCAGATCAAGCAGTCCACCTGCGGGAGCCAGGCGAACCGCACATGCCCggcggc                                                                                                                                                     >  3:280693/1‑130 (MQ=255)
                  cAGAACTCGAACTCGCCGTACTTGGCCACGCTGAACAGGTTGGTACCGGCCAGCAGCAACACCGAAGCCAGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGAc                                                                                                                                <  3:99300/149‑1 (MQ=255)
                    gaactcgaactcGCCGTACTTGGCCACGCTGAACAGGTTGGTACCGGCCAGCAGCAACACCGAAGCCAGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCa                                                                                                                              <  3:162152/149‑1 (MQ=255)
                          gaactcgCCGTACTTGGCCACGCTGAACAGGTTGGTACCGGCCAGCAGCAACACCGAAGCCAGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTAcc                                                                                                                         >  1:242532/1‑148 (MQ=255)
                           aactcgCCGTACTTGGCCACGCTGAACAGGTTGGTACCGGCCAGCAGCAACACCGAAGCCAGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCa                                                                                                                       >  8:329866/1‑149 (MQ=255)
                               cgCCGTACTTGGCCACGCTGACCAGGTTGGTACCGGCCAGCAGCAACACCGAAGCCAGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAac                                                                                                                   >  5:40615/1‑149 (MQ=255)
                               cgCCGTACTTGGCCACGCTGAACAGGTTGGTACCGGCCAGCAGCAACACCGAAGCCAGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAac                                                                                                                   >  4:73397/1‑149 (MQ=255)
                                   gTACTTGGCCACGCTGAACAGGTTGGTACCGGCCAGCAGCAACACCGAAGCCAGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAaccacc                                                                                                               <  6:140356/149‑1 (MQ=255)
                                       ttGGCCACGCTGAACAGGTTGGTACCGGCCAGCAGCAACACCGAAGCCAGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTa                                                                                                           >  4:96823/1‑149 (MQ=255)
                                         ggCCACGCTGAACAGGTTGGTACCGGCCAGCAGCAACACCGAAGCCAGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTAc                                                                                                          <  2:181416/148‑1 (MQ=255)
                                               gCTGAACAGGTTGGTACCGGCCAGCAGCAACACCGAAGCCAGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGc                                                                                                   <  6:20887/149‑1 (MQ=255)
                                                 tGAACAGGTTGGTACCGGCCAGCAGCAACACCGAAGCCAGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACaa                                                                                                        >  1:15369/1‑142 (MQ=255)
                                                 tGAACAGGTTGGTACCGGCCAGCAGCAACACCGAAGCCAGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACaa                                                                                                        <  2:15369/142‑1 (MQ=255)
                                                        gTTGGTACCGGCCAGCAGCAACACCGAAGCCAGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGAACAGTACCCAGAACCACCAGTACAACCAGCCGATGGta                                                                                           >  8:51636/1‑148 (MQ=255)
                                                            gttccGGCCAGCAGCTACACCGAAGCCAGTGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCGAGAACCACCAGTACAACCAGCCGATGGTATAGc                                                                                       <  8:126970/145‑1 (MQ=255)
                                                                 ggCCAGCAGCAACACCGAAGCCAGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGCCGATGGTATAGCCGGccc                                                                                 <  2:114096/149‑1 (MQ=255)
                                                                  gCCAGCAGCAACACCGAAGCCAGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGTCAATGCTTCGATAGGAATGACCAGTACCCa                                                                                                                       <  2:254451/110‑1 (MQ=255)
                                                                  gCCAGCAGCAACACCGAAGCCAGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGTCAATGCTTCGATAGGAATGACCAGTACCCa                                                                                                                       >  1:249424/1‑110 (MQ=255)
                                                                  gCCAGCAGCAACACCGAAGCCAGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGTCAATGCTTCGATAGGAATGACCAGTACCCa                                                                                                                       >  1:254451/1‑110 (MQ=255)
                                                                  gCCAGCAGCAACACCGAAGCCAGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGTCAATGCTTCGATAGGAATGACCAGTACCCa                                                                                                                       <  2:249424/110‑1 (MQ=255)
                                                                  gCCAGCAGCAACACCGAAGCCAGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGCCGATGGTATAg                                                                                        <  4:131074/141‑1 (MQ=255)
                                                                  gCCAGCAGCAACACCGAAGCCAGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGCCGATGGTATAg                                                                                        >  3:131074/1‑141 (MQ=255)
                                                                      gcagcaACACCGAAGCCAGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGCCGATGGTATAGCCGGCCCAGCgg                                                                            <  3:10646/149‑1 (MQ=255)
                                                                       cagcaACACCGAAGCCAGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGCCGATGGTATAGCCGGCCCAGCGGc                                                                           <  2:242532/149‑1 (MQ=255)
                                                                       cagcaACACCGAAGCCAGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGCCGATGGTATAGCCGGCCCAGCGGc                                                                           >  6:52609/1‑149 (MQ=255)
                                                                              acCGAAGCCAGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGCCGATg                                                                                              >  6:257625/1‑123 (MQ=255)
                                                                              acCGAAGCCAGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGCCGATg                                                                                              <  5:257625/123‑1 (MQ=255)
                                                                              acCGAAGCCAGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGCCGATGGTATAGCCGGCCCAGCGGCCAATGGc                                                                    <  6:40615/149‑1 (MQ=255)
                                                                                cGAAGCCAGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAg                                                                                                             >  4:249299/1‑106 (MQ=255)
                                                                                cGAAGCCAGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAg                                                                                                             <  3:249299/106‑1 (MQ=255)
                                                                                 gAAGCCAGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGCCGATGGTATAGCCGGCCCAGCGGCCAATGGCCTg                                                                 <  3:96823/149‑1 (MQ=255)
                                                                                  aaGCCAGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGCCGATGGTATAGCCGGCCCAGCg                                                                             >  5:1238/1‑136 (MQ=255)
                                                                                  aaGCCAGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGCCGATGGTATAGCCGGCCCAGCg                                                                             <  6:1238/136‑1 (MQ=255)
                                                                                     ccAGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGCCGATGGTATAGCCGGCCCAGCGGCCAATGGCCTGGTCg                                                             >  2:263476/1‑149 (MQ=255)
                                                                                       aGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGCCGATGGTATAGCCGGCCCAGCGGCCAATGGCCTGGTCGGc                                                           <  5:52609/149‑1 (MQ=255)
                                                                                         tgCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGCCGATGGTATAGCCGGCCCAGCGGCCAATGGCCTGGTCGGCAt                                                         <  5:95066/148‑1 (MQ=255)
                                                                                            gAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGCCGATGGTATAGCCGGCCCAGCGGCCAATGGCCTGGTCgg                                                            <  6:28151/143‑1 (MQ=255)
                                                                                            gAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGCCGATGGTATAGCCGGCCCAGCGGCCAATGGCCTGGTCgg                                                            >  5:28151/1‑143 (MQ=255)
                                                                                                    aGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGCCGATGGTATAGCCGGCCCAGCGGCCAATGGCCTGGTCGGCATAGGTGGAAAAc                                              >  5:50350/1‑149 (MQ=255)
                                                                                                            aCCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGCCGATGGTATAGCCGGCCCAGCGGCCAATGGCCTGGTCGGCATAGGTGGAAAACGATCCGgt                                      >  7:115026/1‑149 (MQ=255)
                                                                                                            aCCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGCCGATGGTATAGCCGGCCCAGCGGCCAATGGCCTGGTCGGCATAGGTGGAAAACGATCCGgt                                      >  4:291391/1‑149 (MQ=255)
                                                                                                             ccTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGCCGATGGTATAGCCGGCCCAGCGGCCAATGGCCTGGTCGGCATAGGTGGAAAACGATCCGgtg                                     >  2:359028/1‑149 (MQ=255)
                                                                                                                  gggAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGCCGATTGTATAGCCGGCCCAGCGGCCAATGGCCTGGTCGGCATAGGTGGAAAACGATCCGGTGTCggg                                <  7:329866/149‑1 (MQ=255)
                                                                                                                  gggAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGCCGATGGTATAGCCGGCCCAGCGGCCAATGGCCTg                                                                 >  6:123535/1‑116 (MQ=255)
                                                                                                                  gggAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGCCGATGGTATAGCCGGCCCAGCGGCCAATGGCCTg                                                                 <  5:123535/116‑1 (MQ=255)
                                                                                                                        cAGGCGTTCAGCACATGCCCTGCGGCCATTGCTTCGATAGGAAt                                                                                                                                   <  8:108197/44‑1 (MQ=255)
                                                                                                                        cAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGCCGATGGTATAGCCGGCCCAGCGGCCaa                                                                        <  1:101691/103‑1 (MQ=255)
                                                                                                                        cAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGCCGATGGTATAGCCGGCCCAGCGGCCaa                                                                        >  2:101691/1‑103 (MQ=255)
                                                                                                                        cAGGCGTTCAGCACATGCCCGGCGGCCAATGCATCCATAGGAAt                                                                                                                                   >  7:108197/1‑44 (MQ=255)
                                                                                                                        cAGGCGTTCAGCACATGCACGGCAGCCAATGCTTCGATAGGAATGACCAGTACCAAGAACCACAAGTACAACCAGCCGATGGTATAGCCGGCCCAGCCGCCAATGCCCTGGTCGGCAAAGGAGGAAAACTATCCCGTCTCGGtgtgac                           >  8:249587/1‑146 (MQ=255)
                                                                                                                               tCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGCCGATGGTATAGCCGGCCCAGCGGCCAATGGCCTGGTCGGCATAGGTGGAAAACGATCCGGTGTCGGGGTGCGCCACGGc                    <  8:115026/148‑1 (MQ=255)
                                                                                                                                cAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGCCGATGGTATAGCCGGCCCAGCGGCCAATGGCCTGGTCGGCATAGGTGGAAAACGATCCGGTGTCGGGGTGCGCCACGGCCa                  >  8:32820/1‑149 (MQ=255)
                                                                                                                                 aGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGCCGATGGTATAGCCGGCCCAGCGGCCAATGGCCTGGTCGGCATAGGTGGAAAACGATCCGGTGTCGGGGTGCGCCACGGCCAt                 <  3:73397/149‑1 (MQ=255)
                                                                                                                                  gCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGCCGATGGTATAGCCGGCCCAGCGGCCAATGGCCTGGTCGGCATAGGTGGAAAACGATCCg                                        <  8:334918/125‑1 (MQ=255)
                                                                                                                                  gCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGCCGATGGTATAGCCGGCCCAGCGGCCAATGGCCTGGTCGGCATAGGTGGAAAACGATCCg                                        >  7:334918/1‑125 (MQ=255)
                                                                                                                                    acaTGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGCCGATGGTATAGCCGGCCCAGCGGCCAATGGCCTGGTCGGCATAGGTGGAAAACGATCCGGTGTCGGGGTGCGCCACGGCCATTTc              >  6:310386/1‑149 (MQ=255)
                                                                                                                                      aTGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGCCGATGGTATAGCCGGCCCAGCGGCCAATGGCCTGGTCGGCATAGGTGGAAAACGATCCGgtg                                     >  6:339665/1‑124 (MQ=255)
                                                                                                                                      aTGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGCCGATGGTATAGCCGGCCCAGCGGCCAATGGCCTGGTCGGCATAGGTGGAAAACGATCCGgtg                                     <  5:339665/124‑1 (MQ=255)
                                                                                                                                      aTGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGCCGATGGTATAGCCGGCCCAGCGGCCAATGGCCTGGTCGGCATAGGTGGAAAACGATCCGGTGTCGGGGTGCGCCACGGCCATTTCGc            <  3:291391/149‑1 (MQ=255)
                                                                                                                                        gCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGCCGATGGTATAGCCGGCCCAGCGGCCAATGGCCTGGTCGGCATAGGTGGAAAACGATCCGGTGTCGGGGTGCGCCACGGCCATTTCGcc           >  7:29272/1‑148 (MQ=255)
                                                                                                                                          ccGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGCACCACCAGTACAACCAGCCGATGGTATAGCCGGCCCAGCGGCCAATGGCCTGGTCGGCATAGGTGGCAAACGATCCGGGCTCGGGGTGCGCCACGGCCAGTTCGCCGAg        >  3:51905/1‑149 (MQ=39)
                                                                                                                                             gcggcCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGCCGATGGTATAGCCGGCCCAGCGGCCAATGGCCTGGTCGGCATAGGTGGAAAACGATCCGGTGTCGGGGTGCGCCACGGCCATTTCGCCGAGCAt     <  7:249587/149‑1 (MQ=255)
                                                                                                                                                gcCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGCCGATGGTATAGCCGGCCCAGCGGCCAATGGCCTGGTCGGCATAGGTGGAAAACGATCCGGTGTCGGGGTGCGCCACGGCCATTTCGCCGAGCATgcg  <  1:263476/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
ACCTTGAGGATGGCAAACCAGAACTCGAACTCGCCGTACTTGGCCACGCTGAACAGGTTGGTACCGGCCAGCAGCAACACCGAAGCCAGGGCGAAGATCCAGCTGTCCACCTGCGGGAACCAGGCGTTCAGCACATGCCCGGCGGCCAATGCTTCGATAGGAATGACCAGTACCCAGAACCACCAGTACAACCAGCCGATGGTATAGCCGGCCCAGCGGCCAATGGCCTGGTCGGCATAGGTGGAAAACGATCCGGTGTCGGGGTGCGCCACGGCCATTTCGCCGAGCATGCG  >  NC_002947/5415298‑5415590

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