Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,776,699 T→C 7.2% I339T (ATC→ACC)  betA‑II → choline dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,776,6990TC 7.2% 196.9 / 4.7 69I339T (ATC→ACC) betA‑IIcholine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (31/33);  new base C (3/2);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

ACCTTCCGGGCGTTGGTGAAAACCTTCAGGACCACCTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG  >  NC_002947/5776552‑5776847
                                                                                                                                                     |                                                                                                                                                     
aCCTTCCGGGCGTTGGTGAAAACCTTCAGGACCACCTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATc                                                                                                                                                      <  4:87318/149‑1 (MQ=255)
       gggCGTTGGTGAAAACCTTCAGGACCACCTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCGCCGATACCGGTGCCa                                                                                                                                               <  1:401681/149‑1 (MQ=255)
       gggCGTTGGTGAAAACCTTCAGGACCACCTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCa                                                                                                                                               <  8:32034/149‑1 (MQ=255)
       gggCGTTGGTGAAAACCTTCAGGACCACCTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCa                                                                                                                                               >  2:197130/1‑149 (MQ=255)
       gggCGTTGGTGAAAACCTTCAGGACCACCTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCa                                                                                                                                               <  2:311941/149‑1 (MQ=255)
           gTTGGTGAAAACCTTCAGGACCACCTGGAACTGTCCCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTTGATGGTCAATGGCACCGGTATCGGCGCCAGCaa                                                                                                                                           <  1:259929/149‑1 (MQ=255)
                 gAAAACCTTCAGGACCACCTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGtt                                                                                                                                     >  6:325313/1‑149 (MQ=255)
                 gAAAACCTTCAGGACCACCTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGtt                                                                                                                                     >  7:281132/1‑149 (MQ=255)
                  aaaaCCTTCAGGACCACCTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTc                                                                                                                                    <  8:231230/149‑1 (MQ=255)
                      ccTTCAGGACCACCTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAgg                                                                                                                                >  5:363960/1‑149 (MQ=255)
                           aGGAGCACGTGGAAGTGTACCTTCAGTTCGCCTGCACCCTGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCggg                                                                                                                           <  6:270615/149‑1 (MQ=255)
                           aGGACCACCTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCggg                                                                                                                           <  7:240182/149‑1 (MQ=255)
                             gACCACCTGGAACTGTACCTGCAGTACGCCTGCACCCAGAAGGTTTCGCTATACCCTTCGCTGCTGTGGTGGAACCAGCCGGC‑‑TATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTAGAGGCAggggg                                                                                                                         >  5:385096/1‑149 (MQ=255)
                                caccTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGtt                                                                                                                      >  4:144/1‑149 (MQ=255)
                                caccTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGtt                                                                                                                      <  6:363960/149‑1 (MQ=255)
                                  ccTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCa                                                                                                                    >  2:130853/1‑149 (MQ=255)
                                   cTGGAACTGTACCTGCAGTACGCCTGCACCCAGCTGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCAt                                                                                                                   >  8:90470/1‑149 (MQ=255)
                                   cTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCAt                                                                                                                   >  1:234150/1‑149 (MQ=255)
                                     ggAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATcc                                                                                                                 >  4:336373/1‑149 (MQ=255)
                                        aCTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGa                                                                                                              <  8:119289/149‑1 (MQ=255)
                                         cTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGAc                                                                                                             <  2:234150/149‑1 (MQ=255)
                                         cTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGAc                                                                                                             <  8:80144/149‑1 (MQ=255)
                                              ccTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCggg                                                                                                                           <  2:90254/130‑1 (MQ=255)
                                              ccTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCggg                                                                                                                           >  1:90254/1‑130 (MQ=255)
                                               cTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCAGGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGt                                                                                                         >  1:120398/1‑149 (MQ=255)
                                               cTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCggggg                                                                                                                         >  8:122886/1‑131 (MQ=255)
                                               cTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCggggg                                                                                                                         <  7:122886/131‑1 (MQ=255)
                                                  cAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCg                                                                                                      <  5:416458/147‑1 (MQ=255)
                                                       cGCCTGCACCCAGCTGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCa                                                                                              <  7:90470/149‑1 (MQ=255)
                                                             cACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCTGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGc                                                                                        <  5:222700/149‑1 (MQ=255)
                                                                ccAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATcc                                                                                                                 <  4:348812/122‑1 (MQ=255)
                                                                ccAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATcc                                                                                                                 >  3:348812/1‑122 (MQ=255)
                                                                ccAGCCGGCTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCg                                                                                      >  8:145233/1‑148 (MQ=255)
                                                                     ctgTTTCGCTATACCCGTCGCTGCTGTTGTGGAACAAGCCGGC‑‑GATAGGTGCAGAGTGGATGTTCAATGTCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATc                                                                                <  6:385096/147‑1 (MQ=255)
                                                                        cttCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAg                                                                             <  7:145233/148‑1 (MQ=255)
                                                                            gCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTAc                                                                          <  1:327068/148‑1 (MQ=255)
                                                                            gCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTAc                                                                          >  3:332315/1‑148 (MQ=255)
                                                                                  cccGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTc                                                                    <  8:319684/148‑1 (MQ=255)
                                                                                     gTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGt                                                                                           >  8:220497/1‑122 (MQ=255)
                                                                                     gTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGt                                                                                           <  7:220497/122‑1 (MQ=255)
                                                                                         ctgctgTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGt                                                            <  1:197130/149‑1 (MQ=255)
                                                                                          tgctgTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGTCCGGGGGGTTCATCCGGACTCGTGCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGtt                                                           >  6:278532/1‑149 (MQ=255)
                                                                                          tgctgTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCgggg                                                                                                                          >  2:191095/1‑87 (MQ=255)
                                                                                           gctgTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTg                                                          >  6:290737/1‑149 (MQ=255)
                                                                                            ctgTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTg                                                          <  6:76206/148‑1 (MQ=255)
                                                                                                  tggAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTc                                                                    >  8:148308/1‑132 (MQ=255)
                                                                                                  tggAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTc                                                                    <  7:148308/132‑1 (MQ=255)
                                                                                                  tggAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTaa                                                   >  6:108396/1‑149 (MQ=255)
                                                                                                   ggAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAAc                                                  >  8:333328/1‑149 (MQ=255)
                                                                                                      aCCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGAt                                                      <  6:170893/142‑1 (MQ=255)
                                                                                                      aCCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGAt                                                      <  6:164829/142‑1 (MQ=255)
                                                                                                      aCCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGAt                                                      >  5:170893/1‑142 (MQ=255)
                                                                                                      aCCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGAt                                                      >  5:164829/1‑142 (MQ=255)
                                                                                                      aCCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTa                                                >  6:20747/1‑148 (MQ=255)
                                                                                                            cGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGc                                         <  4:186281/149‑1 (MQ=255)
                                                                                                             ggC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCt                                        <  3:336373/149‑1 (MQ=255)
                                                                                                                          ccGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGt                              <  3:144/148‑1 (MQ=255)
                                                                                                                                     ttCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAAt                                   <  4:338008/132‑1 (MQ=255)
                                                                                                                                     ttCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAAt                                   >  3:338008/1‑132 (MQ=255)
                                                                                                                                      tCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAAc                                                  >  3:107950/1‑116 (MQ=255)
                                                                                                                                      tCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAAc                                                  <  4:107950/116‑1 (MQ=255)
                                                                                                                                        aaTGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGa                                                       <  1:283009/110‑1 (MQ=255)
                                                                                                                                        aaTGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGa                                                       >  2:283009/1‑109 (MQ=255)
                                                                                                                                          tGGCGCCGATACCGGTGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGc                                >  6:240050/1‑130 (MQ=255)
                                                                                                                                          tGGCGCCGATACCGGTGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTc              >  6:38654/1‑148 (MQ=255)
                                                                                                                                          tGGCGCCGATACCGGTGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTc              >  1:119725/1‑148 (MQ=255)
                                                                                                                                             cgccGATACCGGTGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAgg          <  5:38654/147‑1 (MQ=255)
                                                                                                                                              aCCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGgc         >  7:280933/1‑149 (MQ=255)
                                                                                                                                                cGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGgcg        >  5:214522/1‑148 (MQ=255)
                                                                                                                                                     tCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATg  <  8:281132/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACCTTCCGGGCGTTGGTGAAAACCTTCAGGACCACCTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG  >  NC_002947/5776552‑5776847

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