Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	127,606	C→G	7.5%	intergenic (+19/‑306)	aceF → / → lpd	pyruvate dehydrogenase, dihydrolipoyltransacetylase component E2/dihydrolipoyl dehydrogenase; E3 component of pyruvate and 2‑oxoglutarate dehydrogenases complexes; glycine cleavage system L protein; dihydrolipoamide dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	127,606	0	C	G	7.5%	267.8 / 11.1	80	intergenic (+19/‑306)	aceF/lpd	pyruvate dehydrogenase, dihydrolipoyltransacetylase component E2/dihydrolipoyl dehydrogenase; E3 component of pyruvate and 2‑oxoglutarate dehydrogenases complexes; glycine cleavage system L protein; dihydrolipoamide dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (37/37);  new base G (2/4);  total (39/41)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

TTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTC  >  NC_000913/127459‑127753                                                                                                                                                    |                                                                                                                                                     ttCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGAGCCCGTTTCATTACCATCATTAACAACACGCGGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAAcg                                                                                                                                                     <  4:83707/149‑1 (MQ=255)     tGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAAcggccg                                                                                                                                                 >  4:196517/1‑149 (MQ=255)             cTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTCACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCtttttt                                                                                                                                         <  1:55007/149‑1 (MQ=255)                aTGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCtttttt                                                                                                                                         <  4:272200/146‑1 (MQ=255)                aTGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTg                                                                                                                                     <  5:31287/149‑1 (MQ=255)                  gctgcCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTgg                                                                                                                                    >  4:11457/1‑148 (MQ=255)                        gATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAAtct                                                                                                                              >  1:25613/1‑148 (MQ=255)                            tctctctcCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATg                                                                                                                          <  7:244926/148‑1 (MQ=255)                              tctctcCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAAt                                                                                                                       >  8:292265/1‑149 (MQ=255)                                  tcCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCCGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGt                                                                                                                     <  4:202658/147‑1 (MQ=255)                                    cTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTAtt                                                                                                                  >  8:246348/1‑148 (MQ=255)                                        gACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAgg                                                                                                              <  3:258967/148‑1 (MQ=255)                                        gACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGt                                                                                                             >  6:328726/1‑149 (MQ=255)                                              cgcgTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACACCACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAg                                                                                                       <  4:158319/149‑1 (MQ=255)                                                   gATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGaa                                                                                                   >  7:122394/1‑148 (MQ=255)                                                    aTCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTACAAGAGCCGGCCCAACGGCCGGCGTTTTT‑CTGTTAATCTCATGAATGTATTGAGGTTATTAGCGAATa                                                                                                 <  6:297129/149‑1 (MQ=255)                                                    aTCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATa                                                                                                 <  7:50999/149‑1 (MQ=255)                                                      cGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGa                                                                                               >  8:298487/1‑149 (MQ=255)                                                        aCGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGAc                                                                                              >  7:158732/1‑148 (MQ=255)                                                            aatttttaaacaacaaacggcaaccgatttgtctattcgctaattacctcaatacattcatgagattaccagaaaaaAGCCGGCCGTTGGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGttatta                                                                                                        <  7:162021‑M1/57‑1 (MQ=255)                                                            aatttttaaacaacaaacggcaaccgatttgtctattcgctaattacctcaatacattcatgagattaccagaaaaaAGCCGGCCGTTGGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGttatta                                                                                                        >  8:162021‑M1/78‑134 (MQ=255)                                                                 aTGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGtt                                                                                     <  7:246348/148‑1 (MQ=255)                                                                 aTGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTg                                                                                    >  4:278348/1‑149 (MQ=255)                                                                     tGCCCGTTTCATTAACATCATTAACAACACGCTGTCTGAAATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCAGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGACGt                                                                                >  7:95958/1‑149 (MQ=255)                                                                           tttCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTttgttg                                                                          >  2:89983/1‑149 (MQ=255)                                                                               aTTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCgg                                                                                       <  7:96376/132‑1 (MQ=255)                                                                               aTTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCgg                                                                                       >  8:96376/1‑132 (MQ=255)                                                                               aTTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGttt                                                                       <  3:196517/148‑1 (MQ=255)                                                                                ttACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGtt                                                                                                            <  6:7697/110‑1 (MQ=255)                                                                                  aCCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTaaa                                                                    <  3:212136/148‑1 (MQ=255)                                                                                    catcatTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAAt                                                                 >  3:296931/1‑149 (MQ=255)                                                                                     atcatTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAAtt                                                                >  6:330099/1‑149 (MQ=255)                                                                                     atcatTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAAtt                                                                <  5:328726/149‑1 (MQ=255)                                                                                        atTAACAACACGCTGTCTGACAATCGCCGTCTGGTCATGTAAGTAAAAGAGCCGGCCCCACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGtt                                                             <  4:158725/149‑1 (MQ=255)                                                                                         ttAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGtt                                                             >  1:11233/1‑148 (MQ=255)                                                                                             caacaCGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACaa                                                        >  6:304094/1‑149 (MQ=255)                                                                                              attcgctaataacctcaatacattcatgagattaccagaaaaaAGCCGGCCGTTGGGCCGGCTCTTTTACTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACaa                                                        <  6:101030‑M1/106‑1 (MQ=255)                                                                                              aacaCGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGGCGTTTGTTGTTTAAAAATTGTTTACAAt                                                       >  3:181552/1‑149 (MQ=255)                                                                                                 acGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAAtttt                                                    <  6:89045/149‑1 (MQ=255)                                                                                                 acGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAAtttt                                                    <  2:290188/149‑1 (MQ=255)                                                                                                   gCTGTCTGACATTCGTCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTg                                                   >  7:15608/1‑148 (MQ=255)                                                                                                   gCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGt                                                  <  3:278348/149‑1 (MQ=255)                                                                                                   ctaataacctcaatacattcatgagattaccagaaaaaAGCCGGCCGTTGGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTg                                                   <  8:2984‑M1/110‑1 (MQ=255)                                                                                                         tGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTTACAATTTTGTAAAATa                                            >  5:222/1‑149 (MQ=255)                                                                                                         tGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTCGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATa                                            >  1:51338/1‑149 (MQ=255)                                                                                                          gACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATAc                                           >  1:270716/1‑149 (MQ=255)                                                                                                           aCATTCGCCGTCTGGTTATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATTTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATAc                                           <  3:107089/148‑1 (MQ=255)                                                                                                             aTTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTa                                                 >  4:79523/1‑140 (MQ=255)                                                                                                             aTTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTa                                                 <  3:79523/140‑1 (MQ=255)                                                                                                             aTTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGa                                        >  7:152620/1‑149 (MQ=255)                                                                                                              ttCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGAc                                       <  7:36978/149‑1 (MQ=255)                                                                                                              ttCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGAc                                       >  3:243031/1‑149 (MQ=255)                                                                                                                cGCCGTCTGGCGATGTAAGTAAACGAGCCGGCCCACCGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGCTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGAc                                       <  4:181552/147‑1 (MQ=255)                                                                                                                 gCCGTCTGTTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGa                                    <  7:292265/149‑1 (MQ=255)                                                                                                                 gCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATAACGACGGa                                    >  4:134158/1‑149 (MQ=255)                                                                                                                 gCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGACGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACgg                                     <  8:95958/148‑1 (MQ=255)                                                                                                                   cGTCTTGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATa                                  <  8:229412/149‑1 (MQ=255)                                                                                                                   cGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATa                                  <  5:297543/149‑1 (MQ=255)                                                                                                                     tCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGa                                <  4:243031/149‑1 (MQ=255)                                                                                                                         gTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGaa                               <  1:332709/146‑1 (MQ=255)                                                                                                                         gTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGaa                               <  8:158732/146‑1 (MQ=255)                                                                                                                              gtaagtaaAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGa                                    >  1:16583/1‑136 (MQ=255)                                                                                                                              gtaagtaaAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGa                                    <  2:16583/136‑1 (MQ=255)                                                                                                                              gtaagtaaAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCgg                       >  3:109517/1‑149 (MQ=255)                                                                                                                                 cagaaaaaAGCCGGCCGTTGGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGaa                               >  6:65506‑M1/9‑138 (MQ=255)                                                                                                                                 cagaaaaaAGCCGGCCGTTGGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGaa                               <  5:65506‑M1/130‑1 (MQ=255)                                                                                                                                 agtaaAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTGTGTTGTTTAAAAATTGTTATCAATTTTGTAAAATACCGACGGATAGAACGACCCggtgg                    >  5:16977/1‑149 (MQ=255)                                                                                                                                  gtaaAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCggtggt                   >  1:270779/1‑149 (MQ=255)                                                                                                                                  gtaaAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCggtggt                   >  3:159138/1‑149 (MQ=255)                                                                                                                                   taaAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCggtggtg                  >  7:102620/1‑149 (MQ=255)                                                                                                                                   taaAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCggtggtg                  >  7:59028/1‑149 (MQ=255)                                                                                                                                   taaAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCggtggtg                  <  8:152620/149‑1 (MQ=255)                                                                                                                                    aaaaGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCggtggtg                  <  2:270716/148‑1 (MQ=255)                                                                                                                                     aaaGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTg                                                               >  2:134040/1‑102 (MQ=255)                                                                                                                                     aaaGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTg                                                               <  1:134040/102‑1 (MQ=255)                                                                                                                                             GGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTAt        >  2:25613‑M1/1‑149 (MQ=255)                                                                                                                                              GCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTAt        <  7:14248‑M1/148‑1 (MQ=255)                                                                                                                                               CCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTa      >  1:158180‑M1/1‑149 (MQ=255)                                                                                                                                               CCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTa      >  4:296931‑M1/1‑149 (MQ=255)                                                                                                                                                CCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTAc     <  3:32282‑M1/149‑1 (MQ=255)                                                                                                                                                  aaCGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCAGGTGGTGGTTAGGGTATTACt    >  1:232404/1‑148 (MQ=255)                                                                                                                                                   aCGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTc  <  7:27812/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     TTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTT‑CTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTC  >  NC_000913/127459‑127753

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