Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,033,827	Δ1 bp	8.0%	coding (94/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,827	0	C	.	8.0%	222.8 / 13.8	75	coding (94/1212 nt)	dfp	bifunctional 4'‑phosphopantothenoylcysteine decarboxylase/phosphopantothenoylcysteine synthetase
Reads supporting (aligned to +/- strand):  ref base C (33/36);  new base . (3/3);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

CGCCGCAGGCATTTCAAACGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATGGGCCATATCGAGCTGGCCAA  >  NC_002947/6033680‑6033975                                                                                                                                                     |                                                                                                                                                      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                                                                                                                                <  1:660969/150‑1 (MQ=255)                      cgtcgCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCgtcgtc                                                                                                                                <  2:798325/150‑1 (MQ=255)                              gCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCg                                                                                                                        >  2:1610784/1‑150 (MQ=255)                                tctcACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCa                                                                                                                               <  1:1379053/141‑1 (MQ=255)                                tctcACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCa                                                                                                                               >  2:1379053/1‑141 (MQ=255)                                tctcACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCggg                                                                                                                      <  2:642034/150‑1 (MQ=255)                                tctcACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCggg                                                                                                                      <  1:1610784/150‑1 (MQ=255)                                tctcACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCggg                                                                                                                      <  1:175092/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                                                                                  >  2:196885/1‑148 (MQ=255)                                                                     tCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAg                                                                                  <  1:196885/148‑1 (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                                                                                  <  2:1056818/140‑1 (MQ=255)                                                                             gcATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAg                                                                                  >  1:1056818/1‑140 (MQ=255)                                                                             gcATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCgg                                                                        <  2:1024416/150‑1 (MQ=255)                                                                             gcATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCgg                                                                        <  1:1096291/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)                                                                                                 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                                                             <  1:1537459/114‑1 (MQ=255)                                                                                                                            gAGCTGATTCTCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTa                                                             >  2:1537459/1‑114 (MQ=255)                                                                                                                               cTGATTCGCCGT‑CCCCTGGAGCACGGCGCGCAGGTGCGCGTCGTAATGACCC‑GGACGACGCGAGAATGCAGCTCCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATg                      <  1:603615/150‑1 (MQ=255)                                                                                                                                tGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATgg                     <  2:500775/150‑1 (MQ=255)                                                                                                                                  aTTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑a                                                                                                             <  2:1554105/60‑1 (MQ=255)                                                                                                                                  aTTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑a                                                                                                             >  1:1554105/1‑60 (MQ=255)                                                                                                                                    gcGCCGTGCTCCAGGAG‑ACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATtc                                                                                                       <  2:968879/63‑1 (MQ=25)                                                                                                                                    gcGCCGTGCTCCAGGAG‑ACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATtc                                                                                                       >  1:968879/2‑64 (MQ=25)                                                                                                                                               cTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCa                        >  2:693807/1‑133 (MQ=255)                                                                                                                                               cTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCa                        >  2:693868/1‑133 (MQ=255)                                                                                                                                               cTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCa                        <  1:693807/133‑1 (MQ=255)                                                                                                                                               cTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCa                        <  1:693868/133‑1 (MQ=255)                                                                                                                                                tGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAg                                                                                  <  1:1548401/74‑1 (MQ=255)                                                                                                                                                tGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAg                                                                                  >  2:1548401/1‑74 (MQ=255)                                                                                                                                                 ggAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATGGGCCATATCGAGCTGGc     >  2:422838/1‑150 (MQ=255)                                                                                                                                                  gAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATGGGCCATATCGAGCTGGcc    >  1:1075220/1‑150 (MQ=255)                                                                                                                                                    gCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATGGGCCATATCGAGCTGGCCaa  >  1:305111/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                      CGCCGCAGGCATTTCAAACGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATGGGCCATATCGAGCTGGCCAA  >  NC_002947/6033680‑6033975

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