Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,757,038	C→G	7.7%	G26A (GGC→GCC)	PP_3320 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,757,038	0	C	G	7.7%	197.9 / 5.6	65	G26A (GGC→GCC)	PP_3320	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (28/31);  new base G (3/2);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.39e-01

GCGTTGGATGCCAGCACGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGATTTGCAACCTGCGATAGATGTTATG  >  NC_002947/3756892‑3757186                                                                                                                                                   |                                                                                                                                                     gCGTTGGAAGACAGCACGCTATGCTGACATTTCAGTCTTCGGCACATCCCGTTTGTCATAATGTAACAAGCGACCAGAGTTGCTGAGCAGCAGGAATGTACTGAAGCAGCGCCTTATCGTTCATCCAGGCCAGCCTGAGGACTTCGCCg                                                                                                                                                    <  8:328339/149‑1 (MQ=255)   gTTGGATGCCAGCACGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGAt                                                                                                                                                  <  5:233113/149‑1 (MQ=255)       gATGCCAGCACGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCa                                                                                                                                              <  5:36262/149‑1 (MQ=255)       gATGCCAGCACGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCACCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTc                                                                                                                                               <  1:267379/148‑1 (MQ=255)         tGCCAGCACGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCAtg                                                                                                                                            <  3:292232/149‑1 (MQ=255)         tGCCAGCACGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCAtg                                                                                                                                            <  7:199076/149‑1 (MQ=255)                aCGATATGCTGACATTTCCCTCTTCTCCACATCGCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTgg                                                                                                                                     <  4:209005/149‑1 (MQ=255)                 cGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTggg                                                                                                                                    >  1:275856/1‑149 (MQ=255)                      tGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGt                                                                                                                                   <  8:172276/145‑1 (MQ=255)                      tGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGt                                                                                                                                   >  7:172276/1‑145 (MQ=255)                      tGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGAt                                                                                                                                <  5:144692/148‑1 (MQ=255)                       gCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATgg                                                                                                                              >  7:322981/1‑149 (MQ=255)                         tGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCg                                                                                                                            >  7:277551/1‑149 (MQ=255)                           aCATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTgg                                                                                                                                     >  1:340726/1‑138 (MQ=255)                           aCATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTgg                                                                                                                                     <  2:340726/138‑1 (MQ=255)                               ttCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGAt                                                                                                                       >  2:298957/1‑148 (MQ=255)                                         tcCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGc                                                                                                                             >  6:365278/1‑132 (MQ=255)                                         tcCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGc                                                                                                                             <  5:365278/132‑1 (MQ=255)                                           cacaTCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGctgctg                                                                                                                                      <  2:324140/121‑1 (MQ=255)                                           cacaTCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGAGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAggggg                                                                                                          >  8:283591/1‑149 (MQ=255)                                                 ccGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAg                                                                                                    <  3:153569/149‑1 (MQ=255)                                                           tAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGc                                                                                          >  8:36846/1‑149 (MQ=255)                                                           tAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGc                                                                                          >  2:164698/1‑149 (MQ=255)                                                               gTAACAAGCGCCCAGGGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTCGGTGATGGCGGGGATCAGCCCCATGGGGTTACAGGGTATGTGGCGAAg                                                                                      >  4:246428/1‑149 (MQ=255)                                                               gTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAg                                                                                      <  8:192057/149‑1 (MQ=255)                                                                  aCAAGCGCCCAGAGTTGCTGAGCACCATCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGcac                                                                                   <  8:322981/149‑1 (MQ=255)                                                                  aCAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGGGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGcac                                                                                   <  3:243179/149‑1 (MQ=255)                                                                    aaGCGCCCAGAGTTGCTGAGCACCAGAAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGaaca                                                                                  >  4:46490/1‑148 (MQ=255)                                                                      gcgcCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGACGATGTCATGCTGCTGGGTGATGGCGGGGAGCAGCCCCCGGGGGTTCCAGGGTATGTGGAGAAGCACAGCg                                                                               >  2:100596/1‑149 (MQ=255)                                                                           cAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTc                                                                            >  6:142075/1‑147 (MQ=255)                                                                                           cAGCAATGTACTGAAGCAGCGCATTTTCGTTCATCCAGGCAAGAATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCACAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTccc                                                          >  6:88813/1‑149 (MQ=255)                                                                                             gCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGc                                                        <  4:122220/149‑1 (MQ=255)                                                                                                 tGTACTGAAGGAGCTCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCtgt                                                    <  4:160913/149‑1 (MQ=255)                                                                                                    aCTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGcc                                                 >  3:109851/1‑149 (MQ=255)                                                                                                     cTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCt                                                <  2:275856/149‑1 (MQ=255)                                                                                                      tGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGAt                                                                                                                                                  <  6:30332/50‑1 (MQ=255)                                                                                                      tGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGAt                                                                                                                                                  >  5:30332/1‑50 (MQ=255)                                                                                                      tGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCCTGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTg                                               >  5:147212/1‑149 (MQ=255)                                                                                                      tGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTg                                               >  2:77814/1‑149 (MQ=255)                                                                                                      tGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTAGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCCGTCAAGGCGAAGCGGCTTCCCGCCTTTGCCTg                                               >  1:150676/1‑149 (MQ=255)                                                                                                        aaGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCt                                             <  5:142075/149‑1 (MQ=255)                                                                                                             gcgcCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAgg                                        <  8:277551/149‑1 (MQ=255)                                                                                                             gcgcCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAgg                                        <  2:5671/149‑1 (MQ=255)                                                                                                               gcCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGt                                      <  7:36846/149‑1 (MQ=255)                                                                                                                  ttttCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGAGAACCGGCTTCCCGCCTGTGCCTGAACAAGGGt                                      >  7:234526/1‑146 (MQ=255)                                                                                                                        ttCATCGAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCTGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACAt                             <  3:246428/149‑1 (MQ=255)                                                                                                                        ttCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACAt                             >  4:246140/1‑149 (MQ=255)                                                                                                                         tCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACAt                             >  4:246767/1‑148 (MQ=255)                                                                                                                           aTCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGAt                          >  6:213896/1‑149 (MQ=255)                                                                                                                           aTCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGAt                          >  2:294822/1‑149 (MQ=255)                                                                                                                            tCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGAtt                         <  1:164698/149‑1 (MQ=255)                                                                                                                            tCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGAtt                         <  7:283591/149‑1 (MQ=255)                                                                                                                            tCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGAtt                         <  4:109851/149‑1 (MQ=255)                                                                                                                            tCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATCTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGAtt                         <  6:152972/149‑1 (MQ=255)                                                                                                                              tcacccAGCAGCATGACATCGGCAAAGTAATGATGATGGGTGATGGCGGGGATCATCCCCAAGGGGTTCCAGGGTATGTGGCGAAGCACATCGTTCAAGGCGAAGCGGCTTACAGCATGTTCCTGCTCAAGGGTAAGGAACATGAttt                        >  7:56964/5‑148 (MQ=255)                                                                                                                                ggCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGATTTGCa                     <  1:298957/149‑1 (MQ=255)                                                                                                                                 cccAGCAGCATGACATCGGCGAAGTCCTGATGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGATTTGc                      >  7:275881/2‑147 (MQ=38)                                                                                                                                 cccAGCAGCATGACATCGGCGAAGTCCTGATGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGATTTGc                      <  8:275881/146‑1 (MQ=38)                                                                                                                                     gcagcaTGACATCGGCGAAGTCCTGATGCT‑GGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCt                                                <  2:76482/116‑1 (MQ=37)                                                                                                                                     gcagcaTGACATCGGCGAAGTCCTGATGCT‑GGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCt                                                >  1:76482/1‑116 (MQ=37)                                                                                                                                               ttCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGATTTGCAACCTGCGATAGatgt      >  4:189153/1‑149 (MQ=255)                                                                                                                                               ttCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGATTTGCAACCTGCGATAGatgt      >  3:87655/1‑149 (MQ=255)                                                                                                                                               ttCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGATTTGCAACCTGCGATAGat        <  8:267245/147‑1 (MQ=255)                                                                                                                                               ttCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGATTTGCAACCTGCGATAGat        >  7:267245/1‑147 (MQ=255)                                                                                                                                               ttCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCATGCTCAAGGGTAAGGAACATGATTTGCAACCTGCGATAGatgt      >  5:150679/1‑149 (MQ=255)                                                                                                                                                   ccGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGTGTATGTGTCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGATTTGCAACCTGCGATAGatgttatg  <  5:88813/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     GCGTTGGATGCCAGCACGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGATTTGCAACCTGCGATAGATGTTATG  >  NC_002947/3756892‑3757186

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