BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	191,892	T→A	21.7%	L13Q (CTG→CAG)	pyrH →	UMP kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	191,892	0	T	A	21.7%	344.2 / 94.0	171	L13Q (CTG→CAG)	pyrH	UMP kinase
Reads supporting (aligned to +/- strand): ref base T (53/80); new base A (11/26); total (64/107)
Fisher's exact test for biased strand distribution p-value = 3.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

Only 100 of 173 total aligned reads displayed.

TGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGAAATCAAAGAACTGGTTGAACTGGGTATTCAGGTTGGTGTGGTGATTGGTGGGGGTAACCTGTTCCGTG  >  NC_000913/191745‑192041

   gCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTa                                                                                                                                                   <  2:242091/150‑1 (MQ=255)

   gCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTa                                                                                                                                                   <  1:1406992/150‑1 (MQ=255)

            gTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTg                                                                                                                                          <  2:366531/150‑1 (MQ=255)

               aaTTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCg                                                                                                                                       >  1:2311395/1‑150 (MQ=255)

                 ttCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGaa                                                                                                                                     <  2:2020723/150‑1 (MQ=255)

                   cAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGc                                                                                                                                   >  1:1136222/1‑150 (MQ=255)

                   cAGCTAACCCTTGTGGGGCCGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGATACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCAGCTTAAGTTGAGTGGCGAAGc                                                                                                                                   >  2:1246614/1‑150 (MQ=255)

                        aaCCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGc                                                                                                                              >  1:1288507/1‑150 (MQ=255)

                          cccTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGATAGAAACATGGCTACCAATGCAAAACCCGTCTATAAAAGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAg                                                                                                                            >  2:2374995/1‑150 (MQ=255)

                          cccTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAg                                                                                                                            >  2:1822882/1‑150 (MQ=255)

                                tgGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAACCATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACt                                                                                                                      >  1:1586180/1‑150 (MQ=255)

                                   ggCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGaa                                                                                                                   <  1:2416730/150‑1 (MQ=255)

                                    gCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAg                                                                                                                  >  2:7339/1‑150 (MQ=255)

                                    gCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAg                                                                                                                  <  2:2152604/150‑1 (MQ=255)

                                     cTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAg‑g                                                                                                                <  2:1751757/150‑1 (MQ=255)

                                        cgcTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCtt                                                                                                             >  2:490514/1‑150 (MQ=255)

                                         gcTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTc                                                                                                            >  1:639374/1‑150 (MQ=255)

                                         gcTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCAGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTc                                                                                                            <  2:1629481/150‑1 (MQ=255)

                                          cTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCAGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCg                                                                                                           <  1:116517/150‑1 (MQ=255)

                                          cTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGAAAAACCAGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCg                                                                                                           >  1:1630973/1‑150 (MQ=255)

                                           tGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCgg                                                                                                          <  2:6473/150‑1 (MQ=255)

                                           tGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCgg                                                                                                          <  2:2229860/150‑1 (MQ=255)

                                           tGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCgg                                                                                                          <  2:1451499/150‑1 (MQ=255)

                                           tGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGC‑CTTCAGTGCCCTGCAGAGCTTCgc                                                                                                          <  2:1403311/150‑2 (MQ=255)

                                            gAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGC‑CTTCAGTGCCCTGCAGAGCTTCgcc                                                                                                         <  2:9798/150‑3 (MQ=255)

                                             aaaaGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTa                                                                                                        <  2:1019466/150‑1 (MQ=255)

                                             aaaaGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGC‑CTTCAGTGCCCTGCAGAGCTTCgcca                                                                                                        <  1:2350609/150‑4 (MQ=255)

                                              aaaGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTAt                                                                                                       >  2:281925/1‑150 (MQ=255)

                                              aaaGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGC‑CTTCAGTGCCCTGCAGAGCTTCgccac                                                                                                       <  1:255581/150‑5 (MQ=255)

                                               aaGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTAtt                                                                                                      <  1:2009857/150‑1 (MQ=255)

                                               aaGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTAtt                                                                                                      <  1:1792540/150‑1 (MQ=255)

                                               aaGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCAGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTAtt                                                                                                      <  1:139177/150‑1 (MQ=255)

                                               aaGCGACGTACAATGTCGCAAGTATTAATTAATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTATGCTTAAGTTGAGTGGCGAAGCTCTCCAGGGCACTGAAG‑GCTTCGGTAtt                                                                                                      >  2:590692/1‑150 (MQ=255)

                                                 gCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGa                                                                                                    <  1:895004/150‑1 (MQ=255)

                                                   gACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATg                                                                                                  >  2:671491/1‑150 (MQ=255)

                                                   gACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATg                                                                                                  <  1:23581/150‑1 (MQ=255)

                                                      gTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCaa                                                                                               <  1:2415372/150‑1 (MQ=255)

                                                       tACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAg                                                                                              >  1:860335/1‑150 (MQ=255)

                                                       tACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAg                                                                                              >  1:250769/1‑150 (MQ=255)

                                                       tACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAg                                                                                              >  1:1527107/1‑150 (MQ=255)

                                                        aCAATGTCGCTAGTATTAATTCATTTCAATCGTTGAGAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTGTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGc                                                                                             <  1:793426/150‑1 (MQ=255)

                                                        aCAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGc                                                                                             <  2:232191/150‑1 (MQ=255)

                                                        aCAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCCGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGc                                                                                             <  2:1373669/150‑1 (MQ=255)

                                                        aCAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCAGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGc                                                                                             <  2:1654181/150‑1 (MQ=255)

                                                         cAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCa                                                                                            <  2:2469000/150‑1 (MQ=255)

                                                         aaaTGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCa                                                                                            <  2:460929/149‑1 (MQ=255)

                                                            tGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCAGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATAc                                                                                         <  1:1171884/150‑1 (MQ=255)

                                                              tCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTg                                                                                       <  1:2441870/150‑1 (MQ=255)

                                                               cGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTgg                                                                                      >  2:41063/1‑150 (MQ=255)

                                                                gCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGa                                                                                     <  2:2286051/150‑1 (MQ=255)

                                                                  tAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATc                                                                                   <  2:860335/150‑1 (MQ=255)

                                                                    gTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCAGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGt                                                                                 <  2:1367316/150‑1 (MQ=255)

                                                                         aaTTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCAGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGc                                                                            >  1:1691697/1‑150 (MQ=255)

                                                                          aTTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCAGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCt                                                                           >  2:1616462/1‑150 (MQ=255)

                                                                            tCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCAGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCa                                                                         <  1:464419/150‑1 (MQ=255)

                                                                             cATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAg                                                                        <  1:720660/150‑1 (MQ=255)

                                                                              aTTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAgg                                                                       <  1:671491/150‑1 (MQ=255)

                                                                                  cAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGaaaa                                                                   >  1:1559312/1‑149 (MQ=255)

                                                                                      cGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGAAATCaaa                                                               <  1:550172/150‑1 (MQ=255)

                                                                                      cGTTGACAGTCTAAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGAAATCaaa                                                               >  1:1045712/1‑150 (MQ=255)

                                                                                        ttGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGAAATCAAAGa                                                             <  1:1061194/150‑1 (MQ=255)

                                                                                         tGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCAGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGAAATCAAAGaa                                                            >  1:1138553/1‑150 (MQ=255)

                                                                                                ctcAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGAAATCAAAGAACTGGTTg                                                     <  1:939173/150‑1 (MQ=255)

                                                                                                    ggAAAGATACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCAGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGAAATCAAAGAACTGGTTGAACt                                                 <  1:1246614/150‑1 (MQ=255)

                                                                                                    ggAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGC‑CTTCAGTGCCCTGCAGAGCTTCgcc                                                                                                         <  2:368500/94‑3 (MQ=25)

                                                                                                    ggAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGC‑CTTCAGTGCCCTGCAGAGCTTCgcc                                                                                                         >  1:368500/1‑92 (MQ=25)

                                                                                                       aagaaaCATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCAGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGAAATCAAAGAACTGGTTGAACTggg                                              <  2:562240/150‑1 (MQ=255)

                                                                                                        agaaaCATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGAAATCAAAGAACTGGTTGAACTGGGt                                             <  1:1271830/150‑1 (MQ=255)

                                                                                                           aaCATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGAAATCAAAGAACTGGTTGAACTGGGTAtt                                          <  2:1483519/150‑1 (MQ=255)

                                                                                                            aCATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGAAATCAAAGAACTGGTTGAACTGGGTATTc                                         <  1:2384599/150‑1 (MQ=255)

                                                                                                             cATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGAAATCAAAGAACTGGTTGAACTGGGTATTCa                                        >  1:1263738/1‑150 (MQ=255)

                                                                                                             cATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGAAATCAAAGAACTGGTTGAACTGGGTATTCa                                        <  1:281925/150‑1 (MQ=255)

                                                                                                             cATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGAAATCAAAGAACTGGTTGAACTGGGTATTCa                                        <  2:915587/150‑1 (MQ=255)

                                                                                                                 gCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGAAATCAAAGAACTGGTTGAACTGGGTATTCAGGtt                                    <  1:1649452/150‑1 (MQ=255)

                                                                                                                  cTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGAGGCGAAGCTCTGCAAGGCACTGAAG‑GCTTCGGTAATGATGCAAGCATACTGGATAGTATCGCTCAGGAAATCAAAAAACTGCTTGAACTGGGTATTAAGGttt                                   >  2:1644939/1‑149 (MQ=255)

                                                                                                                   tACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGAAATCAAAGAACTGGTTGAACTGGGTATTCAGGTtgg                                  <  2:250769/150‑1 (MQ=255)

                                                                                                                   tACCAATGCAAAACCCGTCTATAAACGCATTCAGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGAAATCAAAGAACTGGTTGAACTGGGTATTCAGGTtgg                                  <  1:1240322/150‑1 (MQ=255)

                                                                                                                    aCCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGAAATc                                                                  <  2:16821/117‑1 (MQ=255)

                                                                                                                     ccAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGAAATCAAAGAACTGGTTGAACTGGGTATTCAGGTtggtg                                >  2:889775/1‑150 (MQ=255)

                                                                                                                     ccAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGAAATCAAAGAACTGGTTGAACTGGGTATTCAGGTtggtg                                <  1:1324264/150‑1 (MQ=255)

                                                                                                                     ccAATGCAAAACCCGTCTATAAACGCATTCAGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGAAATCAAAGAACTGGTTGAACTGGGTATTCAGGTtggtg                                <  1:1444796/150‑1 (MQ=255)

                                                                                                                      cAATGCAAAACCCGTCTATAAACGCATTCAGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGAAATCAAAGAACTGGTTGAACTGGGTATTCAGGTtggtgt                               >  1:1684256/1‑150 (MQ=255)

                                                                                                                       aaTGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGAAATCAAAGAACTGGTTGAACTGGGTATTCAGGTtggtgtg                              <  1:185330/150‑1 (MQ=255)

                                                                                                                        aTGCAAATCCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGTCGAAGCTCTGCAGTGCACTTAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGAAATCAAAGAACTGGTTGAACTGGGTATTCAGGTtggtgtgg                             <  2:2208004/150‑1 (MQ=255)

                                                                                                                        aTGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGAAATCAAAGAACTGGTTGAACTGGGTATTCAGGTtggtgtgg                             <  1:69035/150‑1 (MQ=255)

                                                                                                                        aTGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGAAATCAAAGAACTGGTTGAACTGGGTATTCAGGTtggtgtgg                             >  1:755732/1‑150 (MQ=255)

                                                                                                                         tGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGAAATCAAAGAACTGGTTGAACTGGGTATTCAGGTtggtgtggt                            >  2:1205371/1‑150 (MQ=255)

                                                                                                                           cAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGAAATCAAAGAACTGGTTGAACTGGGTATTCAGGTTGGTGTGGTGa                          >  2:754397/1‑150 (MQ=255)

                                                                                                                           cAAAACACGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGAAATCAAAGAACTGGTTGAACTGGGTATTCAGGTTGGTGTGGTGa                          >  1:2039305/1‑150 (MQ=255)

                                                                                                                             aaaCCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGAAATCAAAGAACTGGTTGAACTGGGTATTCAGGTTGGTGTGGTGAtt                        <  1:2206894/150‑1 (MQ=255)

                                                                                                                               aCCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGAAATCAAAGAACTGGTTGAACTGGGTATTCAGGTTGGTGTGGTGATtgg                      <  1:796842/150‑1 (MQ=255)

                                                                                                                                  cGTCTATAAACGCATTCAGCTTAAGTTGAGTGGCTAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGAAATCAAAGAACTGGTTGAACTGGGTATTCAGGTTGGTGTGGTGATttgtgg                   <  1:810277/150‑1 (MQ=255)

                                                                                                                                   gTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGAAATCAAAGAACTGGTTGAACTGGGTATTCAGGTTGGTGTGGTGATTGGTggg                  <  2:1288507/150‑1 (MQ=255)

                                                                                                                                    tCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGAAATCAAAGAACTGGTTGAACTGGGTATTCAGGTTGGTGTGGTGATTGGTgggg                 >  2:323634/1‑150 (MQ=255)

                                                                                                                                    tCTATAAACGCATTCAGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGAAATCAAAGAACTGGTTGAACTGGGTATTCAGGTTGGTGTGGTGATTGGTgggg                 <  2:1192160/150‑1 (MQ=255)

                                                                                                                                      tataAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGAAATCAAAGAACTGGTTGAACTGGGTATTCAGGTTGGTGTGGTGATTGGTGGGGGt               >  1:1286558/1‑150 (MQ=255)

                                                                                                                                      tataAACGCATTCAGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGAAATCAAAGAACTGGTTGAACTGGGTATTCAGGTTGGTGTGGTGATTGGTGGGGGt               <  1:604270/150‑1 (MQ=255)

                                                                                                                                         aaaCGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGAAATCAAAGAACTGGTTGAACTGGGTATTCAGGTTGGTGTGGTGATTGGTGGGGGTAAc            <  1:1011343/150‑1 (MQ=255)

                                                                                                                                          aaCGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGAATACTGGATCGTATGGCTCAGGAAATCAAAGAACTGGTTGAACTGGGTATTCAGGTTGGTGTGGTGATTGGTGGGGGTAAcc           >  2:1818066/1‑150 (MQ=255)

                                                                                                                                                ttCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGAAATCAAAGAACTGGTTGAACTGGGTATTCAGGTTGGTGTGGTGATTGGTGGGGGTAACCTGTTcc     <  2:1876332/150‑1 (MQ=255)

TGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAACGCATTCTGCTTAAGTTGAGTGGCGAAGCTCTGCAGGGCACTGAAG‑GCTTCGGTATTGATGCAAGCATACTGGATCGTATGGCTCAGGAAATCAAAGAACTGGTTGAACTGGGTATTCAGGTTGGTGTGGTGATTGGTGGGGGTAACCTGTTCCGTG  >  NC_000913/191745‑192041

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