Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,804,020	C→G	8.4%	intergenic (‑104/‑40)	PP_3361 ← / → PP_3362	aminotransferase/AsnC family transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,804,020	0	C	G	8.4%	171.9 / 7.4	60	intergenic (‑104/‑40)	PP_3361/PP_3362	aminotransferase/AsnC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (26/28);  new base G (3/2);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.04e-01

TCGCTTGCCGAAAACCGAGGTGCCGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGC‑TGATTGTTTTTCAATAGCT‑TGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCATAGGCCTTTCCCGAAGT  >  NC_002947/3803873‑3804164                                                                                                                                                           |                                                                                                                                                  tCGCTTGCCGAAAACCGAGGTGCCGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGC‑TGATTGTTTTTCAATAGCT‑TGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑cc                                                                                                                                                  <  4:5224/149‑1 (MQ=255)   gCTTGCCGAAAACCGAGGTGCCGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGC‑TGATTGTTTTTCAATAGCT‑TGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACCGGCAT‑A‑‑ATGCCTGTT‑‑GCgc                                                                                                                                                <  3:106131/149‑2 (MQ=255)     ttGCCGAAAACCGAGGTGCCGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGC‑TGATTGTTTTTCAATAGCT‑TGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGAcc                                                                                                                                              >  3:361150/1‑149 (MQ=255)       gCCGAAAACCGAGGTGCCGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGC‑TGATTGTTTTTCAATAGCT‑TGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGAcc                                                                                                                                              <  2:1337/147‑1 (MQ=255)            aaaCCGAGGTGCCGGGTGGGAAAGGTCTGGCATATCATCAGGTTCCGAGCGCGC‑TGATTGTTTTTCAATAGCT‑TGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATaa                                                                                                                                       <  8:271931/149‑1 (MQ=255)             aaCCGAGGTGCCGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGC‑TGATTGTTTTTCAATAGCT‑TGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTC‑‑GGAACCGGCATAAt                                                                                                                                      <  2:114596/149‑1 (MQ=255)               ccGAGGTGCCGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGC‑TGATTGTTTTTCAATAGCT‑TGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGa                                                                                                                                    <  4:165355/149‑1 (MQ=255)                cGAGGTGCCGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGC‑TGATTGTTTTTCAATAGCT‑TGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGAc                                                                                                                                   >  2:110886/1‑149 (MQ=255)                 gAGGTGCCGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGC‑TGATTGTTTTTCAATAGCT‑TGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACg                                                                                                                                  <  2:71989/149‑1 (MQ=255)                               gAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGC‑TGATTGTTTTTCAATAGCT‑TGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAg                                                                                                                    >  1:38050/1‑149 (MQ=255)                                aaaGGTCTGTCATATCATCAGGTTCCGAGCGCGC‑TGATTGTTTTTCAATAGCT‑TGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAgg                                                                                                                   >  6:75152/1‑149 (MQ=255)                                    gTCTGTCATATCATCAGGTTCCGAGCGCGCGTGATTGTTTTTCAATAGCT‑TGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCaa                                                                                                                              >  1:236030/1‑135 (MQ=255)                                    gTCTGTCATATCATCAGGTTCCGAGCGCGC‑TGATTGTTTTTCAATAGCT‑TGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCaa                                                                                                                              <  2:236030/134‑1 (MQ=255)                                           atatCATCAGGTTCCGAGCGCGC‑TGATTGTTTTTCAATAGCT‑TGGCAGGCGCCGCCGAAACGGACACCCTGCATACGATCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGTCAAGCTTATCAAGGGGCTGC‑AATGc                                                                                                       <  6:31231/149‑1 (MQ=255)                                               catcaGTTTTCGATCGGGC‑TGATTGTTTTTCAATAGCTCAGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGt                                                                                                     <  1:329076/148‑1 (MQ=255)                                                 tcaGGTTCCGAGCGCGC‑TGATTGTTTTTCAATAGCT‑TGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGAcc                                                                                                 >  3:292719/1‑149 (MQ=255)                                                 tcaGGTTCCGAGCGCGC‑TGATTGTTTTTCAATAGCT‑TGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGAcc                                                                                                 >  5:112777/1‑149 (MQ=255)                                                  caGGTTCCGAGCGCGC‑TGATTGTTTTTCAATAGCT‑TGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGTTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGACCt                                                                                                >  5:112627/1‑149 (MQ=255)                                                    ggTTCCGAGCGCGC‑TGATTGTTTTTCAATAGCT‑TGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGACCTg                                                                                               >  1:178566/1‑148 (MQ=255)                                                    ggTTCCGAGCGCGC‑TGATTGTTTTTCAATAGCT‑TGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCAA‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGACCTgg                                                                                              <  5:327926/149‑1 (MQ=255)                                                     gTTCCGAGCGCGC‑TGATTGTTTTTCAATAGCT‑TGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGACCTGGa                                                                                             >  4:189768/1‑149 (MQ=255)                                                          gAGCGCGC‑TGATTGTTTTTCAATAGCT‑TGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGACCTGGACGCTa                                                                                        <  1:110886/149‑1 (MQ=255)                                                              gcgc‑TGATTGTTTTTCAATAGCT‑TGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGACCTGGACGCTa                                                                                        >  3:357121/1‑145 (MQ=255)                                                                    gATTGTTTTTCAATAGCT‑TGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGACCTGGACGCTAAAGATcgc                                                                                >  6:47123/1‑148 (MQ=255)                                                                     aTTGTTTTTCAATAGGT‑TGTCATGCGCCGCCGAAACGGACACGATGCATACGGTCGGGCGCAAGAGGCAT‑T‑‑ATGGCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTg                                                                                                              <  2:47799/118‑1 (MQ=255)                                                                     aTTGTTTTTCAATAGCT‑TGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGAACGGCATAATGACGGCAAGCTTATCAAGGGGCTg                                                                                                              >  1:47799/1‑118 (MQ=255)                                                                      ttGTTTTTCAATAGCT‑TGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAAAAGCGTGACCTGGACGCTAAAGATCGCga                                                                              <  1:16708/149‑1 (MQ=255)                                                                      ttGTTTTTCAATAGCT‑TGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGACCTGGACGCTAAAGATCGCga                                                                              <  2:38050/148‑1 (MQ=255)                                                                              cAATAGCT‑TGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTg                                                                      <  6:112777/148‑1 (MQ=255)                                                                               aaTAGCT‑TGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGa                                                                     <  3:282827/148‑1 (MQ=255)                                                                                aTAGCT‑TGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAgg                                                                   >  4:104037/1‑149 (MQ=255)                                                                                   tgt‑tGTCTTGCGCCGCCGAAACGGACACCCTGCTTAGGGTCGGGCGCAACAGGCATCTAGATGCCGGTTCCCCCATCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTg                                                                      <  1:123157/146‑1 (MQ=255)                                                                                   gCT‑TGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTAc                                                                <  4:292719/149‑1 (MQ=255)                                                                                     t‑tGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACtt                                                              >  1:358032/1‑149 (MQ=255)                                                                                         tCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAg                                                                                                                    <  2:265982/93‑1 (MQ=255)                                                                                         tCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAg                                                                                                                    >  1:265982/1‑93 (MQ=255)                                                                                           aTGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑ACGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGaa                                                         >  6:252274/1‑149 (MQ=255)                                                                                              cgccgcCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAgg                                                       <  2:358032/148‑1 (MQ=255)                                                                                               gccgccGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGAc                                                     <  4:361150/149‑1 (MQ=255)                                                                                                 cgccgAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGAcgc                                                   >  1:13511/1‑149 (MQ=255)                                                                                                   ccgAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGAcgcgc                                                 <  5:47123/149‑1 (MQ=255)                                                                                                        aCGGACACCCTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTAtt                                            >  1:283778/1‑149 (MQ=255)                                                                                                             cacCCTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCt                                       <  3:104037/149‑1 (MQ=255)                                                                                                                ccTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGaa                                                         <  7:188246/128‑1 (MQ=255)                                                                                                                ccTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTc                                      >  7:159658/1‑147 (MQ=255)                                                                                                                ccTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCaa                                    <  8:84883/149‑1 (MQ=255)                                                                                                                ccTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGACATGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGaa                                                         >  8:188246/1‑128 (MQ=255)                                                                                                                ccTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCGTATCAAGGGGCTGC‑AATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTc                                      >  1:301291/1‑147 (MQ=255)                                                                                                                         ggTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAgc                           <  3:189768/149‑1 (MQ=255)                                                                                                                              ggCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAgccggcc                      >  2:234064/1‑149 (MQ=255)                                                                                                                              ggCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAgccggcc                      >  8:188527/1‑149 (MQ=255)                                                                                                                                cgcAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAgccggccg                     >  5:120487/1‑148 (MQ=255)                                                                                                                                      cgtcAT‑T‑‑ATGCCGGTC‑‑GGAACCGGCATAATAACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCATAGGc              >  7:358922/4‑149 (MQ=255)                                                                                                                                       ggCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCATAGGc              >  7:160967/1‑148 (MQ=255)                                                                                                                                       ggCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCATAGGc              <  8:160967/148‑1 (MQ=255)                                                                                                                                        gCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCATAGGCCt            >  6:183124/1‑149 (MQ=255)                                                                                                                                         cAT‑T‑‑ATGCCGGTC‑‑GGAACCGGCATAATGCCGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCATAGGCCtt           >  2:218816/1‑149 (MQ=255)                                                                                                                                           t‑t‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCATAGGCCTTTc         <  6:120487/149‑1 (MQ=255)                                                                                                                                             t‑‑ATGCCGGTC‑‑GGAACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCATAGGCCTTTcc        >  4:182834/1‑149 (MQ=255)                                                                                                                                                     ggTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCATAGGCCTTTCCCGAAGt  <  2:178566/149‑1 (MQ=255)                                                                                                                                                     ggTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGACCTGGACGCTAAAGATCCCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCATAGGCCTTTCCCGAAGt  <  7:188527/149‑1 (MQ=255)                                                                                                                                                           |                                                                                                                                                  TCGCTTGCCGAAAACCGAGGTGCCGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGC‑TGATTGTTTTTCAATAGCT‑TGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCAT‑T‑‑ATGCCGGTT‑‑CCGACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGC‑AATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCATAGGCCTTTCCCGAAGT  >  NC_002947/3803873‑3804164

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