Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 1,166,671 T→C 6.7% E136E (GAA→GAG)  hcaR_1 ← Hca operon transcriptional activator HcaR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1431,166,6710TC 6.7% 186.6 / 3.2 62E136E (GAA→GAG) hcaR_1Hca operon transcriptional activator HcaR
Reads supporting (aligned to +/- strand):  ref base T (22/34);  new base C (2/2);  total (25/37)
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.88e-01

ACCCGCCAAGGGGTGTCCAGCGGCGACCGCCAGCAACAAGGGTTCGTCGAACAACTCGTGCAGCACCACGCCTTCATGCTGGCGCAGCACCGGCAAGCGCAACAGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAGCGAAACCCCTGGGTTGTGCTTCAGGTAGTTCCCCAGGCGAGCCAGAAGCCCGCCGGTCAGTGGCACGGTGCTGGAATGGTTCAGGCGCAAG  >  ACB143/1166536‑1166800
                                                                                                                                       |                                                                                                                                 
aCCCGCCAAGGGGTGTCCAGCGGCGACCGCCAGCAACAAGGGTTCGTCGAACAACTCGTGCAGCACCACGCCTTCATGCTGGCGCAGCACCGGCAAGCGCAACAGGCCGATATCCAGCCGTCCGGCAGCAATGTCt                                                                                                                                   >  3:279577/1‑136 (MQ=255)
aCCCGCCAAGGGGTGTCCAGCGGCGACCGCCAGCAACAAGGGTTCGTCGAACAACTCGTGCAGCACCACGCCTTCATGCTGGCGCAGCACCGGCAAGCGCAACAGGCCGATATCCAGCCGTCCGGCAGCAATGTCt                                                                                                                                   >  8:265726/1‑136 (MQ=255)
  ccGCCAAGGGGTGTCCAGCGGCGACCGCCAGCAACAAGGGTTCGTCGAACAACTCGTGCAGCACCACGCCTTCATGCTGGCGCAGCACCGGCAAGCGCAACAGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTc                                                                                                                                 >  7:195636/1‑136 (MQ=255)
       aaGGGGTGTCCAGCGGCGACCGCCAGCAACAAGGGTTCGTCGAACAACTCGTGCCGCACCACGCCGTCATGCTGGCGCAGCACCGGCAAGCGCAACAGGCCGATATCCAGCCGTCCGGCAGCAATGTCGTCCAGCt                                                                                                                            >  1:227781/1‑136 (MQ=255)
          gggTGTCCAGCGGCGACCGCCAGCAACAAGGGTTCGTCGAACAACTCGTGCAGCACCACGCCTTCATGCTGGCGCAGCACCGGCAAGCGCAACAGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCtgtg                                                                                                                         >  4:70737/1‑136 (MQ=255)
               tCCAGCGGCGACCGCCAGCAACAAGGGTTCGTCGAACAACTCGTGCAGCACCACGCCTTCATGCTGGCGCAGCACCGGCAAGCGCAACAGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTcc                                                                                                                    <  6:87677/136‑1 (MQ=255)
                  aGCGGCGACCGCCAGCAACAAGGGTTCGTCGAACAACTCGTGCAGCACCACGCCTTCATGCTGGCGCAGCACCGGCAAGCGCAACAGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCgct                                                                                                                 >  7:213746/1‑134 (MQ=255)
                  aGCGGCGACCGCCAGCAACAAGGGTTCGTCGAACAACTCGTGCAGCACCACGCCTTCATGCTGGCGCAGCACCGGCAAGCGCAACAGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGAt                                                                                                                 >  3:295161/1‑136 (MQ=255)
                  aGCGGCGACCGCCAGCAACAAGGGTTCGTCGAACAACTCGTGCAGCACCACGCCTTCATGCTGGCGCAGCACCGGCAAGCGCAACAGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGAt                                                                                                                 >  3:203300/1‑136 (MQ=255)
                  aGCGGCGACCGCCAGCAACAAGGGTTCGTCGAACAACTCGTGCAGCAACACGCCTTCATGCTGGAGCAGCACCGGCAAGCGCAACAGCCCGATATACAGCCGTCCGGCAGCAAAGTCTTCCAGCTGTGCCTCCGAt                                                                                                                 >  7:71765/1‑136 (MQ=255)
                           cGCCAGCAACAAGGGTTCGTCGAACAACTCGTGCAGCACCACGCCTTCATGCTGGCGCAGCACCGGCAAGCGCAACAGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACtgctgc                                                                                                        <  8:237783/136‑1 (MQ=255)
                            gCCAGCAACAAGGGTTCGTCGAACACCTCGTGCAGCGCCACGCCTTCATGCTGGAGCAGCACCGGCAAGCGCAACAGGCCGATATCCAGCCGTCCGGCAGCACAGTCTTCCAGCTGTGCCTCCGATGACTGCTgcg                                                                                                       >  5:104908/1‑136 (MQ=255)
                            gCCAGCAACAAGGGTTCGTCGAACAACTCGTGCAGCACCACGCCTTCATGCTGGCGCAGCACCGGCCAGCGCAACAGGCCGATATCCAGCCGTCCGGAAGCAATGTCTTCCAGCTGTGCCTCCGATGTCTGCTgcg                                                                                                       >  5:231002/1‑136 (MQ=255)
                            gCCAGCAACAAGGGTTCGTCGAACAACTCGTGCAGCACCACGCCTTCATGCTGGCGCAGCACCGGCAAGCGCAACAGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACTGCTgcg                                                                                                       >  3:131758/1‑136 (MQ=255)
                            gCCAGCAACAAGGGTTCGTCGAACAACTCGTGCAGCACCACGCCTTCATGCTGGCGCAGCACCGGCAAGCGCAACAGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACTGCTgcg                                                                                                       >  2:169034/1‑136 (MQ=255)
                                      aGGGTTCTTCGAACAACTCGTGCAGCACCACTCCTTCATGCTGGCGCAGCACCGTCTAGCGCAACAGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAg                                                                                             <  8:201831/136‑1 (MQ=255)
                                      aGGGTTCGTCGAACAACTCGTGCAGCACCACGCCTTCATGCTGGCGCAGCACCGGCAAGCGCAACAGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAg                                                                                             <  8:10424/136‑1 (MQ=255)
                                      aGGGTTCGTCGAACAACTCGTGCAGCACCACGCCTTCATGCTGGCGCAGCACCGGCAAGCGCAACAGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAg                                                                                             <  5:289031/136‑1 (MQ=255)
                                      aGGGTTCGTCGAACAACTCGTGCAGCACCACGCCTTCATGCTGGCGCAGCACCGGCAAGCGCAACAGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAg                                                                                             <  7:269645/136‑1 (MQ=255)
                                      aGGGTTCGTCGAACAACTCGTGCAGCACCACGCCTTCATGCTGGCGCAGCACCGGCAAGCGCAACAGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAg                                                                                             <  1:185844/136‑1 (MQ=255)
                                          ttCGTCGAACAACTCGTGCAGCACCACGCCTTCATGCTGGCGCAGCACCGGCAAGCGCAACAGTCCGATATCCAGCCGTCCGGCAGCAATGTCCTCCAGCTGTGCCTCCGATGACTGCTGCTCAATCTCCAGCGa                                                                                          >  3:15682/1‑135 (MQ=255)
                                          ttCGTCGAACAACTCGTGCAGCACCACGCCTTCATGCTGGCGCAGCACCGGCAAGCGCAACAGGCCGATATCCAGCCGTCCGGCAGCAATGTCCTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAGCGa                                                                                          >  3:15679/1‑135 (MQ=255)
                                            cgtcgAACAACTCGTGCAGCACCACGCCTTCATGCTGGCGCAGCACCGGCAAGCGCAACAGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAGCGaaa                                                                                        <  3:244527/135‑1 (MQ=255)
                                            cgtcgAACAACTCGTGCAGCACCACGCCTTCATGCTGGCGCAGCACCGGCAAGCGCAACAGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAGCGaaa                                                                                        <  1:229382/135‑1 (MQ=255)
                                            cgtcgAACAACTCGTGCAGCACCACGCCTTCATGCTGGCGCAGCACCGGCAAGCGCAACAGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAGCGaaa                                                                                        <  5:102779/135‑1 (MQ=255)
                                               cgAACAACTCGTGCAGCACCACGCCTTCATGCTGGCGCAGCACCGGCAAGCGCAACAGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAGCGAAAcccc                                                                                    >  7:147534/1‑136 (MQ=255)
                                               cgAACAACTCGTGCAGCACCACGCCTTCATGCTGGCGCAGCACCGGCAAGCGCAACAGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAGCGAAAcccc                                                                                    >  6:59312/1‑136 (MQ=255)
                                               cgAACAACTCGTGCAGCACCACGCCTTCATGCTGGCGCAGCACCGGCAAGCGCAACAGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAGCGAAAcccc                                                                                    >  1:208868/1‑136 (MQ=255)
                                                gAACAACTCGTGCAGCACCACGCCTTCATGCTGGCGCAGCACCGGCAAGCGCAACAGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAGCGAAACCCCt                                                                                   <  6:221140/136‑1 (MQ=255)
                                                  acaacTCGTGCAGCACCACGCCTTCATGCTGGCGCAGCACCGGCAAGCGCAACAGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAGCGAAACCCCTgg                                                                                 <  8:101971/136‑1 (MQ=255)
                                                               caccacGCCTTCATGCTGGCGCAGCACCGGCAAGCGCAACAGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACTGATGCGCAATCTCCAGCGAAACCCCTGGGTTGTGCTTCAgg                                                                    <  3:257166/136‑1 (MQ=255)
                                                                  cacGCCTTCATGCTGGCGCAGCACCGGCAAGCGCAACCGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAGCGAAACCCCTGGGTTGTGCTTCAGGTAg                                                                 >  3:48518/1‑136 (MQ=255)
                                                                     gCCTTCATGCTGGCGCAGCACCGGCAAGCGCAACAGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAGCGAAACCCCTGGGTTGTGCTTCAGGTAGTTc                                                              <  8:147534/136‑1 (MQ=255)
                                                                          cATGCTGGCGCAGCACCGGCAAGCGCAACAGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAGCGAAACCCCTGGGTTGTGCTTCAGGTAGTTcccc                                                           <  7:50283/134‑1 (MQ=255)
                                                                                  cgcAGCACCGGCAAGCGCAACAGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAGCGAAACCCCTGGGTTGTGCTTCAGGTAGTTCCCCAGGCGAGCCa                                                 <  5:107212/136‑1 (MQ=255)
                                                                                    cagcaCCGGCAAGCGCAACAGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAg                                                                                             >  4:95495/1‑90 (MQ=255)
                                                                                    cagcaCCGGCAAGCGCAACAGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAg                                                                                             <  3:95495/90‑1 (MQ=255)
                                                                                    cagcaCCGGCAAGCGCAACAGGCCGATATCAAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAGCGAAACCCCTGGGTAGTGCTTCAGATAGTTCCCCAGGCGAGCCAGa                                               >  7:245822/1‑136 (MQ=255)
                                                                                            gCAAGCGCAACAGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAGCGAAACCCCTGGGTTGTGCTTCAGGTATTTCCCCAGGCGAGCCAGAAGCccgcc                                       <  1:104955/136‑1 (MQ=255)
                                                                                            gCAAGCGCAACAGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAGCGAAACCCCTGGGTTGTGCTTCAGGTAGTTCCCCAGGCGAGCCAGAAGCccgcc                                       <  3:148318/136‑1 (MQ=255)
                                                                                                    aaCAGGCCTATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGGCTCCGATGACTGCTGCGCAATCTCCAGCGAAACCCCTGGGGTGTGCTTCAGGTAGTTCCCCAGGCGAGCCAGAAGCCCGCCGGTCAGTg                               <  8:71765/136‑1 (MQ=255)
                                                                                                    aaCAGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAGCGAAACCCCTGGGTTGTGCTTCAGGTAGTTCCCCAGGCGAGCCAGAAGCCCGCCGGTCAGTg                               <  6:99291/136‑1 (MQ=255)
                                                                                                    aaCAGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAGCGAAACCCCTGGGTTGTGCTTCAGGTAGTTCCCCAGGCGAGCCAGAAGCCCGCCGGTCAGTg                               >  8:89947/1‑136 (MQ=255)
                                                                                                    aaCAGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAGCGAAACCCCTGGGTTGTGCTTCAGGTAGTTCCCCAGGCGAGCCAGAAGCCCGCCGGTCAGTg                               >  3:231512/1‑136 (MQ=255)
                                                                                                    aaCAGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAGCGAAACCCCTGGGTTGTGCTTCAGGTAGTTCCCCAGGCGAGCCAGAAGCCCGCCGGTCAGTg                               >  6:126518/1‑136 (MQ=255)
                                                                                                    aaCAGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAGCGAAACCCCTGGGTTGTGCTTCAGGTAGTTCCCCAGGCGAGCCAGAAGCCCGCCGGTCAGTg                               <  5:71491/136‑1 (MQ=255)
                                                                                                    aaCAGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAGCGAAACCCCTGGGTTGTGCTTCAGGTAGTTCCCCAGGCGAGCCAGAAGCCCGCCGGTCAGTg                               <  5:59312/136‑1 (MQ=255)
                                                                                                    aaCAGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAGCGAAACCCCTGGGTTGTGCTTCAGGTAGTTCCCCAGGCGAGCCAGAAGCCCGCCGGTCAGTg                               <  2:243869/136‑1 (MQ=255)
                                                                                                    aaCAGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAGCGAAACCCATGGGTTGTGCTTCAGGTAGTTCCCCAGGCGAGCCAGAAGCCAGCCGCTCAGTg                               >  6:126522/1‑136 (MQ=255)
                                                                                                         gCCGATATCCCGCCGTCCGGCAGCAATGTCTTCCACCTGTGCCTCCGATGACTGCTGCGCAATCTCCAGCGAAACCCATGGGTTGTGCTTCAGGTAGTTCCCCAGGCGAGCCAGAAGCCAGCCGGTCAGTggaagg                          >  6:184225/1‑132 (MQ=255)
                                                                                                              tatCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCa                                                                                              <  6:78433/63‑1 (MQ=255)
                                                                                                              tatCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCa                                                                                              >  5:78433/1‑63 (MQ=255)
                                                                                                               atCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGAGCCTCCGATGACTGCTGCGCAATCTCCAGCGAAACCCCTGGGTTGTGCTTCAGGTAGTTCCCCAGGCGAGCCAGAAGCCCGCCGGTCAGTGGCACGGTGCTg                    <  5:71366/136‑1 (MQ=255)
                                                                                                                tCCAGCCGTCCGTCATCAATGTCCTCCAGCTGTGCCTCCGATGACTTCTGCGCAATCTCCAGCGAAACCCCTGGGTTGTGCTTCAGGTAGTTCCCCAGGCGAGCCAGAAGCCCGCCGGTCAGTGGCACGGTGCTg                    <  4:15679/135‑1 (MQ=255)
                                                                                                                tCCAGCCGTCCGGCAGCAATTTCCTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAGCGAAACCCCTGGGTTGTGCTTCAGGTAGTTCCCCAGGCGAGCCAGAAGCCCGCCGGTCAGTGGCACGGTGCTg                    <  4:15682/135‑1 (MQ=255)
                                                                                                                tCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAGCGAAACCCCTGGGTTGTGCTTCAGGTAGTTCCCCAGGCGAGCCAGAAGCCCGCCGGTCAGTGGCACGGTGCTg                    <  6:267515/135‑1 (MQ=255)
                                                                                                                tCCAGCCGTCCGGCAGCAATGTCATCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAGCGAAACCCCTGGGTTGTGCTTCAGGTAGTTCCCCAGGCGAGCCAGAAGCCCGCCGGTCAGTGGCACGGTGCTg                    <  2:216518/135‑1 (MQ=255)
                                                                                                                    gCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAGCGAAACCCCTGGGTTGTGCTTCAGGTAGTTCCCCAGGCGAGCCAGAAGCCCGCCGGTCAGTGGCACGGTGCTGGAATg               <  5:237502/136‑1 (MQ=255)
                                                                                                                    gCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAGCGAAACCCCTGGGTTGTGCTTCAGGTAGTTCCCCAGGCGAGCCAGAAGCCCGCCGGTCAGTGGCACGGTGCTGGAATg               <  2:46839/136‑1 (MQ=255)
                                                                                                                            gcagcaATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAGCGAAACCCCTGGGTTGTGCTTCAGGTAGGTCCCCAGGCGAGCCAGAAGCCCGCCGGTCAGTGGCACGGTGCTGGAATGGTTCAGgc       <  2:142613/136‑1 (MQ=255)
                                                                                                                              agcaATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAGCGAAACCCCTGGGTTGTGCTTCAGGTAGTTCCCCAGGCGAGCCAGAAGCCCGCCGGTCAGTGGCACGGTGCTGGAATGGTTCAGgcgc     <  5:161002/136‑1 (MQ=255)
                                                                                                                              agcaATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAGCGAAACCCCTGGGTTGTGCTTCAGGTAGTTCCCCAGGCGAGCCAGAAGCCCGCCGGTCAGTGGCACGGTGCTGGAATGGTTCAGgcgc     <  5:182510/136‑1 (MQ=255)
                                                                                                                              agcaATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAGCGAAACCCCTGGGTTGTGCTTCAGGTAGTTCCCCAGGCGAGCCAGAAGCCCGCCGGTCAGTGGCACGGTGCTGGAATGGTTCAGgcgc     <  2:189203/136‑1 (MQ=255)
                                                                                                                                 aaTGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAGCGAAACCCCTGGGTTGTGCTTCAGGTAGTTCCCCAGGCGAGCCAGAAGCCCGCCGGTCAGTGGCACGGTGCTGGAATGGTTCAGGCGCAAg  <  5:18415/136‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                 
ACCCGCCAAGGGGTGTCCAGCGGCGACCGCCAGCAACAAGGGTTCGTCGAACAACTCGTGCAGCACCACGCCTTCATGCTGGCGCAGCACCGGCAAGCGCAACAGGCCGATATCCAGCCGTCCGGCAGCAATGTCTTCCAGCTGTGCCTCCGATGACTGCTGCGCAATCTCCAGCGAAACCCCTGGGTTGTGCTTCAGGTAGTTCCCCAGGCGAGCCAGAAGCCCGCCGGTCAGTGGCACGGTGCTGGAATGGTTCAGGCGCAAG  >  ACB143/1166536‑1166800

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