Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,938,223	T→A	5.4%	F75L (TTT→TTA) ‡	PP_3471 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,938,223	0	T	A	5.4%	229.5 / 2.5	74	F75L (TTT→TTA) ‡	PP_3471	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (30/40);  new base A (2/2);  total (32/42)
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.30e-01

ATAAAACATCCGATGGATATTTACTCGGGAATGGTCGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAAC  >  NC_002947/3938075‑3938370                                                                                                                                                     |                                                                                                                                                    aTAAATCATCCGATGGATATTTACTCGGGAATGGTCGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACttt                                                                                                                                                     <  5:240552/149‑1 (MQ=255) aTAAAACATCCGATGGATATTTACTCGGGAATGGTCGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACttt                                                                                                                                                     >  8:86290/1‑149 (MQ=255) aTAAAACATCCGATGGATATTTACTCGGGAATGGTCGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACttt                                                                                                                                                     <  5:364420/149‑1 (MQ=255)  tAAAACATCCGATGGATATTTACTCGGGAATGGTCGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTttg                                                                                                                                                    <  4:61436/149‑1 (MQ=255)    aaaCATCCGATGGATATTTACTCGGGAATGGTCGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTttgt                                                                                                                                                   <  3:208128/148‑1 (MQ=255)        aTCCGATGGATATTTACTCGGGAATGGTCGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATc                                                                                                                                              >  2:185413/1‑149 (MQ=255)             aTGGATATTTAGTCGGGTATGGTCGGCGCGGATACAAGCCGCGCTTGTTGCGTTTTTGTTCGCCTGACAATTTAAGTCCATTTCAATCAGGCGTCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCATCTTTGTTGATCCAAGt                                                                                                                                         <  8:299219/149‑1 (MQ=255)             aTGGATATTTACTCGGGAATGGTCGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGt                                                                                                                                         <  6:223677/149‑1 (MQ=255)             aTGGATATTTACTCGGGAATGGTCGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACAAAGTTGATcggatc                                                                                                                                         <  6:13338/149‑6 (MQ=255)                gATATTTACTCGGGAATGGTCGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGc                                                                                                                                      <  8:252142/149‑1 (MQ=255)                    tttACTCGGGAATGGTCGGCGCTGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCAtt                                                                                                                                  >  7:55951/1‑149 (MQ=255)                      tACTCGGGAATGGTCGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTc                                                                                                                                 <  1:185413/148‑1 (MQ=255)                        cTCGGGAATGGTCGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCg                                                                                                                               <  7:86290/148‑1 (MQ=255)                        cTCGGGAATGGTCGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCg                                                                                                                               <  2:215938/148‑1 (MQ=255)                         tCGGGAATGGTCGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGcc                                                                                                                                     >  6:2650/1‑141 (MQ=255)                         tCGGGAATGGTCGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGcc                                                                                                                                     <  5:2650/141‑1 (MQ=255)                         tCGGGAATGGTCGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGc                                                                                                                             >  4:11988/1‑149 (MQ=255)                         tCGGGAATGGTCGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGc                                                                                                                             >  8:88741/1‑149 (MQ=255)                           gggAATGGTCGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCggn                                                                                                                             <  2:309782/147‑2 (MQ=255)                                tGGTCGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGc                                                                                                                      <  4:104053/149‑1 (MQ=255)                                   tCGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCAc                                                                                                                   >  7:82046/1‑149 (MQ=255)                                     ggCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATcc                                                                                                                                             >  8:156111/1‑121 (MQ=255)                                       cgcgGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTAc                                                                                                                >  8:283241/1‑148 (MQ=255)                                        gcgGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGAGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACaa                                                                                                              <  5:369383/149‑1 (MQ=255)                                           gATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAAt                                                                                                            >  8:251969/1‑148 (MQ=255)                                           gATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACAAAGTTGATcgg                                                                                                                                            <  8:249719/116‑3 (MQ=255)                                           gATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACAAAGTTGATcgg                                                                                                                                            >  7:249719/1‑114 (MQ=255)                                             tACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCt                                                                                                         <  4:295707/149‑1 (MQ=255)                                               cAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTgg                                                                                                       <  7:88741/149‑1 (MQ=255)                                               cAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTACATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTgg                                                                                                       <  8:100852/149‑1 (MQ=255)                                                     gcgcTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTgg                                                                                                                                       <  2:97977/111‑1 (MQ=255)                                                     gcgcTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTgg                                                                                                                                       >  1:97977/1‑111 (MQ=255)                                                         ttGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCaaa                                                                                              >  4:156507/1‑148 (MQ=255)                                                            aTGCGTTTTAGTTCGCCGGACAATTTAAGTCCATTTAAATCAGGCGGCATCAACAGCTACGCCTACTGCAAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGCCGATCTGCCAATACAAATTCTGACCGCGCAAAAAt                                                                                           >  1:265773/1‑148 (MQ=255)                                                             tGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTAGCCGCGCAAAAATTa                                                                                         >  1:361847/1‑149 (MQ=255)                                                             tGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTAGCCGCGCAAAAATTa                                                                                         <  2:361847/149‑1 (MQ=255)                                                              gCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCcgc                                                                                                   >  6:167258/1‑138 (MQ=255)                                                              gCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCcgc                                                                                                   <  5:167258/138‑1 (MQ=255)                                                               cggttttGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCGAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTAcc                                                                                       <  7:291318/146‑1 (MQ=255)                                                               cGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTAcc                                                                                       <  4:92916/149‑1 (MQ=255)                                                               cGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTAcc                                                                                       >  1:208467/1‑149 (MQ=255)                                                                      gTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATGAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTc                                                                                 <  5:320221/148‑1 (MQ=255)                                                                           ccGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAg                                                                           <  4:13880/149‑1 (MQ=255)                                                                             ggAAAATTTAAGTCCATTTCAATCAGGGGGCATCAACATCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGc                                                                         <  2:265773/149‑1 (MQ=255)                                                                               aCAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGc                                                                         >  1:166894/1‑147 (MQ=255)                                                                                  aTTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGa                                                                    <  8:82046/149‑1 (MQ=255)                                                                                   tttAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGa                                                                    >  5:374883/1‑148 (MQ=255)                                                                                       aGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAAtt                                                                <  7:251969/148‑1 (MQ=255)                                                                                        gTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAAttt                                                               >  4:274949/1‑148 (MQ=255)                                                                                         tCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTaa                                                             <  6:105282/149‑1 (MQ=255)                                                                                           cATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCACAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAag                                                            <  7:270449/148‑1 (MQ=255)                                                                                             tttCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAagag                                                          >  8:59704/1‑148 (MQ=255)                                                                                                 aaTCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATcc                                                     <  2:55196/149‑1 (MQ=255)                                                                                                   tCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCaa                                                   <  7:283241/149‑1 (MQ=255)                                                                                                     aGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAAc                                                  <  8:55951/148‑1 (MQ=255)                                                                                                          gcATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCgg                                            >  2:137371/1‑149 (MQ=255)                                                                                                               aaCAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATc                                       >  5:157082/1‑149 (MQ=255)                                                                                                                aCAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATcc                                      <  7:16386/149‑1 (MQ=255)                                                                                                                aCAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATcc                                      >  2:129974/1‑149 (MQ=255)                                                                                                                aCAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATcc                                      <  7:22007/149‑1 (MQ=255)                                                                                                                       cGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTaa                               >  5:184638/1‑149 (MQ=255)                                                                                                                          cTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGt                            >  4:61323/1‑149 (MQ=255)                                                                                                                             cTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGttt                          <  1:137371/148‑1 (MQ=255)                                                                                                                             cTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGtt                           <  1:129974/147‑1 (MQ=255)                                                                                                                                ccAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCTGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTtgat                       <  7:312229/148‑1 (MQ=255)                                                                                                                                ccAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTtga                        <  3:274949/147‑1 (MQ=255)                                                                                                                                ccAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCAAATTAATCTTtgat                       >  1:67297/1‑148 (MQ=255)                                                                                                                                    ggAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGAtt                   >  3:292942/1‑148 (MQ=255)                                                                                                                                       gaTCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACa                >  2:171610/1‑148 (MQ=255)                                                                                                                                        cttggATCAACAAAGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACAt               >  6:86999/5‑148 (MQ=255)                                                                                                                                           cGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTtgatga                     <  8:4587/139‑1 (MQ=255)                                                                                                                                           cGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTtgatga                     >  7:4587/1‑139 (MQ=255)                                                                                                                                           cGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATaaaa           <  2:208467/149‑1 (MQ=255)                                                                                                                                           cGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATaaaa           <  6:374883/149‑1 (MQ=255)                                                                                                                                            gATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTTGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAaga                                                           >  4:266619/1‑100 (MQ=255)                                                                                                                                            gATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTTGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAaga                                                           <  3:266619/100‑1 (MQ=255)                                                                                                                                              tCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTa        <  7:59704/149‑1 (MQ=255)                                                                                                                                                 aCTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTAtt      >  4:8182/1‑148 (MQ=255)                                                                                                                                                  cTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCTCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTaa    <  5:327583/149‑1 (MQ=255)                                                                                                                                                   tttGTTGATCCAAGTGGCGATTACGGCGTTCTGCCACTACAAATTCTGGCAGCGCAAAAATTACCGACTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAAACTGCGGTCATCCCAATTAATGATGGATGATTACATAAAACtaataa    >  4:314697/1‑148 (MQ=255)                                                                                                                                                   tttGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTaaa   <  3:11988/149‑1 (MQ=255)                                                                                                                                                     tgttgATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAaac  <  3:156507/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    ATAAAACATCCGATGGATATTTACTCGGGAATGGTCGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAAC  >  NC_002947/3938075‑3938370

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