Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,757,046	A→C	7.8%	H23Q (CAT→CAG)	PP_3320 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,757,046	0	A	C	7.8%	189.0 / 3.3	61	H23Q (CAT→CAG)	PP_3320	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (28/28);  new base C (3/2);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGCCAGCACGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGATTTGCAACCTGCGATAGATGTTATGTTAT  >  NC_002947/3756898‑3757190                                                                                                                                                     |                                                                                                                                                 gATGCCAGCACGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCa                                                                                                                                                  <  5:36262/149‑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                                                                                                                                       >  7:172276/1‑145 (MQ=255)                tGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGt                                                                                                                                       <  8:172276/145‑1 (MQ=255)                tGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGAt                                                                                                                                    <  5:144692/148‑1 (MQ=255)                 gCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATgg                                                                                                                                  >  7:322981/1‑149 (MQ=255)                   tGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCg                                                                                                                                >  7:277551/1‑149 (MQ=255)                     aCATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTgg                                                                                                                                         <  2:340726/138‑1 (MQ=255)                     aCATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTgg                                                                                                                                         >  1:340726/1‑138 (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                                                                                              >  2:164698/1‑149 (MQ=255)                                                     tAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGc                                                                                              >  8:36846/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)                                                                                                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                                            <  2:5671/149‑1 (MQ=255)                                                                                                       gcgcCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAgg                                            <  8:277551/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                              >  2:294822/1‑149 (MQ=255)                                                                                                                     aTCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGAt                              >  6:213896/1‑149 (MQ=255)                                                                                                                      tCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGAtt                             <  7:283591/149‑1 (MQ=255)                                                                                                                      tCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGAtt                             <  1:164698/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          >  3:87655/1‑149 (MQ=255)                                                                                                                                         ttCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGATTTGCAACCTGCGATAGatgt          >  4:189153/1‑149 (MQ=255)                                                                                                                                         ttCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGATTTGCAACCTGCGATAGat            >  7:267245/1‑147 (MQ=255)                                                                                                                                         ttCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGATTTGCAACCTGCGATAGat            <  8:267245/147‑1 (MQ=255)                                                                                                                                         ttCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCATGCTCAAGGGTAAGGAACATGATTTGCAACCTGCGATAGatgt          >  5:150679/1‑149 (MQ=255)                                                                                                                                             ccGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGTGTATGTGTCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGATTTGCAACCTGCGATAGatgttatg      <  5:88813/149‑1 (MQ=255)                                                                                                                                                 tGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGATTTGCAACCTGCGATAGATgttatgttat  <  3:246767/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 GATGCCAGCACGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGATTTGCAACCTGCGATAGATGTTATGTTAT  >  NC_002947/3756898‑3757190

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