Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,341,586 G→C 5.6% V231V (GTG→GTC panB → 3‑methyl‑2‑oxobutanoate hydroxymethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,341,5860GC 5.6% 229.4 / 4.1 71V231V (GTG→GTCpanB3‑methyl‑2‑oxobutanoate hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base G (32/35);  new base C (2/2);  total (34/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 = 8.45e-01

CGTACCCAGCGAACTGGCGGCGGAAATCACCAATGCCGTTGGCATCCCGGTAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAGCGAACACGGGTTCAGTGCATGAATACAGTCAAGACCGTCCGCGAACTGCGCG  >  NC_002947/5341439‑5341724
                                                                                                                                                   |                                                                                                                                            
cGTACCCAGCGAACTGGCGGCGGAAATCACCAATGCCGTTGGCATCCCGGTAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTg                                                                                                                                              >  2:100818/1‑148 (MQ=255)
cGTACCCAGCGAACTGGCGGCGGAAATCACCAATGCCGTTGGCATCCCGGTAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTg                                                                                                                                              <  3:26414/148‑1 (MQ=255)
     ccAGCGAACTGGCGGCGGAAATCACCAATGCCGTTGGCATCCCGGTAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAAc                                                                                                                                       >  2:344238/1‑149 (MQ=255)
        gCGAACTGGCGGCGGAAATCACCAATGCCGTTGGCATCCCGGTAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTc                                                                                                                                    <  1:344238/149‑1 (MQ=255)
         cGAACTGGCGGCGGAAATCACCAATGCCGTTGGCATCCCGGTAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCTTCACGAACTTcg                                                                                                                                   <  4:144391/149‑2 (MQ=255)
                gcggcgGAAATCACCAATGCCGTTGGCATCCCGGTAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAg                                                                                                                              <  3:136417/147‑1 (MQ=255)
                  ggcggAAATCACCAATGCCGTTGGCATCCCGGTAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGc                                                                                                                           >  1:124919/1‑148 (MQ=255)
                     ggAAATCACCAATGCCGTTGGCATCCCGGTAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAAc                                                                                                                       >  7:347413/1‑149 (MQ=255)
                     ggAAATCACCAATGCCGTTGGCATCCCGGTAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGAGCCAAc                                                                                                                       <  8:335760/149‑1 (MQ=255)
                         aTCACCAATGCCGTTGGCATCCCGGTAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGa                                                                                                                   <  4:193028/149‑1 (MQ=255)
                           cACCAATGCCGTTGGCATCCCGGTAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGata                                                                                                                >  2:39181/1‑149 (MQ=255)
                              cAATGCCGTTGGCATCCCGGTAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATcc                                                                                                             <  3:31452/149‑1 (MQ=255)
                                 tGCCGTTGGCATCCCGGTAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATa                                                                                                          <  8:347413/149‑1 (MQ=255)
                                    cGTTGGCATCCCGGTAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATTCAGGGCCAACCGGATATCCATAgcg                                                                                                       >  7:55953/1‑149 (MQ=255)
                                     gTTGGCATCCCGGTAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGc                                                                                                                                >  3:311782/1‑124 (MQ=255)
                                     gTTGGCATCCCGGTAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGc                                                                                                                                <  4:311782/124‑1 (MQ=255)
                                         gCATCCCGGTAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGcc                                                                                                     >  4:101071/1‑146 (MQ=255)
                                         gCATCCCGGTAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGcc                                                                                                     <  3:101071/146‑1 (MQ=255)
                                           aTCCCGGTAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCtcgt                                                                                                >  2:334601/1‑149 (MQ=255)
                                           aTCCCGGTAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGTCGCTGAGCGGCCGCGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCtcgt                                                                                                >  5:209750/1‑149 (MQ=255)
                                              ccGGTAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGTCGCTGAGCGGCCGTGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCtcgtcg                                                                                              <  1:100818/148‑1 (MQ=255)
                                              ccGGTAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCtcgtcg                                                                                              <  3:23502/148‑1 (MQ=255)
                                                      catcggtGGCGCTGCCGGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTc                                                                                      >  6:252983/7‑148 (MQ=255)
                                                       ggCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCTGATATCCATAGCGCCCTCGTCGCTTACGTCg                                                                                     <  1:39181/148‑1 (MQ=255)
                                                        gCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAg                                                                                   >  6:270213/1‑149 (MQ=255)
                                                        gCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAg                                                                                   >  5:187528/1‑149 (MQ=255)
                                                        gCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAg                                                                                   >  8:58296/1‑149 (MQ=255)
                                                           ttGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCg                                                                                <  1:334601/149‑1 (MQ=255)
                                                               cgcCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGTca                                                                            <  5:275791/149‑1 (MQ=255)
                                                                  cGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGTGGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGTcaagc                                                                         <  1:137201/149‑1 (MQ=255)
                                                                      aGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGTCAAGCAAGt                                                                     <  7:58296/149‑1 (MQ=255)
                                                                        cgcCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGTCAAGCAAGTCa                                                                   <  4:279456/149‑1 (MQ=255)
                                                                           cACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGt                                                                              <  3:150896/135‑1 (MQ=255)
                                                                           cACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGt                                                                              >  4:150896/1‑135 (MQ=255)
                                                                             ccGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTc                                                              >  5:287057/1‑149 (MQ=255)
                                                                                 tGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAg                                                          >  5:373261/1‑149 (MQ=255)
                                                                                     cAGGTACTGGTACTGCACGACATTCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAg                                                                                   >  4:228380/1‑120 (MQ=255)
                                                                                     cAGGTACTGGTACTGCACGACATTCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAg                                                                                   <  3:228380/120‑1 (MQ=255)
                                                                                     cAGGTACTGGTACTGCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCATGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAg                                                      >  6:266927/1‑149 (MQ=255)
                                                                                           ctggtactgCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAGCGAac                                                 <  4:341906/148‑1 (MQ=255)
                                                                                            tggtactgCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAGCGAacac                                               >  6:73927/1‑149 (MQ=255)
                                                                                              gtactgCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAGCGAACACgg                                             >  7:219435/1‑149 (MQ=255)
                                                                                              gtactgCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAGCGAACACgg                                             >  6:100974/1‑149 (MQ=255)
                                                                                               tactgCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAGCGAACACgg                                             <  7:45028/148‑1 (MQ=255)
                                                                                                 ctgCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAATTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAGCGAACACGGGtt                                          <  8:328284/149‑1 (MQ=255)
                                                                                                 ctgCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAGCGAACACGGGt                                           <  6:373261/148‑1 (MQ=255)
                                                                                                  tgCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAAGCTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAGCGAACACGGGtt                                          >  1:101493/1‑149 (MQ=255)
                                                                                                    cACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAGCGAACACGGGTTCAg                                       >  3:84078/1‑149 (MQ=255)
                                                                                                       gACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAGTGAACACgg                                             <  2:209781/140‑1 (MQ=255)
                                                                                                       gACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAGTGAACACgg                                             >  1:209781/1‑140 (MQ=255)
                                                                                                       gACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAGCGAACACGGGTTCAGTg                                     >  3:269124/1‑148 (MQ=255)
                                                                                                         cATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAGCGAACACGGTTTCAGTGCAt                                  <  5:73927/149‑1 (MQ=255)
                                                                                                         cATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAGCGAACACGGGTTCAGTGCAt                                  <  6:287057/149‑1 (MQ=255)
                                                                                                             cTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAGCGAACACGGGTTCAGTGCATGAAt                              <  2:124919/149‑1 (MQ=255)
                                                                                                                ggCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATGATCCAAACCGCCGACGTCGCTTACGACCAGGCGGGCAAGCAAGTCAGTTTCCCAGGAACCGAACACGGGTTAAGTGCATGATTa                             >  1:111027/1‑148 (MQ=255)
                                                                                                                ggCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAGCGAACACGGGTTCAGTGCATGAAt                              <  5:84366/146‑1 (MQ=255)
                                                                                                                ggCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAGCGAACACGGGTTCAGTGCATGAAt                              >  6:84366/1‑146 (MQ=255)
                                                                                                                ggCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAGCGAACACGGGTTCAGTGCATGAATa                             <  5:252983/147‑1 (MQ=255)
                                                                                                                   cTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAGCGAACACGGGTTCAGTGCATGAATACAGTc                        <  5:270213/149‑1 (MQ=255)
                                                                                                                         cTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAGCGAACACgg                                             <  8:191608/122‑1 (MQ=255)
                                                                                                                         cTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAGCGAACACgg                                             >  7:191608/1‑122 (MQ=255)
                                                                                                                         cTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAGCGAACACGGGTTCAGTGCATGAATACAGTCAAGAcc                  <  2:101493/149‑1 (MQ=255)
                                                                                                                             gCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAGCGAACACGGGTTCAGTGCATGAATACAGTCAAGACCGTcc              >  1:242173/1‑149 (MQ=255)
                                                                                                                                   gggTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAGCGAACACGGGTTCAGTGCATGAATACAGTCAAGACCGTCCGCGacc        >  3:58487/1‑147 (MQ=255)
                                                                                                                                    ggTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAGCGAACACGGGTTCAGTGCATGAATACAGTCAAGACCGTCCGCGAACt       <  6:187528/149‑1 (MQ=255)
                                                                                                                                      tACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAGCGAACACGGGTTCAGTGCATGAATACAGTCAAGACCGTCCGCGAACTgc     <  2:242173/149‑1 (MQ=255)
                                                                                                                                        atgAAGTTCTTCACGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGCCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAGCGAACACGGGTTCAGTGCATGaa                               <  6:266934/119‑1 (MQ=255)
                                                                                                                                        atgAAGTTCTTCACGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGACGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAGCGAACACGGGTTCAGTGCATGaa                               >  5:266934/3‑121 (MQ=255)
                                                                                                                                         tgAAGTTCTTCACGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAGCGAACACGGGTTCAGTGCATGAATGCAGTCAAGACCGTCCGCGAACTgcgcg  >  4:133847/2‑149 (MQ=255)
                                                                                                                                         cGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAGCGAACACGGGTTCAGt                                      >  6:66223/1‑113 (MQ=255)
                                                                                                                                         cGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAGCGAACACGGGTTCAGt                                      <  5:66223/113‑1 (MQ=255)
                                                                                                                                         cGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAGCGAACACGGGTTCAGTGCATGTATACAGTCAAGACCGTCCGCGAACt       <  4:221368/144‑1 (MQ=255)
                                                                                                                                         cGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAGCGAACACGGGTTCAGTGCATGTATACAGTCAAGACCGTCCGCGAACt       >  3:221368/1‑144 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
CGTACCCAGCGAACTGGCGGCGGAAATCACCAATGCCGTTGGCATCCCGGTAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGTCGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAGCGAACACGGGTTCAGTGCATGAATACAGTCAAGACCGTCCGCGAACTGCGCG  >  NC_002947/5341439‑5341724

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: