Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,004,148	G→A	8.4%	Q38Q (CAG→CAA) ‡	PP_3533 →	ornithine cyclodeaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,004,148	0	G	A	8.4%	200.1 / 6.2	71	Q38Q (CAG→CAA) ‡	PP_3533	ornithine cyclodeaminase
Reads supporting (aligned to +/- strand):  ref base G (27/38);  new base A (3/3);  total (30/41)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

GCAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACGGTGATGGCC  >  NC_002947/4004006‑4004295                                                                                                                                                |                                                                                                                                                    gCAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGttt                                                                                                                                               <  1:16734/149‑1 (MQ=255) gCAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGttt                                                                                                                                               >  7:150228/1‑149 (MQ=255)  cAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTACTTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGttt                                                                                                                                               <  1:85059/149‑1 (MQ=255)   aaGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGa                                                                                                                                             <  4:293773/149‑1 (MQ=255)     gcgcACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCCTGCCAGCGTTTGaaa                                                                                                                                           <  1:301433/149‑3 (MQ=255)     gcgcACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCCTGCCAGCGTTTGaaa                                                                                                                                           <  2:12881/149‑3 (MQ=255)      cgcACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACaa                                                                                                                                          <  8:150228/149‑1 (MQ=255)      cgcACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACaa                                                                                                                                          <  2:207734/149‑1 (MQ=255)        cacaATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGt                                                                                                                                        >  6:322254/1‑149 (MQ=255)                    gTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTg                                                                                                                            >  1:106050/1‑149 (MQ=255)                       cAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTgcca                                                                                                                         <  2:95596/149‑1 (MQ=255)                        aaGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGAc                                                                                                                                            <  1:307848/129‑1 (MQ=255)                        aaGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGAc                                                                                                                                            >  2:307848/1‑129 (MQ=255)                        aaGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTgccag                                                                                                                        <  8:221609/149‑1 (MQ=255)                        aaGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTgccag                                                                                                                        >  6:156467/1‑149 (MQ=255)                                gACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCg                                                                                                                 >  8:71638/1‑148 (MQ=255)                                   gTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAg                                                                                                             <  6:321128/149‑1 (MQ=255)                                          aTTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGtt                                                                                                                             >  3:12482/1‑126 (MQ=255)                                          aTTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGtt                                                                                                                             <  4:12482/126‑1 (MQ=255)                                          aTTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCg                                                                                                       <  2:295357/148‑1 (MQ=255)                                            tGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTgccagcca                                                                                                                     <  4:65223/132‑1 (MQ=255)                                            tGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTgccagcca                                                                                                                     >  3:65223/1‑132 (MQ=255)                                              aTGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTc                                                                                                                  >  3:215014/1‑133 (MQ=255)                                              aTGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTc                                                                                                                  <  4:215014/133‑1 (MQ=255)                                                   ccAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAAc                                                                                              <  4:187883/148‑1 (MQ=255)                                                   ccAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCCTGCCAGCGTTTg                                                                                                                                              <  5:302900/100‑1 (MQ=255)                                                   ccAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCCTGCCAGCGTTTg                                                                                                                                              >  6:302900/1‑100 (MQ=255)                                                      aCCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGAt                                                                                          >  6:207431/1‑149 (MQ=255)                                                       ccATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATg                                                                                         <  2:106050/149‑1 (MQ=255)                                                       ccATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATg                                                                                         <  5:207431/149‑1 (MQ=255)                                                        cATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGc                                                                                        <  5:156467/149‑1 (MQ=255)                                                        cATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGc                                                                                        >  4:59228/1‑149 (MQ=255)                                                          tGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCg                                                                                      <  5:308269/149‑1 (MQ=255)                                                           gTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCgg                                                                                     >  4:49201/1‑149 (MQ=255)                                                              ggACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCTAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTc                                                                                   <  3:39472/148‑1 (MQ=255)                                                                aCTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGcc                                                                                <  7:71638/149‑1 (MQ=255)                                                                        cATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACaaaa                                                                         <  8:186175/148‑1 (MQ=255)                                                                            aCATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCg                                                                     <  7:185753/148‑1 (MQ=255)                                                                                          ccGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGATAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATa                                                      >  6:217262/1‑149 (MQ=255)                                                                                             tttA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACaaaa                                                                         >  6:233956/1‑127 (MQ=255)                                                                                             tttA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACaaaa                                                                         <  5:233956/127‑1 (MQ=255)                                                                                             tttA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGt                                                   <  3:49201/149‑1 (MQ=255)                                                                                                  ttGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTAGGCCTTCAAATACGTCAAc                                               >  6:224453/1‑149 (MQ=255)                                                                                                     gCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGc                                            >  3:171769/1‑149 (MQ=255)                                                                                                      cGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGc                                            >  4:187048/1‑148 (MQ=255)                                                                                                          cTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCaaa                                                        >  6:279060/1‑132 (MQ=255)                                                                                                          cTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCAccc                                       >  8:169540/1‑149 (MQ=255)                                                                                                          cTTGCCGCTGCCCTGCGGGACGATATCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCaaa                                                        <  5:279060/132‑1 (MQ=255)                                                                                                           ttGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCg                                      >  5:131720/1‑149 (MQ=255)                                                                                                               cGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGcc                                   >  2:238094/1‑148 (MQ=255)                                                                                                                 cTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCaa                                 >  8:298560/1‑148 (MQ=255)                                                                                                                    cccTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATa                               <  3:212827/147‑1 (MQ=255)                                                                                                                         cGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACg                        <  3:59228/149‑1 (MQ=255)                                                                                                                             aCGATTTCAAACGCTGGCAGGCGTTTGATAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAAc                    <  5:217262/149‑1 (MQ=255)                                                                                                                             aCGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAAc                    <  4:171769/149‑1 (MQ=255)                                                                                                                             aCGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGCGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAAc                    <  7:169540/149‑1 (MQ=255)                                                                                                                                 tttCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTg                 <  8:50133/148‑1 (MQ=255)                                                                                                                                 tttCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTg                 >  7:50133/1‑148 (MQ=255)                                                                                                                                 tttCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGc                <  1:238094/149‑1 (MQ=255)                                                                                                                                  gtCAAACGCCTGCCAGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCa               >  5:44534/2‑149 (MQ=255)                                                                                                                                   tCAAACGCTGGCAGGCGTTTGACAATTCCTCGCGCGATGCCAGCCAATAGTAAGTGTGAGTGATCGAACTTATGCCGGTCGCCGAAAAAATCCGCTATGCCTTCAAATACGTCAACGGAAACCAGGCCTATACTGAACTCAAcattca               >  3:264729/1‑143 (MQ=255)                                                                                                                                   tCAAACGCCTGCCAGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCa               >  2:172906/1‑148 (MQ=255)                                                                                                                                    cAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGc                                                                                        >  6:40474/1‑74 (MQ=255)                                                                                                                                    cAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGc                                                                                        <  5:40474/74‑1 (MQ=255)                                                                                                                                       aCGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACgg          <  3:187048/149‑1 (MQ=255)                                                                                                                                             gCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACGGTGAt      <  7:298560/147‑1 (MQ=255)                                                                                                                                             gCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACGGTGATgg    <  5:307561/149‑1 (MQ=255)                                                                                                                                              cAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACTGCCACCCGGCCAATACTGCACGCAACCTGCATACGGTGATgg    <  5:227267/148‑1 (MQ=255)                                                                                                                                              cAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACGGTGATGGc   <  5:224453/149‑1 (MQ=255)                                                                                                                                              cAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTAGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACACGGCCAATACTGCACGCAACCTGCATACGGTGATgga   >  2:292598/1‑148 (MQ=255)                                                                                                                                               aGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACGGTGATGGcc  >  8:214181/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                    GCAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTA‑TTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACGGTGATGGCC  >  NC_002947/4004006‑4004295

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