Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,757,030 G→T 6.4% L29M (CTG→ATG)  PP_3320 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,757,0300GT 6.4% 200.5 / 3.5 62L29M (CTG→ATG) PP_3320hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (29/29);  new base T (2/2);  total (31/31)
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.35e-01

TAGGCGTCAGCGTTGGATGCCAGCACGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGATTTGCAACC  >  NC_002947/3756883‑3757170
                                                                                                                                                    |                                                                                                                                            
taGGCGTCAGCGTTGGATGCCAGCACGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAg                                                                                                                                              <  5:188337/148‑1 (MQ=255)
       cAGCGTTGGATGCCAGCACGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGCAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGGGGGGAAAACAGCGCCTTTTCGTTAATCCAGCCCAGCATCAGGACTTCg                                                                                                                                       >  1:136348/1‑149 (MQ=255)
         gCGTTGGATGCCAGCACGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGcc                                                                                                                                     >  8:7970/1‑148 (MQ=255)
               gATGCCAGCACGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCa                                                                                                                              <  7:7970/149‑1 (MQ=255)
                 tGCCAGCACGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCAtg                                                                                                                            <  1:388860/149‑1 (MQ=255)
                  gCCAGCACGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTG‑AAGCAGCGCATTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCAtgc                                                                                                                           >  8:155268/1‑149 (MQ=255)
                       cACGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGctgctg                                                                                                                      <  1:21785/149‑1 (MQ=255)
                           atatGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCAGCATGACATCGGCGAAGTCCTGATGCTggcct                                                                                                                  <  3:219203/149‑4 (MQ=38)
                              tGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATg                                                                                                               <  4:63004/149‑1 (MQ=255)
                               gCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTG‑AAGCAGCGCCTGTTCGTTCATCCAGGCCTGCATCAGGAGTTCGCCGATGTCATGCAGCGGAGTGATgg                                                                                                              >  4:353364/1‑149 (MQ=255)
                                   aCATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCggg                                                                                                          >  5:158079/1‑149 (MQ=255)
                                   aCATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCgg                                                                                                           <  1:246627/148‑1 (MQ=255)
                                    gaTTTCCCTCTTCTCCTCATCCGGTTTGTCATAACGTAACAATCGCCCAGAGTTGCTGAGCACCATCAATGTACTG‑AAGCAGCGGCTTTTCGTTCATCCAGGCCAGCATCAg                                                                                                                                              <  8:220269/111‑1 (MQ=255)
                                    cATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAg                                                                                                                                              >  7:220269/1‑112 (MQ=255)
                                         cccTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAg                                                                                                    >  8:18360/1‑149 (MQ=255)
                                          ccTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGAt                                                                                                                                  <  7:223377/118‑1 (MQ=255)
                                          ccTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGAt                                                                                                                                  >  8:223377/1‑118 (MQ=255)
                                          ccTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCAtg                                                                                                                            >  6:57298/1‑124 (MQ=255)
                                          ccTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCAtg                                                                                                                            <  5:57298/124‑1 (MQ=255)
                                          ccTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGc                                                                                                   >  4:124518/1‑149 (MQ=255)
                                                   cacaTCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAggggg                                                                                          >  8:154283/1‑149 (MQ=255)
                                                       tCCCTTTTTTCATAACGCAACAACCGCCCATATTTGCTGAGCACCTGAAATGT‑CTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCa                                                                                     <  2:136348/149‑1 (MQ=255)
                                                        cccGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTccg                                                                                     >  8:338092/1‑148 (MQ=255)
                                                         ccGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAg                                                                                    <  7:154283/149‑1 (MQ=255)
                                                         ccGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAg                                                                                    <  6:158079/149‑1 (MQ=255)
                                                           gTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAggg                                                                                  >  6:317157/1‑149 (MQ=255)
                                                             ttGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAg                                                                                    >  3:162927/1‑145 (MQ=255)
                                                             ttGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAg                                                                                    <  4:162927/145‑1 (MQ=255)
                                                                 cATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCAtg                                                                                                                            <  7:318391/101‑1 (MQ=255)
                                                                 cATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCAtg                                                                                                                            >  8:318391/1‑101 (MQ=255)
                                                                    aaCGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCgg                                                                                                           <  6:38739/115‑1 (MQ=255)
                                                                    aaCGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCg                                                                         <  5:130187/149‑1 (MQ=255)
                                                                       gTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAg                                                                      >  2:32704/1‑149 (MQ=255)
                                                                                  ccAGAGTTGCTGAGCACCAGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATc                                                                                                      >  3:235138/1‑106 (MQ=255)
                                                                                  ccAGAGTTGCTGAGCACCAGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATc                                                                                                      <  4:235138/106‑1 (MQ=255)
                                                                                    agagTTGCTGAGCACCAGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTc                                                                                       <  7:319157/119‑1 (MQ=255)
                                                                                    agagTTGCTGAGCACCAGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTc                                                                                       >  8:319157/1‑119 (MQ=255)
                                                                                    agagTTGCTGAGCACCAGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAg                                                         <  7:18360/149‑1 (MQ=255)
                                                                                       gTTGCTGAGCACCAGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGc                                                       >  1:30684/1‑148 (MQ=255)
                                                                                             gAGCACCAGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGaa                                                                       <  4:239334/126‑1 (MQ=255)
                                                                                             gAGCACCAGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGaa                                                                       >  3:239334/1‑126 (MQ=255)
                                                                                              aGCACCAGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCgg                                                <  2:30684/148‑1 (MQ=255)
                                                                                                    aGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGTTCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCg                                                 <  7:27127/141‑1 (MQ=255)
                                                                                                    aGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCg                                                 >  8:27127/1‑141 (MQ=255)
                                                                                                           tACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGc                                  >  7:357068/1‑149 (MQ=255)
                                                                                                             cTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAg                                                         >  1:349301/1‑124 (MQ=255)
                                                                                                             cTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAg                                                         <  2:349301/124‑1 (MQ=255)
                                                                                                              tG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTg                               >  8:152701/1‑149 (MQ=255)
                                                                                                              tG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCA‑GGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAg                                                         <  4:346940/122‑1 (MQ=255)
                                                                                                              tG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCA‑GGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAg                                                         >  3:346940/1‑122 (MQ=255)
                                                                                                                 acgaagcGCCTTTTCGTGAAGCCAGCCCCGCACCAGGACTTCGCCCATGTCAGGCTCCTGGGGGAGGGCGGGCTTCCGCCCCAGGGGGTTCC‑GGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGc                              <  7:338092/143‑1 (MQ=255)
                                                                                                                 aaGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATc                                                                                                      >  7:293440/1‑76 (MQ=255)
                                                                                                                 aaGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATc                                                                                                      <  8:293440/76‑1 (MQ=255)
                                                                                                                    cagcGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCaa                          <  2:283829/149‑1 (MQ=255)
                                                                                                                      gcgcCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAgg                        <  2:30819/149‑1 (MQ=255)
                                                                                                                                 ttCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACAt             >  5:133584/1‑149 (MQ=255)
                                                                                                                                     tCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGAtt         <  1:32704/149‑1 (MQ=255)
                                                                                                                                           ccAGCAGCATGACATCGGCGAAGTCCTGATGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGAtt         <  7:269751/143‑1 (MQ=38)
                                                                                                                                           ccAGCAGCATGACATCGGCGAAGTCCTGATGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGAtt         >  8:269751/1‑143 (MQ=38)
                                                                                                                                           ccAGCAGCATGACATCGGCGAAGTCCTGATGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGATTTGCAAc   >  6:425291/1‑149 (MQ=39)
                                                                                                                                            cAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCt                                <  3:331610/119‑1 (MQ=255)
                                                                                                                                            cAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCt                                >  4:331610/1‑119 (MQ=255)
                                                                                                                                            cAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGATTTGCAAcc  >  4:200181/1‑149 (MQ=255)
                                                                                                                                                 tCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGAtt         <  3:313759/137‑1 (MQ=255)
                                                                                                                                                 tCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGAtt         >  4:313759/1‑137 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
TAGGCGTCAGCGTTGGATGCCAGCACGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTG‑AAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGATTTGCAACC  >  NC_002947/3756883‑3757170

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