Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,338,632	G→A	7.8%	noncoding (356/395 nt)	PP_mr53 →	CrcZ ncRNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,338,632	0	G	A	7.8%	185.9 / 4.5	64	noncoding (356/395 nt)	PP_mr53	CrcZ ncRNA
Reads supporting (aligned to +/- strand):  ref base G (33/26);  new base A (3/2);  total (36/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.08e-01

ATCAGCGGCCAAAGCAATCCGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTCACACAGTGCATCCCATG  >  NC_002947/5338485‑5338775                                                                                                                                                    |                                                                                                                                                aTCAGCGGCCAAAGCAATCCGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAgg                                                                                                                                                <  7:148602/149‑1 (MQ=255)    aGCGGCCAAAGCAATCCGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTc                                                                                                                                             <  1:111295/149‑1 (MQ=255)     gCGGCCAAAGCAATCCGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTcc                                                                                                                                            >  7:64642/1‑149 (MQ=255)        gCCAAAGCAATCCGTTTGCTCTTTACCCCCGGTTTGGGGTTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCCTAAGCTcccca                                                                                                                                         <  3:272046/149‑2 (MQ=255)         ccAAAGCAATCCGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCCTAAGCTccccaa                                                                                                                                        <  3:48833/149‑3 (MQ=255)           aaaGCAATCCGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGt                                                                                                                                      <  8:64642/149‑1 (MQ=255)           aaaGCAATCCGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGt                                                                                                                                      >  2:192610/1‑149 (MQ=255)            aaGCAATCCGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTc                                                                                                                                     <  3:124898/149‑1 (MQ=255)                aaTCCGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGttt                                                                                                                                 <  5:169075/149‑1 (MQ=255)                     gTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACCATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTAcc                                                                                                                            >  8:144341/1‑149 (MQ=255)                      tttGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTAccc                                                                                                                           <  6:369712/149‑1 (MQ=255)                                aCCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTggcag                                                                                                                  <  6:70260/148‑1 (MQ=255)                                  ccccGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTggcaggca                                                                                                               >  1:8669/1‑149 (MQ=255)                                  ccccGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTggcaggca                                                                                                               >  2:204980/1‑149 (MQ=255)                                   cccGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTc                                                                                                                                <  4:63855/131‑1 (MQ=255)                                   cccGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTc                                                                                                                                >  3:63855/1‑131 (MQ=255)                                    ccGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATc                                                                                                             >  8:292117/1‑149 (MQ=255)                                           gggggTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTa                                                                                                       <  1:192610/148‑1 (MQ=255)                                                   tcCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCATa                                                                                               >  1:346092/1‑148 (MQ=255)                                                     cacaAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTaca                                                                                             <  8:131020/148‑1 (MQ=255)                                                     cacaAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGAACGCAACGACTTCTTGGGGAGATTAGGCTCCACTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTacac                                                                                            >  1:86976/1‑149 (MQ=255)                                                           cTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTcc                                                                                      >  6:343439/1‑149 (MQ=255)                                                            ttCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTa                                                                                               >  4:56725/1‑139 (MQ=255)                                                            ttCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTa                                                                                               <  3:56725/139‑1 (MQ=255)                                                             tCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTAc                                                                                                      >  8:11952/1‑131 (MQ=255)                                                              cGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTAcc                                                                                                                            <  3:32606/108‑1 (MQ=255)                                                              cGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTAcc                                                                                                                            >  4:32606/1‑108 (MQ=255)                                                               gagaTTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTa                                                                                               <  4:184877/136‑1 (MQ=255)                                                               gagaTTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTa                                                                                               >  3:184877/1‑136 (MQ=255)                                                                agaTTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAg                                                                                 >  4:292032/1‑149 (MQ=255)                                                                     tgtgGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTaaa                                                                             <  5:343439/148‑1 (MQ=255)                                                                       tgGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGc                                                                                                  >  7:74138/1‑125 (MQ=255)                                                                       tgGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGc                                                                                                  <  8:74138/125‑1 (MQ=255)                                                                         gACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGc                                                                          >  8:242666/1‑147 (MQ=255)                                                                              ggCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATc                                                                   >  4:367001/1‑149 (MQ=255)                                                                              ggCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATc                                                                   >  6:239363/1‑149 (MQ=255)                                                                                  cctatctaaaaCAAGAAGCCCAGCAAAAAACCAATACTAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCATTGTCGTTTCTACCCTTCTGGCATGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATcccc                                                                <  6:226833/141‑1 (MQ=255)                                                                                         aaaaCAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGCGCTTAGGGTCCCCTTGTCGTTTCTACGGTTCTGGCAGGGATGCTGCTAGCCGCCTAGACCATTCCCCGAGTTAATGTTAGAATCGCCGCCcatga                                                        >  4:290528/1‑147 (MQ=255)                                                                                         aaaaCAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTcccc                                                                                    >  3:143551/1‑121 (MQ=255)                                                                                         aaaaCAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTcccc                                                                                    <  4:143551/121‑1 (MQ=255)                                                                                          aaaCAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACACGCCTACACCATTCCCCGAGAAAATGCTAGAATcccc                                                                >  5:226833/1‑140 (MQ=255)                                                                                                gaagCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGc                                                  >  7:210767/1‑148 (MQ=255)                                                                                                     ccAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATtct                                            >  4:78998/1‑149 (MQ=255)                                                                                                        gCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCt                                         <  7:292117/149‑1 (MQ=255)                                                                                                            aaaaaCAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTa                                                                        <  6:2646/114‑1 (MQ=255)                                                                                                            aaaaaCAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTa                                                                        >  5:2646/1‑114 (MQ=255)                                                                                                              aaaCAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGc                                                              >  4:68285/1‑122 (MQ=255)                                                                                                              aaaCAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGc                                                              <  3:68285/122‑1 (MQ=255)                                                                                                                       gagCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACa                          >  7:183965/1‑149 (MQ=255)                                                                                                                           aCGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTcc                      >  4:102128/1‑149 (MQ=255)                                                                                                                           aCGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTcc                      >  4:96747/1‑149 (MQ=255)                                                                                                                              cAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTc                   <  1:204980/149‑1 (MQ=255)                                                                                                                               aaCGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAg                                                                                 >  6:354469/1‑86 (MQ=255)                                                                                                                               aaCGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAg                                                                                 <  5:354469/86‑1 (MQ=255)                                                                                                                                aCGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTcac                 <  7:144341/149‑1 (MQ=255)                                                                                                                                   aCTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTcacaca              <  3:78998/149‑1 (MQ=255)                                                                                                                                     ttcttGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTCACACAGt            >  6:134890/1‑149 (MQ=255)                                                                                                                                     ttcttGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTCACACAGt            >  6:284575/1‑149 (MQ=255)                                                                                                                                     gacaaggggAGCCTAAGCTCCCCTTGTCGTTTCTGCCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTCACACAGt            >  6:374221/6‑149 (MQ=255)                                                                                                                                     gacaaggggAGCCTAAGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTCACACAGt            >  4:155448/6‑149 (MQ=255)                                                                                                                                       caaggggAGCCTAAGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTCACACAGTg           >  2:137111/4‑148 (MQ=255)                                                                                                                                           gggAGCTTAGGCGCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTCACACAGTGCATc       <  3:274585/148‑1 (MQ=255)                                                                                                                                               gCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTCACACAGTGCATCCCATg  <  3:292032/149‑1 (MQ=255)                                                                                                                                               gCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTCACACAGTGCATCCCATg  >  1:305275/1‑149 (MQ=255)                                                                                                                                               gCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTCACACAGTGCATCCCATg  <  1:151707/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                ATCAGCGGCCAAAGCAATCCGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTCACACAGTGCATCCCATG  >  NC_002947/5338485‑5338775

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