Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,446,494 G→C 8.1% L357L (CTC→CTG)  nicC ← 6‑hydroxynicotinate 3‑monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,446,4940GC 8.1% 191.8 / 6.5 62L357L (CTC→CTG) nicC6‑hydroxynicotinate 3‑monooxygenase
Reads supporting (aligned to +/- strand):  ref base G (32/25);  new base C (2/3);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

GCCGCCATAGCGCAGGTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGCCGCGTCCTCGATGGCC  >  NC_002947/4446348‑4446638
                                                                                                                                                   |                                                                                                                                                
gccgccATAGCGCAGGTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGG‑ATCTTCCTGGCTGTagag                                                                                                                                                >  4:83934/1‑149 (MQ=255)
  cgccATAGCGCAGGTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGG‑ATCTTCCTGGCTCTACAGc                                                                                                                                               >  1:224863/1‑148 (MQ=255)
    ccATAGCGCAGGTAATGCTGGCGGAGGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTATCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAg                                                                                                                                            <  2:91102/149‑1 (MQ=255)
     cATAGCGCAGGTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGG‑ATCTTCCTGGCTCTACAGCCAgg                                                                                                                                           <  4:53227/149‑1 (MQ=255)
      aTAGCGCAGGTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGgt                                                                                                                                          >  2:47389/1‑149 (MQ=255)
       tAGCGCAGGTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGG‑ATCTTCCTGGCTCTACAGCCAggaa                                                                                                                                         <  2:224863/149‑3 (MQ=255)
           gcAGGTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAACCCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTgg                                                                                                                                     >  8:159839/1‑149 (MQ=255)
           gcAGGTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGCTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGCGGTCGTAGCCGGAAACCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTgg                                                                                                                                     <  7:278283/149‑1 (MQ=255)
               gTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGtt                                                                                                                                 >  1:220720/1‑149 (MQ=255)
               gTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGtt                                                                                                                                 <  3:275241/149‑1 (MQ=255)
                tAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTg                                                                                                                                <  2:42497/149‑1 (MQ=255)
                         cGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGAt                                                                                                                       <  1:305215/149‑1 (MQ=255)
                                 aTTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTg                                                                                                               >  1:187854/1‑149 (MQ=255)
                                 aTTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTg                                                                                                               >  3:185261/1‑149 (MQ=255)
                                              ttACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAATCCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGctct                                                                                                  <  3:201570/149‑1 (MQ=255)
                                                   accaGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCg                                                                                             >  7:114373/1‑149 (MQ=255)
                                                   accaGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCg                                                                                             >  1:219639/1‑149 (MQ=255)
                                                     caGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGGCGTAAACCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGt                                                                                           <  7:349997/149‑1 (MQ=255)
                                                     caGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGt                                                                                           >  8:320241/1‑149 (MQ=255)
                                                     caGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGt                                                                                           >  1:345536/1‑149 (MQ=255)
                                                      aGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGtt                                                                                          <  7:320241/149‑1 (MQ=255)
                                                         gacgaAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGc                                                                                       <  7:159839/149‑1 (MQ=255)
                                                             aaGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCa                                                                                   <  3:120096/149‑1 (MQ=255)
                                                             aaGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCa                                                                                   <  7:235379/149‑1 (MQ=255)
                                                              aGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCa                                                                                   >  6:372456/1‑148 (MQ=255)
                                                              aGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCAt                                                                                  <  1:47389/149‑1 (MQ=255)
                                                                gTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATa                                                                                 >  4:171387/1‑148 (MQ=255)
                                                                gTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATa                                                                                 >  3:28655/1‑148 (MQ=255)
                                                                gTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATAc                                                                                <  3:83934/149‑1 (MQ=255)
                                                                 tGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATAc                                                                                >  6:347444/1‑148 (MQ=255)
                                                                      aTGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGAATATTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAAcgc                                                                           >  1:297132/1‑149 (MQ=255)
                                                                              tCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAAcgcg                                                                          >  1:151839/1‑141 (MQ=255)
                                                                                 ccGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCg                                                               <  1:7135/149‑1 (MQ=255)
                                                                                  cGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCgg                                                              >  2:252288/1‑149 (MQ=255)
                                                                                   gCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCggt                                                             >  4:52087/1‑149 (MQ=255)
                                                                                     tttCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCggtgg                                                           >  6:244509/1‑149 (MQ=255)
                                                                                                      tAGAGGTCGTAGCCGTAAACCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGtctc                                          <  5:347444/149‑1 (MQ=255)
                                                                                                      tAGAGGTCGTAGCCGTAAACCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGtctc                                          <  1:357536/149‑1 (MQ=255)
                                                                                                        gagGTCGTAGCCGTAAACCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCt                                        >  8:212613/1‑149 (MQ=255)
                                                                                                        gagGTCGTAGCCGTAAACCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCt                                        >  5:36401/1‑149 (MQ=255)
                                                                                                              gTAGCCGTAAACCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTgcagg                                   >  4:90054/1‑148 (MQ=255)
                                                                                                                   cGTAAACCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCgg                             >  6:131339/1‑149 (MQ=255)
                                                                                                                    gTAAACCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCggg                            >  2:308288/1‑149 (MQ=255)
                                                                                                                      aaaCCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTg                                                                                               >  3:152711/1‑80 (MQ=255)
                                                                                                                      aaaCCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTg                                                                                               <  4:152711/80‑1 (MQ=255)
                                                                                                                       aaCCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTc                          <  2:219639/148‑1 (MQ=255)
                                                                                                                          ccAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCa                      >  3:66374/1‑149 (MQ=255)
                                                                                                                            aGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATg                    <  3:227376/149‑1 (MQ=255)
                                                                                                                             ggCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATgg                   <  3:178933/149‑1 (MQ=255)
                                                                                                                              gCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGc                  <  5:131339/149‑1 (MQ=255)
                                                                                                                              gCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGc                  >  8:193495/1‑149 (MQ=255)
                                                                                                                                 gg‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGCcgc               >  5:32177/1‑149 (MQ=255)
                                                                                                                                     tCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGCCGCGTc            >  1:62830/1‑149 (MQ=255)
                                                                                                                                       ttCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGCCGCGTCCt          >  4:147924/1‑149 (MQ=255)
                                                                                                                                          cTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGCCGCGTCCTCGa       <  4:66374/149‑1 (MQ=255)
                                                                                                                                          cTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGCCGCGTCCTCGa       >  6:175235/1‑149 (MQ=255)
                                                                                                                                          cTGGCTCTACAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGcc                                                 >  5:233344/1‑107 (MQ=255)
                                                                                                                                          cTGGCTCTACAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGcc                                                 <  6:233344/107‑1 (MQ=255)
                                                                                                                                           tGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGCCGCGTCCTCGAt      >  6:223440/1‑149 (MQ=255)
                                                                                                                                            ggCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGCCGCGTCCTCGATg     <  2:297132/149‑1 (MQ=255)
                                                                                                                                               tGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGCCGCGTCCTCGATGGc   >  2:268210/1‑148 (MQ=255)
                                                                                                                                                 tAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCTGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGCCGCGTCCTCGATGGcc  <  1:308288/147‑1 (MQ=255)
                                                                                                                                                 tAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGCCGCGTCCTCGATGGcc  <  2:345536/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GCCGCCATAGCGCAGGTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGG‑ATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGCCGCGTCCTCGATGGCC  >  NC_002947/4446348‑4446638

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