Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	660,589	G→A	9.2%	V545V (GTG→GTA)	speA →	biosynthetic arginine decarboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	660,589	0	G	A	9.2%	182.1 / 6.7	65	V545V (GTG→GTA)	speA	biosynthetic arginine decarboxylase
Reads supporting (aligned to +/- strand):  ref base G (26/33);  new base A (2/4);  total (28/37)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.65e-01

GCCGATGCGGCGTGCGGTATTGCAGGACCTGACTTGCGATTCCGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATCGAAACCCACGACACCATCGAAGACATGCTGCGCTACGTGCA  >  NC_002947/660442‑660735                                                                                                                                                    |                                                                                                                                                   gCCGATGCGGCGTGCGGTATTGCAGGACCTGACTTGCGATTCCGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCtggtgg                                                                                                                                                   <  2:151436/149‑1 (MQ=255)         ggcgTGCGGTATTGCAGGACCTGACTTGCGATTCCGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTAc                                                                                                                                           <  8:145655/149‑1 (MQ=255)           cgTGCGGTATTGCAGGACCTGACTTGCGATTCCGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCa                                                                                                                                         <  1:64956/149‑1 (MQ=255)             tGCGGTATTGCAGGACCTGACTTGCGATTCCGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTAGGCACCCACCAgg                                                                                                                                       <  3:322918/149‑1 (MQ=255)               cGGTATTGCAGGACCTGACTTGCGATTCCGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTAGGCACCCACCAGGaa                                                                                                                                     <  8:38880/149‑1 (MQ=255)               cGGTATTGCAGGACCTGACTTGCGATTCCGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTAGGCACCCACCAGGaa                                                                                                                                     <  4:140795/149‑1 (MQ=255)                ggTATTGCAGGACCTGACTTGCGATTCCGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTAGGCACCCACCAGGaac                                                                                                                                    <  4:16816/149‑2 (MQ=255)                      gccggACCAGACGCGCGATTCCGCCCGCCACATCAACCACCACCTTGACGACCAGAGCCTCGAAACCAGCCTGCCGGGGCATGCGGTGAAGGAGGGGGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTgg                                                                                                                              <  8:258088/146‑1 (MQ=255)                      gCAGGACCTGACTTGCGATTCCGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTg                                                                                                                               <  5:133531/148‑1 (MQ=255)                       cAGGACCTGACTTGCGATTCCGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTggg                                                                                                                             >  8:99506/1‑149 (MQ=255)                       cAGGACCTGACTTGCGATTCCGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTggg                                                                                                                             >  5:358476/1‑149 (MQ=255)                       cAGGACCTGACTTGCGATTCCGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTgg                                                                                                                              <  4:365715/148‑1 (MQ=255)                       cAGGACCTGACTTGCGATTCCGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTgg                                                                                                                              <  3:226855/148‑1 (MQ=255)                            ccTGACTTGCGATTCCGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACa                                                                                                                        >  7:144176/1‑149 (MQ=255)                            ccTGACTTGCGATTCCGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACa                                                                                                                        <  3:277298/149‑1 (MQ=255)                               gACTTGCGATTCCGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGt                                                                                                                            >  5:343068/1‑142 (MQ=255)                               gACTTGCGATTCCGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGt                                                                                                                            <  6:343068/142‑1 (MQ=255)                                aCTTGCGATTCCGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACacg                                                                                                                      <  2:285387/147‑3 (MQ=255)                                aCTTGCGATTCCGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATg                                                                                                                      <  6:8887/147‑1 (MQ=255)                                  ttGCGATTCCGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGcac                                                                                                                   >  5:145134/1‑148 (MQ=255)                                     cGATTCCGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAATTCCTGGGt                                                                                                                            >  1:151969/1‑136 (MQ=255)                                     cGATTCCGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGt                                                                                                                            <  2:151969/136‑1 (MQ=255)                                          ccGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGtt                                                                                                           <  4:170892/148‑1 (MQ=255)                                           cGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCg                                                                                                         >  5:77235/1‑149 (MQ=255)                                           cGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCg                                                                                                         >  3:92466/1‑149 (MQ=255)                                             aCGGCAAGATCTACCAGTACGTTGAAGAGCAGAGCATCGAAACCATGATGGCAGTGCATGGGGTGAAGGATGGTGATGATTACCTGCTTGGCGTGTTCCTGGTGGGTGCCTACGAGGAAATCCTGGGTGACATGCACAACCTGTTCGGt                                                                                                       <  1:358961/149‑1 (MQ=255)                                                 cAAGATCAACCAGTACGTTGACGAGAAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGaca                                                                                                   <  6:358476/149‑1 (MQ=255)                                                  aaGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGATGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGacac                                                                                                  >  2:352367/1‑149 (MQ=255)                                                    gATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACAcc                                                                                                 >  4:51687/1‑148 (MQ=255)                                                      tCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGAt                                                                                              >  3:87206/1‑149 (MQ=255)                                                            aGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCATGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTg                                                                                        <  8:205052/149‑1 (MQ=255)                                                             gTACGTTGAGGAGCATAGCTTCGAAACCTGCATGCCGGTGCATGCTGTGAAGGAGGGTGTGGATTAGCTGCTGGGCGTGTTGCTGGTGGGTGCCTACCATTTAATCCTTTGTGAAATGCTCATGCTGTTGGGTGACACCGATTCGGTGa                                                                                       <  6:77235/149‑1 (MQ=255)                                                             gTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGa                                                                                       <  3:215693/149‑1 (MQ=255)                                                              tACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGATGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGaa                                                                                      <  7:157075/149‑1 (MQ=255)                                                                cGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAAc                                                                                     <  2:156779/148‑1 (MQ=255)                                                                 gTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACa                                                                                    >  5:57786/1‑148 (MQ=255)                                                                 gTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACAt                                                                                   <  3:51687/149‑1 (MQ=255)                                                                     aCGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTAc                                                                               >  2:361947/1‑149 (MQ=255)                                                                         gcagagcaTCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAg                                                                            >  4:232651/1‑148 (MQ=255)                                                                         gcagagcaTCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAg                                                                            >  4:133587/1‑148 (MQ=255)                                                                             agcaTCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGc                                                                       >  4:132030/1‑149 (MQ=255)                                                                               caTCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGt                                                                                         <  4:162301/129‑1 (MQ=255)                                                                               caTCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGt                                                                                         >  3:162301/1‑129 (MQ=255)                                                                                 tCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCg                                                                                           <  1:359403/125‑1 (MQ=255)                                                                                 tCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCg                                                                                           >  2:359403/1‑125 (MQ=255)                                                                                 tCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGAc                                                                   <  3:132030/149‑1 (MQ=255)                                                                                 tCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGAc                                                                   <  6:57786/149‑1 (MQ=255)                                                                                    aaaCCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCGTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGc                                                                >  4:222394/1‑149 (MQ=255)                                                                                     aaCCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCa                                                               <  6:145134/149‑1 (MQ=255)                                                                                       ccAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGc                                                             <  8:144176/149‑1 (MQ=255)                                                                                           cATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCgtgt                                                         <  7:99506/149‑1 (MQ=255)                                                                                                                     gTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATc                                           <  4:249716/137‑1 (MQ=255)                                                                                                                     gTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATc                                           >  3:249716/1‑137 (MQ=255)                                                                                                                        aGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATCGAAACCCACGacac                             >  2:255480/1‑148 (MQ=255)                                                                                                                          gATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGc                                                                       >  5:115304/1‑104 (MQ=255)                                                                                                                          gATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGc                                                                       <  6:115304/104‑1 (MQ=255)                                                                                                                                ctgctgGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATCGAAACCCACGACACCATCGAAGa                    <  4:92466/149‑1 (MQ=255)                                                                                                                                  gctgGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATCGAAACCCACGACACCATCGAAGACa                  >  7:189022/1‑149 (MQ=255)                                                                                                                                  gctgGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATCGAAACCCACGACACCATCGAAGACa                  >  2:305662/1‑149 (MQ=255)                                                                                                                                   ctgGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATCGAAACCCACGACACCATCGAAGAc                   <  1:361947/147‑1 (MQ=255)                                                                                                                                     gggCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTAGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATCGAAACCCACGACACCATCGAAGACAtgc               >  3:191920/1‑149 (MQ=255)                                                                                                                                     ccaggatttCCTGGTAGGCACCCACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATCGAAACCCACGACACCATCGAAGACAtgc               >  8:303893/8‑149 (MQ=255)                                                                                                                                      ggCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATCGAAACCCACGACACCATCGAAGACAtgct              <  3:367469/149‑1 (MQ=255)                                                                                                                                      caggatttCCTGGTAGGCACCCACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATCGAAACCCACGACACCATCGAAGACAtgct              >  5:378064/7‑149 (MQ=255)                                                                                                                                            ttCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATCGAAACCCACGACACCATCGAAGACATGCTGCGCTa        >  8:322929/1‑149 (MQ=255)                                                                                                                                               cTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATCGAAACCCACGACACCATCGAAGACATGCTGCGCTACg      <  1:255480/148‑1 (MQ=255)                                                                                                                                                tggtggGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATCGAAACCCACGACACCATCGaa                      >  6:219751/1‑131 (MQ=255)                                                                                                                                                tggtggGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATCGAAACCCACGACACCATCGaa                      <  5:219751/131‑1 (MQ=255)                                                                                                                                                  gtgggtgCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATCGAAACCCACGACACCATCGAAGACATGCTGCGCTACGTGCa  <  4:87206/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GCCGATGCGGCGTGCGGTATTGCAGGACCTGACTTGCGATTCCGACGGCAAGATCAACCAGTACGTTGACGAGCAGAGCATCGAAACCAGCATGCCGGTGCATGCGGTGAAGGAGGGTGAGGATTACCTGCTGGGCGTGTTCCTGGTGGGTGCCTACCAGGAAATCCTGGGTGACATGCACAACCTGTTCGGTGACACCGATTCGGTGAACATCTACCAGAATGCCGACGGCAGCGTGTACCACGCCGGTATCGAAACCCACGACACCATCGAAGACATGCTGCGCTACGTGCA  >  NC_002947/660442‑660735

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