Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	1,105,994	T→C	T12T (ACT→ACC)	pdhD →	dihydrolipoyl dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,105,994	0	T	C	98.6%	235.8 / ‑5.2	70	T12T (ACT→ACC)	pdhD	dihydrolipoyl dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (34/34);  total (35/35)
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.58e-01

TTGACCACAGACAAATTGATGGTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACTATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGTAGGATGTATTCCTTCAAAAGCATT  >  CP000730/1105858‑1106131                                                                                                                                         |                                                                                                                                          ttGACCACAGACAAATTGATGGTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAg                                                                                                                                       >  1:509469/1‑141 (MQ=255) ttGACCACAGACAAATTGATGGTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAg                                                                                                                                       <  1:103222/141‑1 (MQ=255)  tGACCACAGACAAATTGATGGTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGt                                                                                                                                      <  2:434162/141‑1 (MQ=255)  tGACCACAGACAAATTGATGGTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTCATAATAAATACCAGAAATTGTACTAGGAGAGGGGAATAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGt                                                                                                                                      >  1:463655/1‑141 (MQ=255)   gACCACAGACAAATTGATGGTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTa                                                                                                                                     >  1:958870/1‑141 (MQ=255)    aCCACAGACAAATTGATGGTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTaa                                                                                                                                    <  1:32351/141‑1 (MQ=255)          gACAAATTGATGGTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAg                                                                                                                              >  2:861877/1‑141 (MQ=255)            cggcatcgttgTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATTGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCa                                                                                                                            <  2:958750/131‑1 (MQ=255)             aaaTTGATGGTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAg                                                                                                                           <  2:996700/141‑1 (MQ=255)             aaaTTGATAGCGCAACTGGCGAAACTGCAATGAAACACAGTATACCTTCATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAg                                                                                                                           <  1:970530/141‑1 (MQ=255)              aaTTGATGGTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAgg                                                                                                                          >  1:505603/1‑141 (MQ=255)                ttGATGGTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTc                                                                                                                        >  1:267939/1‑141 (MQ=255)                 aaatcagtttcatGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTcc                                                                                                                       <  2:63522/129‑1 (MQ=255)                       tGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCtggtgg                                                                                                                 <  1:733680/141‑1 (MQ=255)                        gccaCTGGCCAAAATGCAATGAATCACATTAAACGATAATTCAATCATCCCTGATTATTAAATACTGAGGGTTAAATTATCGTAGTTTGAGATTGCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGa                                                                                                                <  1:358586/138‑1 (MQ=255)                        gCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCCGACAAATTATTAATGGAGGGGTAAAACATGATAGTTGGAGATTTACCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGa                                                                                                                >  1:235969/1‑141 (MQ=255)                        gCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGa                                                                                                                >  2:122718/1‑141 (MQ=255)                            cTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACg                                                                                                            >  2:988613/1‑141 (MQ=255)                                  aaaaTGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTgcag                                                                                                      >  1:417424/1‑141 (MQ=255)                                  aaaaTGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTgcag                                                                                                      <  1:644013/141‑1 (MQ=255)                                   aaaTGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTgcagc                                                                                                     <  2:87944/141‑1 (MQ=255)                                   aaaTGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTgcagc                                                                                                     <  2:121797/141‑1 (MQ=255)                                        cAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACTATAGTAATCGGAGCAGATCCTGGTGGATACGTTGCAGCAATTc                                                                                                >  2:814856/1‑141 (MQ=255)                                        cAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTc                                                                                                >  1:315494/1‑141 (MQ=255)                                          aTGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGt                                                                                              >  2:832679/1‑141 (MQ=255)                                          aTGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGt                                                                                              >  2:575177/1‑141 (MQ=255)                                           tGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCCGCAATTCTTg                                                                                             >  2:557751/1‑141 (MQ=255)                                           tGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTg                                                                                             >  1:500657/1‑141 (MQ=255)                                            gAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGc                                                                                            <  2:638132/141‑1 (MQ=255)                                            gAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGc                                                                                            <  1:1018909/141‑1 (MQ=255)                                            gAATCACATTAAACGGGTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTACTGGTGGATACGTTGCAGCAATTCGTGc                                                                                            <  2:638370/141‑1 (MQ=255)                                               tCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGc                                                                                         <  1:68797/141‑1 (MQ=255)                                               tCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGc                                                                                         <  1:343505/141‑1 (MQ=255)                                                  caTTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCa                                                                                      >  2:600418/1‑141 (MQ=255)                                                            ttattaAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACaa                                                                            >  1:123973/1‑141 (MQ=255)                                                               ttaAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACaaaaa                                                                         >  1:1016432/1‑141 (MQ=255)                                                               ttaAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACaaaaa                                                                         >  1:346312/1‑141 (MQ=255)                                                                  aataatCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTa                                                                      <  2:737199/141‑1 (MQ=255)                                                                  aataatCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTa                                                                      <  1:1022617/141‑1 (MQ=255)                                                                   ataatCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTaa                                                                     <  2:169774/141‑1 (MQ=255)                                                                    taatCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAAc                                                                    >  1:720460/1‑141 (MQ=255)                                                                    taatCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAAc                                                                    <  2:912086/141‑1 (MQ=255)                                                                    taatCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAAc                                                                    >  1:696403/1‑141 (MQ=255)                                                                    aaatCCAGAATTATTATTAATGGAGAGGTAAAACATGGTAGTTGGAGATTTACCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAAc                                                                    <  2:240263/140‑1 (MQ=255)                                                                        ccAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTCGGACAAAAAGTAACAATc                                                                >  2:868173/1‑141 (MQ=255)                                                                         cAGAATTATTATTAATGGAGTGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCg                                                               <  1:941185/141‑1 (MQ=255)                                                                           gAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGtt                                                             >  1:288568/1‑141 (MQ=255)                                                                                attattaATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGaa                                                        >  2:157645/1‑141 (MQ=255)                                                                                  aattaATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAg                                                      <  2:53416/140‑1 (MQ=255)                                                                                     taATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTa                                                   >  2:897019/1‑141 (MQ=255)                                                                                            ggggTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTtgg                                            >  2:572718/1‑141 (MQ=255)                                                                                                  aaaCATGGTAGTTGGAGATTTCGCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCATTTCGTGCAGATAAAAAAGGAAAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGgtgt                                      <  2:1009928/141‑1 (MQ=255)                                                                                                   aaCATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTAAAAAAAAAAAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGtt                                     <  2:314268/141‑1 (MQ=255)                                                                                                      aTGGTAGTTGGAGATTTCCCAATTGAAACGGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAAACTTGGTGGTGTTTGc                                  >  1:612000/1‑141 (MQ=255)                                                                                                        ggTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGAAAAAAACTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCtt                                <  2:884807/141‑1 (MQ=255)                                                                                                              tGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGt                          >  1:81602/1‑141 (MQ=255)                                                                                                                   aTTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGTAGGAt                     <  2:500657/141‑1 (MQ=255)                                                                                                                       cccAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGTAGGATGTAt                 >  2:452519/1‑141 (MQ=255)                                                                                                                        ccAATTGAAGCAGATACGAGAGTAATCGGAGCAGGTCCTGGTGGGTACGTTCCAGCAATTCGTGCAGCTCAATTAGGACAGAAAGTAACAAGCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGTAGGATGTAtt                <  2:492021/141‑1 (MQ=255)                                                                                                                         cAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGTAGGATGTAt                 <  1:596856/139‑1 (MQ=255)                                                                                                                         cAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGTAGGATGTAt                 >  2:596856/1‑139 (MQ=255)                                                                                                                          aaTTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGTAGGATGTATTcc              <  1:545779/141‑1 (MQ=255)                                                                                                                            ttGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGTAGGATGTATTCCtt            >  1:70431/1‑141 (MQ=255)                                                                                                                            ttGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGTAGGATGTATTCCtt            <  1:709994/141‑1 (MQ=255)                                                                                                                              gAATCAGGTACCATACTAAAGGGAGCGGGACCTGGTGGATACGTTGCAGCAATTCGTGCAGATCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGTAGGATGTATTCCTTCa          <  2:463655/141‑1 (MQ=255)                                                                                                                                 aCAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGtt                                     >  1:177055/1‑111 (MQ=255)                                                                                                                                 aCAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGtt                                     <  2:177055/111‑1 (MQ=255)                                                                                                                                 aCAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGTAGGATGTATTCCTTCaaaa       >  2:838899/1‑141 (MQ=255)                                                                                                                                 aCAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGTAGGATGTATTCCTTCaaaa       >  2:838910/1‑141 (MQ=255)                                                                                                                                      tACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGTAGGATGTATTCCTTCAAAAGCAtt  <  2:123973/141‑1 (MQ=255)                                                                                                                                         |                                                                                                                                          TTGACCACAGACAAATTGATGGTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACTATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGTAGGATGTATTCCTTCAAAAGCATT  >  CP000730/1105858‑1106131

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