Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,033,817:1	+G	8.8%	coding (84/1212 nt)	dfp →	bifunctional 4'‑phosphopantothenoylcysteine decarboxylase/phosphopantothenoylcysteine synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,033,817	1	.	G	8.8%	258.0 / 18.5	68	coding (84/1212 nt)	dfp	bifunctional 4'‑phosphopantothenoylcysteine decarboxylase/phosphopantothenoylcysteine synthetase
Reads supporting (aligned to +/- strand):  ref base . (28/34);  new base G (3/3);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

CATGCGACGTTCCGCCGCAGGCATTTCAAACGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATGG  >  NC_002947/6033668‑6033956                                                                                                                                                       |                                                                                                                                             cATGCGACGTTCCGCCGCAGGCATTTCAAACGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGt                                                                                                                                               >  2:611264/1‑150 (MQ=255) cATGCGACGTTCCGCCGCAGGCATTTCAAACGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGt                                                                                                                                               >  2:408737/1‑150 (MQ=255)  aTGCGACGTTCCGCCGCAGGCATTTCAAACGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGt‑c                                                                                                                                             <  1:432501/150‑1 (MQ=255)             cgccgcAGGCATTTCAAACGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCAc                                                                                                                                  >  1:169250/1‑150 (MQ=255)                   aGGCATTTCAAACGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CCCCTGGAGCACGgcgcg                                                                                                                            >  2:175092/1‑150 (MQ=255)                     gCATTTCAAACGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGTGCTCCAGGAG‑ACGgcgaat                                                                                                                           >  2:671911/1‑146 (MQ=255)                     gCATTTCAAACGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGTGCTCCAGGAG‑ACGgcgaat                                                                                                                           <  1:671911/149‑4 (MQ=255)                     gCATTTCAAACGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCa                                                                                                                          <  2:169250/150‑1 (MQ=255)                       aTTTCAAACGTCGTCGCCCGCTCTCACTGGAAAAAGGCTAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATGCGCCGTGCTCCAGGAG‑ACGgcgaac                                                                                                                           <  2:1126264/147‑4 (MQ=255)                       aTTTCAAACGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATGCGCCGTGCTCCAGGAG‑ACGgcgaat                                                                                                                           >  1:1126264/1‑144 (MQ=255)                        tttCAAACGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGt                                                                                                                       >  1:642034/1‑150 (MQ=255)                        tttCAAACGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGt                                                                                                                       >  2:1080830/1‑150 (MQ=255)                            aaaCGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTgcgc                                                                                                                   >  2:660969/1‑150 (MQ=255)                            aaaCGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTAGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTgcgc                                                                                                                   >  1:627936/1‑150 (MQ=255)                              aCGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGcgt                                                                                                                 >  1:798325/1‑150 (MQ=255)                                  cgtcgCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCgtcgtc                                                                                                             <  2:798325/150‑1 (MQ=255)                                  cgtcgCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCgtcgtc                                                                                                             <  1:660969/150‑1 (MQ=255)                                          gCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCg                                                                                                     >  2:1610784/1‑150 (MQ=255)                                            tctcACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCa                                                                                                            >  2:1379053/1‑141 (MQ=255)                                            tctcACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCa                                                                                                            <  1:1379053/141‑1 (MQ=255)                                            tctcACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCggg                                                                                                   <  2:642034/150‑1 (MQ=255)                                            tctcACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCggg                                                                                                   <  1:175092/150‑1 (MQ=255)                                            tctcACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCggg                                                                                                   <  1:1610784/150‑1 (MQ=255)                                              tcACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGggt                                                                                                 <  1:1080830/150‑1 (MQ=255)                                               cacggtaaaaaGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGggtg                                                                                                <  2:1579029/144‑1 (MQ=255)                                                  tGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGgtggtg                                                                                             >  1:500775/1‑150 (MQ=255)                                                  tGGAAAAAAGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGgtggtg                                                                                             <  1:1040201/150‑1 (MQ=255)                                                   ggAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGc                                                                                            <  1:1238729/150‑1 (MQ=255)                                                    gAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑a                                                                                          >  2:603615/1‑150 (MQ=255)                                                      aaaaGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCggg                                                                                                   >  1:8499/1‑140 (MQ=255)                                                      aaaaGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCggg                                                                                                   <  2:8499/140‑1 (MQ=255)                                                      aaaaGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGa                                                                                        >  2:1096291/1‑150 (MQ=255)                                                       aaaGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTTGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGaa                                                                                       <  2:1186726/150‑1 (MQ=255)                                                       aaaGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGggtg                                                                                                >  1:665586/1‑142 (MQ=255)                                                       aaaGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGggtg                                                                                                <  2:665586/142‑1 (MQ=255)                                                         aGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAAtt                                                                                     <  1:1068958/150‑1 (MQ=255)                                                          ggCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATtc                                                                                    <  2:1209922/150‑1 (MQ=255)                                                           gCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATtca                                                                                   >  1:805075/1‑150 (MQ=255)                                                                              gTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCa                                                                <  2:884976/150‑1 (MQ=255)                                                                              gTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCa                                                                <  2:1131343/150‑1 (MQ=255)                                                                                 tCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAg                                                               <  1:196885/148‑1 (MQ=255)                                                                                 tCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAg                                                               >  2:196885/1‑148 (MQ=255)                                                                                   gCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGACGCGCACCTGCGcgccgt                                                                                                              <  2:1563914/100‑1 (MQ=255)                                                                                   gCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGACGCGCACCTGCGcgccgt                                                                                                              >  1:1563914/1‑100 (MQ=255)                                                                                   gCAAGCGCATCATTCTCGCCGTGGGAGGCCGAATTGCCGCCTCCGAAAGCGCAGAGCTGATTCGCCGT‑CCCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGGATTCATCACCCCGCTGACCCTGCGGgcgc                                                           >  2:525706/1‑150 (MQ=255)                                                                                         gcATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAg                                                               >  1:1056818/1‑140 (MQ=255)                                                                                         gcATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAg                                                               <  2:1056818/140‑1 (MQ=255)                                                                                         gcATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCgg                                                     <  1:1096291/150‑1 (MQ=255)                                                                                         gcATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCgg                                                     <  2:1024416/150‑1 (MQ=255)                                                                                          cATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGAGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCggg                                                    <  2:805075/150‑1 (MQ=255)                                                                                           aTCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGc                                                   <  2:132726/150‑1 (MQ=255)                                                                                              gTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCAc                                                >  1:340225/1‑150 (MQ=255)                                                                                               ttCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATTACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCAcc                                               <  1:611264/150‑1 (MQ=255)                                                                                                  tcGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCg                                            >  2:1608082/1‑150 (MQ=255)                                                                                                      cGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTg                                                                                                                                        <  1:557234/55‑1 (MQ=255)                                                                                                      cGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTg                                                                                                                                        >  2:557234/1‑55 (MQ=255)                                                                                                             ggcggcATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGAc                                 <  2:340225/150‑1 (MQ=255)                                                                                                                ggcATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGTGCCACCCCGTACACATGGACCTg                              >  1:1376674/1‑150 (MQ=255)                                                                                                                 gcATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGc                             <  1:408737/150‑1 (MQ=255)                                                                                                                      gccgccTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGATGTCGGGCCACCCCGTACACATGGACCTGCTCGAc                        >  2:877428/1‑150 (MQ=255)                                                                                                                                     gcCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCt         >  1:1169672/1‑150 (MQ=255)                                                                                                                                      ccGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTg        >  2:602510/1‑150 (MQ=255)                                                                                                                                       cGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCCGc       <  1:877428/150‑1 (MQ=255)                                                                                                                                        gAGCTGATTCTCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTa                                          >  2:1537459/1‑114 (MQ=255)                                                                                                                                        gAGCTGATTCTCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTa                                          <  1:1537459/114‑1 (MQ=255)                                                                                                                                           cTGATTCGCCGT‑CCCCTGGAGCACGGCGCGCAGGTGCGCGTCGTAATGACCC‑GGACGACGCGAGAATGCAGCTCCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATg   <  1:603615/150‑1 (MQ=255)                                                                                                                                            tGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATgg  <  2:500775/150‑1 (MQ=255)                                                                                                                                              aTTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑a                                                                                          >  1:1554105/1‑60 (MQ=255)                                                                                                                                              aTTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑a                                                                                          <  2:1554105/60‑1 (MQ=255)                                                                                                                                                gcGCCGTGCTCCAGGAG‑ACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATtc                                                                                    >  1:968879/2‑64 (MQ=25)                                                                                                                                                gcGCCGTGCTCCAGGAG‑ACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATtc                                                                                    <  2:968879/63‑1 (MQ=25)                                                                                                                                                       |                                                                                                                                             CATGCGACGTTCCGCCGCAGGCATTTCAAACGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATGG  >  NC_002947/6033668‑6033956

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑